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I. INTRODUCTION 


Antennas have been a part of our lives since most of us can remember. Whether the 
simple car whip antenna, the old Ultra-High Frequency (UHF) “rabbit ears,” or the 
television antenna that seemed to blow off the roof every winter, we have come to accept 
them as a normal part of life. Within the military, our daily use of antennas is a constant 
which causes us to take for granted the process of optimizing an antenna design to the 
specific task/weapon system. Truly, we are accustomed to seeing the effects of an antenna, 
and never the causes for what we see on the screen or hear through a speaker. It is only 
when the effects are not what we desire that we may actually think of the antenna as an 
integral part of the system. 

As technology continues to move forward, though, there are increasing applications 
of very specialized antennas for weapon systems and communications, alike. This work 
centers around just such an application. Specifically, it stems from a requirement to design 
an antenna that is capable of receiving and accurately reproducing a short-duration impulse 
of electromagnetic energy, or an impulse receiving antenna. From signal analysis theory, it 
can be shown that a single pulse of electromagnetic energy that is of extremely short time 
duration (an impulse) will infer a frequency characteristic of very wide bandwidth [1]. 
Specifically, if one were able to create a pulse of infinitely short duration, the resulting 
frequency bandwidth of that impulsive signal would be infinitely wide. In reality, this is 
impossible; however, that inverse relationship between signals in the time- and frequency- 
domain is still valid. For that reason, alone, if we look to design and construct an antenna 
which is able to transmit or receive such a short duration, impulsive signal, then we must 
take into consideration the requirement for the antenna to operate flawlessly across that 
entire ultra-wide bandwidth. 

This research will investigate the theory behind ultra-wideband (UWB) antennas 
and develop a prototype design of an UWB antenna that is capable of receiving and 


replicating an impulse at near-perfect fidelity over the frequency range of 10 MHz to 10 


GHz (1000:1 bandwidth ratio). While there exists a number of antennas which are 
advertised to be ultra-wideband in nature, the reality is that most do not have the linear 
phase characteristic necessary for replication of time-domain signals. Specific examples 
(and their maximum bandwidths) are spiral antennas (50:1), quadraphase antennas (36:1), 
and log-periodic antennas (40:1). 

In addition, while there are certainly numerous "tools" available to the antenna 
designer, the unique requirement to preserve the fidelity of the signal in the time domain 
poses the question: How does one best acknowledge an antenna's performance in the time 
domain without actually analyzing it there? To that end, the time-domain representation of 
an antenna and its near electromagnetic fields (near fields) is, to say the least, a dynamic 
environment. But the requirement to measure performance in the time domain is a 
necessity in determining the ability of such an antenna to meet its design specifications. 

As such, a specific goal within the design process, and the primary focus of this 
work, is the creation, testing, and application of software that would aid in the design 
process through a visual display of the dynamic near fields in the time domain. Through 
this, an individual of even minimal background in antenna design could then be afforded 
the opportunity to “see and learn," to observe the changes in the near fields, and 
heuristically understand the physics of what is occurring. 

To begin, we will first discuss the theory behind UWB antennas, explore the 
possible designs which would fulfill our requirement, and choose a prototype design for 
analysis. Following that, this work will explain in detail the combination of the Analysis of 
Wire Antennas and Scatterers (AWAS) and Near-Field Visualization (NFLDVIS. M) 
software packages for use in the visualization and animation of the antenna's near fields. 
From there, this work will trace the validation/verification of the software as an accurate 
method for displaying the near fields, and then conclude with the actual analysis of a 


prototvpe UWB antenna design. 


Il. ULTRA-WIDEBAND ANTENNA THEORY 


As stated earlier, our search for an antenna that is able to accurately reproduce the 
time-domain waveform of a received impulse signal requires us to look for an antenna 
that is able to operate over a very wide range of frequencies. The shorter the time 
duration of the impulse waveform, the wider the frequency range over which our antenna 
must operate. Our discussion, now, turns to the definition of wideband and the antennas 


that truly meet this definition. 


A. BACKGROUND 


A broadband, or wideband, antenna is generally defined as one whose impedance 
and pattern “do not change significantly over about an octave (f/f, =2 )” [2], where f, 
and fọ are the lower and upper limits on the bandwidth, or frequency range, respectively. 
Using this definition, then, it is easy to see that our desired bandwidth ratio 
(fy/f, =10GHz/10 MHz =1000) implies a need for an ultra-wideband antenna. In 
general, a typical wideband antenna must be a structure that exhibits a smoothness 
throughout its dimensions and does not tend to have abrupt physical changes, even at the 
boundaries. This structure will then produce constant, smooth field amplitudes and input 
impedances that either remain unchanged or change linearly (smoothly) with frequency 
[2]. | 

Long-wire antennas are well known as having excellent broadband properties, 
especially in the High Frequency (2 MHz to 30 MHz) range. As the name implies, these 
antennas are simply long wires supported by towers at either end. They may be operated 
as either resonant antennas, where the current on the wire will be in the form of a 
standing wave, or non-resonant (traveling wave) antennas. The physical difference 
between the two is the introduction of a matched load at the end of the non-resonant long- 
wire antenna which creates a traveling current wave in much the same manner as a 


transmission line terminated at a matched load. Also, the traveling-wave antennas are 


u) 


more directive, lacking the presence of a major lobe in the reverse direction of antenna 
orientation. As an example, a long-wire antenna oriented along the positive z-axis will 
radiate main lobes in both the positive and negative z-direction if the antenna is a 
resonant one, and only in the positive direction if it is non-resonant. 

The traveling-wave antennas have an advantage over their resonant cousins in that 
their input impedance is predominantly resistive (real), just as the impedance of a pure 
traveling wave on a low-loss transmission line is equal to the characteristic, real 
impedance of that line. Because of this, and the fact that the traveling-wave antenna’s 
input impedance is also relatively independent of frequency, these non-resonant antennas 
will result in both an improved radiation pattern and a wider bandwidth as compared to 
the equivalent resonant structure. [2] 

One drawback to both types of long-wire antennas is their size. An antenna 
designed for a 10 MHz operating frequency would be at least 15 meters long (A,/2), and 
can easily exceed six wavelengths in length. As a result, these antennas can in no way be 
considered even remotely portable, or trainable in any but one direction. A second 
drawback is that except for those frequencies where the antenna’s main lobe is nearly 
parallel to the antenna structure, there will always be a portion of the radiation reflected 
off the ground. This reflected field is then out of phase with the direct radiation. For a 
receiver, the reflected field energy will have a phase delay in comparison with the direct 
field. Although this does not create much of a problem for normal continuous wave 
(CW) applications, we will see that a phase delay will have significant impact upon an 
impulse receiving antenna. Lastly, although the input impedance is frequency 
independent, it is very much dependent upon antenna length, and so the longer the 
antenna, the greater input impedance (resistance) required. This serves to severely 
restrict the upper frequency limit, fj, to less than 1 GHz. 

Other broadband antennas, such as helical antennas, biconical antennas, sleeve 
antennas and spiral antennas can generate larger bandwidth ratios than the long-wire 


antennas; however, they tend to suffer problems of scale when the bandwidth ratio begins 


to reach 20:1. For example, the spiral antenna, which is typically referred to as a 
frequency independent antenna due to the fact its radiation pattern and impedance remain 
constant over a bandwidth ratio of 10:1 or higher, is "designed to minimize finite lengths 
and maximize angular dependence" [2]. However, its outer radius, R, is still equal to 
\,/4, or 7.5 meters for f,=10 MHz. While providing a better bandwidth ratio than the 
long-wire antenna, at over 40 feet in diameter, it is even less portable than the long-wire 
antenna, and more difficult to construct. In addition, the phase response of the spiral 
antenna is highly nonlinear, resulting in dispersion of the time-domain pulse [4]. This, 
alone, makes the spiral antenna unsuitable as an impulse receiving antenna. 

Another antenna with general wideband applications is the log-periodic antenna, 
which is seen quite often for television reception. These antennas tend to maintain simple 
geometrical shapes or involve merely straight wire segments, and thus are more easy to 
construct than the spiral antenna. The underlying concept in their construction is that the 
log-periodic antenna is scaled into itself periodically as the frequency (or wavelength) 
changes. This stems from the fact that an antenna which operates well at a wavelength A, 
will operate just as well at a wavelength A, ıf the antenna’s dimensions are changed by a 
factor of A,/A,. It can be shown that every dimension from one element to the next in the 


log-periodic antenna is related by a factor z, which 1s called the log period. This can be 


determined simply from the fact that for a log-periodic antenna, log(f/fi)-log(7) [2], 
where f,2fj, and f, is the next lower design frequency where optimal performance is 
required. | 

This antenna, while providing greater bandwith flexibility, still can be hindered by 
size as f, approaches 10 MHz. In addition, the highest frequencies of a broadband pulse 
applied to a log-periodic antenna will occur first, and the lower frequencies will occur last 
when the antenna response is viewed in the time domain. This dispersion in the time- 
domain is driven by the poor phase response of the antenna. Like the highly dispersive 


spiral antenna, the nonlinear frequency response of the log-periodic antenna has little 


impact on CW operations, but will be detrimental to the antenna's operation as an 


impulse receiver. 


B. IMPULSE RECEIVING ANTENNAS 


Traveling-wave wire antennas, and helical, spiral, and log-periodic antennas are 
but a few of the antennas which support broadband CW applications. However, when the 
signal changes from a CW one to an impulse, we find the dynamics of the antenna design 
process can change, too. Specifically, there tends to arise additional parameters that must 
be taken into consideration when an antenna is no longer being used for CW operations. 
As previously discussed, the requirement driving this research is to design an UWB 
antenna to be used for receiving time-domain signals with near-perfect fidelity. Thus, it 
is extremely critical that the shape of the impulse waveform be preserved throughout the 
entire length of the antenna. Specifically, then, this antenna will have an additional 
requirement to maintain a linear phase change across the entire frequency range. This is 
because if any group of frequencies is either favored or rejected over the others, the time- 
domain waveform as seen by the antenna will be distorted in comparison to the original 
incident waveform [3]. From [1], it can be shown that for any signal, a delay in the time 
domain is caused by a frequency shift in the frequency domain. In the case of a single 
Gaussian impulse, a frequency shift in one group of frequencies compared to the rest will 
cause a time delay in those frequency components, only. This will then cause dispersion 
in the output waveform. If this frequency shift (phase change) is linear across all 
frequencies, however, the entire signal will be delayed in time, and the precise shape of 
the original waveform (signal) will be left intact. 

Additionally, since our antenna is intended to only act as a receiver, the effective 
length of the antenna will prove to be critical. Effective length is determined by 


V (0) 
E' (0) ' 





Log (0) = (2.1) 


where o — 2n/ is the frequency in radians, V,,(@) is the open circuit voltage produced at 
the terminals of an antenna with an electric field E(c) incident upon it [2],[4]. This is a 
key design parameter for an UWB antenna because in order to ensure the output voltage 
waveform is a near-perfect replication of the original incident pulse, the effective length 
must remain constant through the entire frequency range. The relationship between gain, 


G, and L,« is given by [4] 


Crs 


a Laf n (2.2) 


A 
where 7, is the intrinsic impedance of free space, f is the frequency in Hz, c is the speed 
of light, and R, is the radiation resistance of the antenna. Through this relationship, it is 
clear that an UWB receiving antenna requiring a constant effective length across its entire 
bandwidth will not have a constant gain when the radiation resistance is held constant, 
too. Rather, the gain will be proportional to the square of the frequency, which is 
different from standard CW antennas and even impulse transmitting antennas. 

All of the previously discussed antennas can severely distort the received 
waveform, predominantly due to the nonlinear phase change mentioned above. Because 
of this and their overall inflexible size constraints, more specialized antennas must be 
found that will give us the bandwidth ratio and level of fidelity we desire. Reference [4] 
is the most comprehensive exploration of all the possible UWB impulse receiving 
antennas which may meet our bandwidth ratio requirement. The remaining portion of 
this section will summarize the pros and cons for each antenna type investigated in [4]. 

Resistively loaded linear antennas are "similar in concept to the standard 
traveling-wave antennas except, instead of a single resistive termination at the end of the 
device, the antenna is continuously loaded along its entire length" [4]. Although not 
highly directional, a 14cm long dipole was found to exhibit little dispersion and achieve a 
relatively constant effective length from 100 kHz to 1 GHz; however, this antenna's 
ability to reach frequencies higher than 1 GHz is hampered by the fact its input 


impedance is over 20022. Due to this, a definite impedance mismatch between the 


antenna and load will exist. To compensate for this and ensure balanced operation, a 
transformer will be required. Since no transformer is currently capable of operating at 
frequencies above 1 GHz, any effort to scale the antenna for use within our desired 
bandwidth would surely be moot. [4] 

An orthogonally polarized quadraphase antenna has the advantages of both low 
dispersion and a low Voltage Standing Wave Ratio (VSWR) over its bandwidth of 500 
MHz to 18 GHz. It is also easy and inexpensive to build. However, it lacks a frequency 
independent effective length, which will most certainly yield distortion of the time- 
domain waveform. In addition, the same structural size dimensions required of the long- 
wire antennas to reach f; as discussed earlier make this antenna impractical, too. [4] 

Modified log-periodic antennas differ from standard log-periodic antennas in that 
their spacing from element to element is changed to produce a more linear phase 
response. This, in turn, provides for less distortion of the time-domain waveform, and 
better overall performance. Unfortunately, though, the size of the elements required to 
receive the lowest frequency components of the impulse signal are still on the order of a 
half-wavelength, so even modified, the log-periodic antenna is not a wise choice for this 
application. [4] 

Short axial length broadband horns and reflector antennas with wideband impulse 
feeds also show promise as impulse receiving antennas and have demonstrated solid 
wideband applications [4]. As with nearly every antenna mentioned previously, though, 
these, too, suffer from size constraints required to receive the low frequency components. 
The fact that they can receive an ultra-wideband pulse effectively is far overshadowed by 
their required size to get the job done. Thus, when we impose any type of size restriction 
or directional/mobility requirement upon our antenna, it is clear that these antennas will 
not make “the cut.” 

The only known antenna type which can be classified as truly ultra-wideband and 
Is both relatively easy to manufacture and somewhat portable is the TEM horn. Named 


after the fact that it allows the fields within its structure to remain in the transverse 


electromagnetic (TEM) mode, the TEM horn combines constant, frequency-independent 
effective length, good directionality, and a simple design to provide near-perfect fidelity 
in the reception of impulsive signals [3]. It has, without question, the best potential for 


meeting the design requirements previously set forth. 


E: TEM HORN ANTENNA 


The TEM horn, as illustrated in Figure 1, consists of two conductive triangular 
plates which gradually slope along an angle of no more than ten degrees, and then are 
attached to two resistively-loaded parallel plates. For the prototype that will be analyzed 
later, the exact dimensions are 72in x 6in x 4in (1.8288m x 0.1524m x 0.1026m). The 
triangular sections extend 36in (0.9144m) from the feed point along the y-axis, and the 
included angle between them is approximately 9.5?. The parallel plates are 36in x 4in 
(0.9144m x 0.1026m), and have resistive loads running parallel to the antenna boresight 
axis (the positive y-direction). These resistive loads increase roughly exponentially with 
distance from the flare-parallel-plate interface. It should be apparent that the TEM horn 


can be modeled, in a very approximate sense, as a parallel plate transmission line. 





Figure 1. Example of a TEM Horn 


To discover the benefits of the TEM horn in UWB antenna design, we first turn 
our attention to the concepts behind its development. As stated earlier, the time-domain 


antenna is constructed to accurately reproduce the time waveform of the incoming 


electromagnetic fields. Unfortunately, most common receiving antennas severely distort 
these incident fields, allowing any number of modes to propagate from a single incident 
wave. These modes, in turn, cause undesirable distortion in the received waveform, 
restricting the operational bandwidth where accurate fidelity can be maintained. The 
TEM horn is designed to minimize these modes, and thus minimize distortion by 
allowing the electromagnetic fields within the structure to remain in the transverse 
electromagnetic (TEM) mode. 

To provide for near-perfect fidelity in the reproduction of a received waveform, it 
is critical that the antenna maintain a constant (or linear) phase response. To accomplish 
this, the TEM horn is constructed so that the currents induced in the conductors by the 
incident fields move in the same direction as the fields, themselves. Because of this, the 
induced currents and incident fields stay in good time synchronization, thereby preserving 
the phase response. 

Contrarily, a long dipole antenna impinged by those same fields has induced 
currents which would move perpendicular to the direction of propagation. This leads to a 
time lag between arrival of the incident fields and the induced currents at the dipole’s 
feed point. Because of this time lag, the resultant waveform (taken as output from the 
feed) would exhibit dispersion, and thus distortion, in comparison to the original incident 
wave. 

An additional phase shift can result from a time delay caused by the difference in 
path length between the conductors along the centerline axis of the TEM horn and those 
along the outside edge of the antenna. To minimize any potential time delay (phase 
change) caused by the longer outside path length, the TEM horn is always several times 
longer than it is wide. Additionally, the two triangular plates of the TEM horn depicted 
in Figure 1 are always very acute-angled and isosceles. As long as the path length 
difference is less than one-half wavelength for all operating frequencies, then the resultant 


phase shift will be negligible. To find the frequency above which this will present a 


problem for our prototype TEM horn, we must first determine the path length difference 


from 


Al = T +(W/21) zi (2.3) 


where / 1s the length of the triangular plates, and W is the width of the TEM horn at the 
wide end of the triangular plates [4]. For our prototype, / = 0.9144m and W=0.1026m: 
therefore, AJ = 1.43789x10°m. The frequency at which A/ is one-half wavelength is then 
given by 


C 
=—, 2.4 


where c is the speed of light. For our prototype, this frequency is 104.32 GHz, which is 
well beyond the desired operating frequency range of this TEM horn. However, this is 
not the only potential frequency limit associated with the TEM horn. 

As long as the width-to-height ratio of the TEM horn is maintained throughout the 
entire length of the antenna, the antenna can be modeled approximately as a microstrip 
transmission line. Owing to this, we are then able to apply transmission line theory to 
accurately predict the upper and lower bandwidth limits of the TEM horn. This proves to 
be a tremendous advantage to the designer, since the only required data that must be 
known is the antenna's overall length, L, its width at the aperture, W, and total aperture 
height, 2h. 

To begin the process of determining the theoretical bandwidth limits of our 
prototype TEM horn, we return to the discussion of modes. As previously stated, the 
TEM mode is the only one desired in order to preserve the integrity of the incident fields. 
Thus, if any higher order transverse electric (TE) or transverse magnetic (TM) modes are 
excited, some distortion of the output waveform will occur. Fortunately, placement of 
conductors along the longitudinal axis of the antenna structure will prevent any induced 
currents from flowing across the width of the antenna structure (perpendicular to the 
antenna’s boresight axis). This, alone, will serve to sufficiently deter the excitation of 


any TE modes within the TEM horn. 


The same does not hold true for the TM modes, since there are no physical 
barriers, or mode filters, preventing their occurrence. From [3], only the TM,, mode, 
where the electric field has a component in the direction of propagation, has a significant 
probability of being excited, and thus, potentially causing distortion of the incident 
waveform. This mode can occur at any frequency higher than the limit governed by the 
equation 
ie 

2h 
where c is the speed of light, and 2h ıs the total aperture height. For our prototype TEM 


Fm, ~ (2.5) 


horn, this high frequency limitation caused by excitation of the TM,, mode will occur 
anywhere above 6.225 GHz. It is immediately evident that this is nearly a 40% reduction 
from the bandwidth we desire. Obviously, when we model our prototype TEM horn, this 
will be a frequency range to which we will need to pay close attention. If the TM,, mode 
is, in fact, excited at or near this frequency, then there will be a noticeable reduction in 
antenna efficiency (and possible distortion of the input signal) for all frequencies above 
this cutoff frequency, since power and signal transmission will now be passed through 
more than one mode. 

Unfortunately, simply lowering the value of 2h in Equation 2.5 may not be the 
best idea. Certainly, this will increase the cutoff frequency for the TM,, mode; however, 
decreasing 2h will also serve to decrease the overall sensitivity of the TEM horn [4]. A 
reduction in the sensitivity for a receiving antenna will more quickly degrade its overall 
effectiveness than will excitation of the TM,, mode. Thus, there is a trade-off between 
maintaining a high sensitivity and a high TM,, cutoff frequency, with sensitivity being 
the more critical factor. For this reason; modifying the prototype TEM horn to increase 
the TM,, cutoff frequency would not be a wise decision. 

On the low frequency end of our desired bandwidth, we find that the sensitivity is 
again the limiting factor. Specifically, it can be shown that under the assumptions given 


in [4], the low frequency limit can be expressed by 


E 


a (2.6) 


where fj «a, is the frequency where the sensitivity is 6dB less than its nominal (passband) 
value, Z is the overall length of the antenna, and c is the speed of light. For the given 
values of the prototype TEM horn, fisas = 13.054 MHz. This would imply the prototype 
has a 3dB frequency of 26 MHz. These will be areas that will need to be viewed with 
interest when we analyze the prototype TEM horn using the software package described 
later. Although the above is a theoretical calculation, we should expect some degradation 


in performance beginning in the 15-20 MHz range. 


D. LORENTZ RECIPROCITY THEOREM 


One final element that will prove critical not so much in the design portion of this 
research, but rather in the analysis of the near fields that follows, is the Lorentz 
Reciprocity Theorem. One result of reciprocity is that the radiation pattern of a particular 
antenna is the same whether it is used as a transmitter or receiver [2]. This is an 
important fact to remember, since it is sometimes easier to analyze an antenna simply as a 
transmitter than a receiver. Due to reciprocity, then, the radiation pattern and behavior of 
the antenna in a transmitting mode will be no different (save for the direction of 
propagation) than the corresponding antenna acting as a pure receiver. As stated, this will 
apply more to when we begin to analyze actual antennas. It is provided now to 
emphasize that even as we begin our discussion of the impulse receiving antennas, any 
discussion of their field patterns would be the same, even if they were then operated as 


transmitting antennas. 





II. ANALYSIS OF WIRE ANTENNAS AND SCATTERERS (AWAS) 


As critical as the VFLDVIS. M function is to the visualization of the near electric and 
magnetic fields (near fields), so, too, is the AWAS program critical to the actual calculation 
of those fields. It is this, and VFLDVIS. M's inherent dependence on AWAS output data, that 
makes the following discussion on AWAS necessary. While no substitute for the AWAS 
User's Manual, reference [5], this chapter should provide the user of NFLDVIS.M with 
enough detail and applicable background to create antenna geometries in AWAS and obtain 
valid output data. The information contained in this chapter is taken, for the most part, 
directly from the AWAS User's Manual. Included, though, is insight into AWAS as it 
applies to this thesis and to the user of NFLDVIS.M. 


A. BACKGROUND 


AWAS ıs a self-contained program which calculates the electric and magnetic fields 
(E- and H-fields) generated by a wire antenna or scattering structure in the near-field or far- 
field regions of space. As its name implies, AWAS is only able to calculate the fields for 
wire structures like the one shown in Figure 2. It generates the E- or H-fields from the 
currents that are induced along the wires. Wire antennas and scatterers are structures made 
of wire-like conductors, whose radii are much smaller than both their lengths and the 
wavelength of the operating frequency. Antennas which can easily be modeled by such 
wire structures are simple dipoles, V-antennas, wire loop antennas, tower broadcast 
antennas, Yagi-Uda antennas, log-periodic arrays, and antennas made by printed-circuit 
techniques. The wire structure can be placed either in free space or above a perfectly 
conducting ground plane. It can have concentrated or distributed resistive loads, and it can 
be driven at one or more ports in the transmitting mode or be excited by a plane wave in the 
receiving mode. AWAS is capable of analyzing these wire antennas operating in either 


mode individually or both modes simultaneously. [5] 
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Figure 2. eiit of Wire Antenna 


AWAS consists of a windows-based user interface for data input and output and a 
powerful kernel for numerical analysis of the wire antennas. A typical computation cycle 
consists of data input, calculation of the electric and magnetic fields, and creation of the 
data output files. The program solves the following integro-differential equation for the 
current distribution along the wires, referred to as the two-potential equation, 


N hm 
l dI u E, 
y» ju, le) fe, -— BDI 
0 k ds JOU 


m 0 





where m is the index of the wire segment; p is the index of the wire segment for which the 
boundary condition is imposed (p-1,...,N); /,, 1s the length of the m-th wire segment; /,(s,,) 
is the current on the m-th segment along the local axis, s, of that segment; g(r,) 1s the thin- 
wire kernel derived from the free-space Green’s function; r, is the average distance between 
the observation point and the points along the wire circumference; k = o£, u, ; and up & 


uj, are unit vectors along the wire segment's axis. The equation is based on the boundary 
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conditions for the electric field on the wire axes. It is solved numerically using the Method 
of Moments, with a power series approximation for the current distribution and pulse basis 
functions. [5] 

AWAS 1s available in three versions: AWAS for DOS, AWAS for Windows. and the 
Extended DOS versions. While AWAS for Windows is the most popular and better than the 
original DOS version, the Extended DOS version is by far the fastest and most robust of the 
three. Computation times using the Extended DOS version were found to be between 5 and 
10 times faster than the Windows version, regardless of complexity. The following 
discussion of all files and data entry blocks, however, 1s transparent to the particular version 
used, and as such any version can use the files generated by the other. In the case of this 
thesis, construction of all input files using the Windows version and actual calculation of the 
near fields using the Extended DOS version was found to be the most user-friendly. As 
such, all steps for proper data input contained in the next section are based on the AWAS for 


Windows version. 


B. DATA INPUT 


To begin the input phase of the program, first execute the AW AS program as per 
your particular operating system. Once open, you should see the main AWAS window as 
displayed in Figure 3. Either click on the “INPUT” icon on the toolbar or select the “Input” 


menu from the menu bar at the top of the window. 
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Figure 3. Main AWAS Window 
by 


The Input Data window should then appear in the lower portion of the screen, and 
the Show window should simultaneously open in the upper portion of the screen, as shown 
in Figure 4. The Input Data window is used to define all data necessary to create and 
specify the geometry of the wire structure to be analyzed. In addition, it is herein that the 
user will specify the frequency range and near-field group(s) over which the analysis will be 


conducted. The Show window provides a 3-dimensional (3-D) real-time view of the wire 


antenna as it is being constructed. 
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Figure 4. AWAS Input Data & Show Windows at Start-Up 








To begin the input process, select “Yes” on the Load Input Data window to modify 
any input files that have already been created, or select “No” to create a new antenna 


geometry, its operating frequencies, and the location(s) of the associated near-field group(s). 
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All data input and modification is accomplished within the Input Data window through 
either direct entry (as for the nodes and segments) or through a series of pop-up windows 
(as for the plane wave groups, near-field groups, impedances, and ports). Also, all input 
data should be in SI units (meters, amps, etc). 

In order to create the wire antenna structure, remember that it consists of a set of 
straight-line wire segments. Each segment connects two points in space (nodes), and each 
node is defined by its Cartesian coordinates (x,y,z). Figure 5 shows an example of how the 


input data is entered and assembled. The Show window displays a real-time sketch of the 


structure defined by the given nodes and segments. 
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Figure 5. Example of AWAS Data Input & Antenna Display 


Due to the specific relationship between the AWAS input and output files and 
NFLDVIS.M, certain restrictions have been placed on the construction of both the wire 
structure, itself, and the near-field groups for any analysis that will be carried over to 
NFLDVIS.M. First, the antenna must be oriented along the y-axis, i.e. the y-axis is the 
boresight axis of the antenna. Second, node 1 in the input .GEO file must be located at the 
source (feed) of the antenna, and the angle between the y-axis and a line drawn between 
nodes 1 and 3 must be equivalent to the angular dimension of the flared section of the TEM 
horn. Third, for all transmitting antennas, the input voltage at the feed must be +1 Volt to 
provide the "steady state" fields as output, and for all receiving antennas, the input plane 
wave should have an amplitude of 1 Volt/meter. These fields are then convolved with a 
Gaussian input impulse in NFLDVIS.M to provide the corresponding near fields due to the 
input impulse. Fourth, there must be an odd number of points (or none at all) along each 
axis of a near-field group. Fifth, all near-field groups must be symmetric about the y-axis 
and lie within either the xy-plane or the yz-plane. And lastly, no point within a near-field 
group may touch the antenna feed point or, to the greatest extent possible, any wires on the 
antenna structure. Because all near-field groups ın AWAS must be rectangular in shape, and 
in order to see the near fields within the flared section of the TEM horn, some points in a 
near field group created in the yz-plane will have to either lie on or outside of the wire 
structure. This will not create any problems as long as the restriction on nodes 1 and 3 (as 
discussed above) 1s met, since VFLDVIS. M will then be able to “filter” these data points out 
of the final near field calculations. Further data input limitations due to AWAS, itself, are 
discussed later in Section 3.D. 

Once the structure is complete, select “Compute” at the bottom of the window, and 
AWAS will check all nodes, segments, frequency steps, etc., for proper values and ranges. 
After checking the continuity of the wire antenna structure, AWAS then saves all input data, 


and calculates the near fields. 
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C. INPUT AND OUTPUT FILES 


All input data is saved within five separate files whose name extensions are GEO, 
FRO, PWE, NFP, and FFP. The antenna geometry data (node coordinates, segments, etc) 
are saved in files with the default .GEO extension. An example of one is shown in Figure 6. 
The first three lines of the file are the header, which can contain any description the user 
supplies. Each line is 76 characters long, including blanks. The entry in the fourth line 
indicates the operating mode selected for the particular structure. In this example, the “-1” 
indicates it is operating as a transmitter only, and only one port (feed) at a time is excited. 
A “+1” indicates a transmitting antenna with all ports excited, and a “2” would indicate a 
receiving antenna. The fifth line indicates the presence of a ground plane with a “1”. A “0” 
in the fifth line indicates the antenna is in free space. Following this, line six contains the 
number of nodes; line seven is the number of wire segments; and line eight is the number of 
ports (feeds). 

The next major group of entries defines the location of each node in Cartesian 
coordinates (x,y,z). In this example, the first node is located at point (0,0,0), the origin. All 
distances here are in meters, as 1s the same in the Data Input window. 

The next major group describes the location of each segment. The first two 
characters in each line represent the start and termination nodes, respectively, of the 
segment. The third character is the polynomial degree chosen for calculating the two- 
potential equation, and the fourth number is the radius of the wire segment (in meters). The 
final three complex number sets indicate the concentrated impedance at the start node (Zg), 
the concentrated impedance at the termination node (Zp), and the impedance uniformly 
distributed along the entire segment (Zq), respectively. 

The final group corresponds to the ports on the antenna. Each line represents the 
segment number. complex driving voltage, and nominal impedance of each particular port. 
If the port is located at the origin of the segment (as 1s the case for the example shown), then 
the first integer will be positive. If the terminal end of segment #1 were driven for this 


example, then the first integer would be a ^-1". 
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Figure 6. Example of AWAS Antenna Geometry (.GEO) Input File 





The structure of the .GEO file is critical to NFLDVIS.M, as the function reads and 
extracts the pertinent data necessary for reconstructing and visualizing the near fields using 
MATLAB and its associated graphic tools. 

The .FRO file extension indicates those files which contain the data defining every 
frequency for which AWAS will analyze the antenna structure. This file exists solely to pass 
the frequency data to the AWAS program kernel and is not used by the NFLDVIS.M 
function. Likewise, the .PWE, .NFP, and .FFP files contain the data which define the 
incident plane wave(s), the near field zone(s) of interest, or the far field zone(s) of interest to 


the user. They, too, are used strictly by the AWAS program kernel. 


The results of the AWAS analysis are written to seven separate default output files. 
At the user's discretion, these files are then also saved under a user-defined name with 
seven different file extensions: OUT, SPM, YPM, ZPM, CUR, NFL, and FFL. The only file 
extension that will be discussed 1s the .NFL file, since this is the one which contains all the 
near E- and H-field data for the desired frequency range, and thus, is the primary source file 
for NFLDVIS.M for visualization of the near fields. An example is shown in Figure 7. 

It should be noted that the .VFL file does not contain a header as its first line. This 
was inserted using a text editor exclusively for use within this discussion. The actual first 
line of every .NFL file begins with an integer representing the number of near field groups 
analyzed by AWAS (2 for this example). 

Following this line, the second line contains the start and stop frequencies, f, and f, 
respectively, followed by the number of frequency steps, M. It should be noted that the 
frequency step size, df, should always be set equal to the starting frequency, fi, and the 
number of frequency steps must be at least 16, if the user wishes to view the fields in the 
time domain. For the example, the structure was analyzed over the frequency range of 
SMHz through 30MHz in 5MHZz steps (5 steps total). 

The third line indicates the antenna type in the same manner as indicated in the 
.GEO file. Our example here represents a transmitting antenna exciting one port at a time. 
The fourth line indicates the number of ports present on the antenna structure (1 for this 
example), and the next line is the number of plane-wave groups incident upon the structure 
(0 for our pure transmitter, as common sense would dictate). 

The next block of data is the first array of near-field values. All of the following 
blocks will look similar to this one. This first block of near-field data applies to the first 
near-field group and represents the fields calculated at frequency ff. The next block would 
represent the data for the second near-field group at frequency fi, and so on until the near 
fields are calculated for every near-field group at this first frequency. After that, the data 


blocks cycle through all the near-field groups again at each successive frequency. 
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-6 377 1E-07 -2.7520E-05 -1.0013E-14 -1.4970E-12 2.5240E-15 -4.5394E-14 
-2.0350E-10 -5.4072E-12 -1.6062E-02 7.2755E-05 -1.6864E-02 -1.0114E-03 
-6.6637E-07 -3.3259E-05 -1.1727E-14 -1.2931E-12 3.4174E-15 -5.7776E-13 
-19071E-10 -4.3130E-12 -13193E-03 1.8652E-06 -2.1435E-02 -1.0011E-03 
-6. 7691 E-07 -3.5852E-05 -1.2837E-14 -6.2150E-13 3.0011E-15 -8.6237E-13 
-4.2832E-11 -3.6877E-12 -1.3618E-02 9.3978E-05 -9.5490E-04 -1.0658E-03 
-7.0899E-07 -1.5349E-05 -1.1133E-14 -6.8850E-13 1.5925E-15 1.7261E-13 
-9 399 1E-11 -2.2524E-12 -1.2271E-02 7.5159E-05 -5.4630E-03 -1.0479E-03 
-7 3289E-07 -1.9021E-05 -1.1834E.14 -7.3688E-13 2.4303E-15 -5.4179E-14 
-1.0624E-10 -4.9027E-12 -77 557E-03 4.2078E-05 -9.5170E-03 -1.0341E-03 
-7 4947 E-07 -2.1879E-05 -1.1970E-14 -6.3460E-13 2.5153E-15 -2.5416E-13 
-7 3385E-11 -1.7550E-12 -7.8768E-04 1.1382E-06 -1.1435E-02 -1.0289E-03 
-7.5557E-07 -2.3104E-05 -1.3475E-14 -3.3662E-13 3.6603E-15 -3.7553E-13 





(x1,y1,z3) 
(x1,y1,z4) 
(x1,y2,z1}) 


(21.72.12) 


-4 2206E-11 -1.1084E-12 3.5108E-03 -2.5929E-05 -6.0404E-03 -1.0334E-03 
-8.5918E-07 -1.5868E-05 -1.3279E-14-6 4365E-14 3.1607E-15 -1.8757E-13 

-6.6279E-11 -2.4087E-12 6.5167E-03 -4.8680E-05 -4.1696E-03 -1.0395E-03 Narren 
-8.4867E-07 -1.4442E-05 -1.4462E-14 3.9219E-14 3.2370E.15 - 1.39 56E-13 near field group #2 
-45475E-11-10587E-12 8.0195E-03 -6.4194E-05 -1.7180E-03 -1.0469E-03 at frequency fl 

-8 3278E-07 -1.2384E-05 .1.5501E-14 9.4771E-14 2.0157E-15 -6.9944E-14 
-39989E.11-17764E-12 2.0474E-03 -1.8961E-05 -3.8032E-03 -1.0190E-03 

-9 7 782E-07 -1.2004E-05 -1.3663E-14 -6 9347E-14 2.0817E-15 -1.0847E-13 


Figure 7. Example of AWAS Near-Field (.NFL) Output File 
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Within each data block, the first two lines represent the E- and H-fields, 
respectively, for the first planar (spatial) coordinate (x,,y,.z,). Lines 3 and 4 of that block 
represent the E- and H-fields for the spatial coordinate (x,,y,,z,), and so on through each 
spatial coordinate of that near field group. It should be evident that the data cycles through 
the z-coordinates, then the y-coordinates, and finally the x-coordinates, in that specific order. 

Looking across the first line of near field data, the E-field is presented as six 
floating-point numbers. These represent the real and imaginary components of the E-field 
in the x-direction (Ey), the real and imaginary components of the E-field in the y-direction 
(Ey), and those of the E-field in the z-direction (£7), respectively. Similarly, the second line 
contains the H-field components along the three Cartesian axes, Hy, Hy, and Hz. It should 
then be evident that for every block of near field data, all odd rows contain the E-field data 
and the even rows contain the H-field data. This specific structure is critical to the function 
NFLDVIS.M's ability to access and manipulate all the near field data for graphical analysis. 
Consistent here, too, is the fact that all quantities are given in SI units (meters, Hz, 
volts/meter, amps/meter, etc). 

For the case of the receiving antenna analysis, the .VFL output file differs only 
slightly from its counterpart for the transmitting antenna case, Figure 7. As noted in Figure 
8, line 3 will have a “2” representing the receiving antenna, while line 5 indicates the 
presence of a singular incident plane wave. In addition, prior to the block(s) of data 
representing each near-field group, an extra block of data 1s inserted which describes the 


incident plane wave information, as shown. 
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Receiving 


T 
1.0000E+07 


Antenna ~ 2 


# of Plane 
Waves 


Plane Wave 
Info 


—; 


1.0000E+00 
90.00 
90.00 
-2.0000E+00 
1:00005-07 
0.0000E+00 


0.0000E+00 
-=7.2 713 E09 
ZOSI SZET 
>Il SORE 
0.0000E+00 
—128786E-0B 
0.0000E+00 
z2515325-08 
z2999132E-10 
—3202 0E-08B 
=9, 3132-10 
z3 5998508 
0.0000E+00 
=4,1715E-08 
0.0000E+00 
-4.7419E-08 
0.0000E+00 
295 41095-08 
0.0000E+00 


90, 
907,00 0 


1.0000E+09 


0.0000E+00 
00 


2.0000E+00 
2.0000E+00 
0.0000E+00 


0.0000E+00 
3.021207 
0.0000E+00 
-5240 7856-07 
0.0000E+00 
=7.10625-07 
0.0000E+00 
=92 9096607 
2.9104E-11 
=12 1999E-06 
0.0000E+00 
ie 394 IE DE 
0.0000E+00 
1257 23E 06 
0.0000E+00 
-1. 7364E-06 
0.0000E+00 
-1.8822E-06 
9.B20BE-II 


99 


0.0000E+00 
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0.0000E+00 
= 4551-07 
3.8208 11 
-1.4548£-07 
0.0000E+00 
-1.4544E-07 
0.0000E+00 
=1.4537E-07 
2.3283E-10 
1215252 0, 
2.32835-10 
-21:24520E-07 
0.0000E+00 
-1.4509E-07 
0.0000E+00 
-1.4497E-07 
0.0000E+00 
-1.4484E-07 
0.0000E+00 


RCVR_EXAMPL.NFL 


0.0000E+00 


0 k— (Biar > Osrop ; # of steps) 
«— (o > stop > # of > steps) 
50 


0.0000E+00 
=8,.0823E-06 
0.0000E+00 
6 035665806 
0.0000E+00 
22 96elE 06 
0.0000E+00 
-5.8724B8-06 
SAEZ 
za ME SUE 
0.0000E+00 
=5.62545-U06 
0.0000E+00 
=5,4729E06 
0.0000E+00 
==. 086506 
0.0000E+00 
z8:1332E-06 
Z2 B28E-LDIT 


-——— (Es real > En , rca , En) 


-08278-03 
.0000E+00 
.0510E-03 

. 0O000E+00 
-00215-03 

. 0O000E+00 
93 10203 

. OOOOE+00 
285 115203 
0. 
E 
0. 
m 
0. 
zer 
. OOOOE+00 
„42105803 
.0000E+00 
»2929E 03 


0000E+00 
7638E-03 
0000E+00 
6589E-03 
0000E+00 
5442E-03 


JOLOBE-Ud 
. 0O000E+00 
.9164E-04 
.0000E+00 
5153-04 
.0000E+00 
.5149E-04 
.0000E+00 
2:51395-04 
. OOOOE+00 
25125E-04 
. 0O000E+00 
.9110E-04 
. DO0O0OE+OO 
:9092E5E-04 
.0000E*00 
. 5071E-04 
.0000E+00 
.5048E- 04 





Figure 8. Example of AWAS Near-Field (.NFL) Output File 


PROGRAMMING NOTES 


Table 1 lists the upper limits of the various input data for each of the three AWAS 


versions. It should be noted that besides being slower than the Extended DOS version, 
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AWAS for Windows only allows a maximum of 199 frequency steps, a limit that severely 


handcuffs the UWB antenna designer. 


Table 1. Upper Limits for All AWAS Versions 


aaraa m e m mn t t Ma m Mm HH ARAB nm Raph FRI ARP nn) ne Sn an Al Ad 


Item | AWASfor | Extended | AWAS for 

Windows DOS DOS 
Max + of Nodes | 50 | 99 | 50 — 
S tu Seien A A RA E 
Max # of Ports — 


Max # of Freq Steps — 


Max # of Plane Wave Groups 
Max # of Near Field Groups 
Max # of Far Field Groups 


Max 8 of Coordinate Steps per Group 
Max Polynomial Degree 
Max # of Unknowns — 0 





As stated in the AWAS User’s Manual, all input data supplied to AWAS is further 
restricted to certain ranges in addition to the above maximum limits. These restrictions on 
the input data are as follows: 

e The number of nodes must be at least 2. 

e The number of segments must be at least 1. 

e The number of ports must be between 1 and 16 for an antenna. 

e The z-coordinate for a node must be > 0 if the ground plane is present. 


e The node indices for a segment must be between | and the total number of 
nodes. 


Zi 


The segment radius must be positive. 


The polynomial degree must be 0 for automatic segmentation or between 0 and 
9 for manual segmentation. 


The real part of a loading impedance must be 2 0. The imaginary part can be 
positive, negative, or zero. 


The segment index for a port must be between 1 and the total number of 
segments. 


At least one port’s electromotive force must be > 0 for a transmitting antenna. 
The port nominal impedance must be > 0. 

The start and stop frequencies must be positive. 

The number of frequency steps must be at least 1. 


The number of plane wave groups must be between 1 and 8 for receiving 
antennas. 


The angle @ for the plane waves must lie between -180° and 180° (or -90°and 
90° if the ground plane is present). 


The angle q for the plane waves must lie between -180° and 180°. 
The number of near field groups must be at least 1. 


The upper and lower z-coordinate limits in the near field groups must be > 0 if 
the ground plane is present. 


The number of coordinate steps in the near field groups must be = 0. 


The information presented within this chapter is by no means meant as a 


substitution for the AWAS User’s Manual; however, it should give the antenna designer 


enough information to quickly create, load, and compute the near E- and H-fields for a 


particular antenna structure. Once those near fields are calculated, though, it 1s the function 


NFLDVIS.M which takes that data and transforms it into a visual display for complete 


graphical analysis of the antenna’s effectiveness at any given frequency. 


IV. VISUALIZATION OF A PROPAGATING IMPULSE 


In order to visualize the transient near-fields produced by an UWB antenna, and 
more specifically, by the TEM horn, a graphical user interface (GUI) was developed using 
the MATLAB programming language as its base. The GUI and the entire algorithm 
developed to conduct the necessary matrix calculations were written as the function 


NFLDVIS.M, included here as Appendix A. 


A. SOFTWARE OVERVIEW 


NFLDVIS.M, heretofore referred to as “the software,” was developed as the key 
analysis tool for the UWB antenna designer. Its system requirements are a 486DX or faster 
processor with Windows 95 or Windows NT as the operating system, an 8-bit graphics 
adapter, and 16 MB of Random Access Memory (all MATLAB 5.x constraints). The actual 
processor speed, amount of Random Access Memory (RAM), and graphics requirements, 
though, will be dictated by the size of the near-field arrays calculated within AWAS. The 
system used in the development and testing of the software was a 200MHz Pentium Pro 
with 132MB of RAM. This proved adequate for manipulation of arrays as large as 
101x101x1024, where the first two dimensions represent the physical size of the near-field 
array in the xy- or yz-plane, and the third dimension is the number of time steps (directly 
derived as twice the number of frequencies analyzed). 

The software must be run under MATLAB 5.0 or higher, due to the extensive use of 
3-dimensional (3-D) matrices and also the integration of new graphics and user-interface 
controls unavailable in previous MATLAB versions. Although a complete User's Manual 
Is not included, the highlights as they apply to this thesis will be presented and reinforced 
throughout the main body of this document. 

The software was designed to specifically work with TEM Rome or wire V-antennas 
which were constructed and tested using AWAS and physically oriented along the y-axis in 


3-D space. Due to this physical orientation, only the y- and z-components of the electric 


field and the x-component of the magnetic field are ever made available for full analysis and 
visualization (Ez and Hy are the major near-field components, while any values of Ey 
represent noise in the direction of propagation). While this may at first appear to be the 
most simple or trivial scenario for an antenna alignment, it should be noted that all analysis 
here can be easily manipulated for application to any direction of propagation through the 
use of "direction cosines". Simply put, what you actually visualize here is a true 
representation of the near fields along the antenna's boresight axis. 

The software does, however, present certain limitations to the designer. As 
mentioned in Chapter III, AWAS will only construct rectangular near-field grids. This 
presents a problem any time the designer wishes to view the near fields in the yz-plane, 
since the TEM horn is not rectangular throughout its entire shape. In order, then, to filter 
out the points on or outside of the antenna structure (especially near the feed), the software 
calculates the half-angle dimension of the TEM horn, then simply replaces those values 
outside the designated angle with zeroes. To do this, it is assumed that the angle between 
the y-axis and a line drawn through the 15% and 374 nodes in the geometry (.GEO) file 
accurately duplicates the angular dimension of the flared portion of the TEM horn. This 
limitation should pose little in the way of additional effort on the part of the designer, but 
this knowledge is essential prior to construction of any antenna geometry in AWAS. This 
limitation does not apply for the case of a near-field array lying in the xy-plane due to the 
fact the TEM hom has no physical sides, and thus, field points within the rectangular array 
will not lie on or near any physical part of the antenna. 

Throughout the course of this chapter, the terms “near-field groups” and “near-field 
arrays” will be used. As a matter of clarification, “near-field groups” are the physical xy- or 
yz-plane limits imposed by the user within which AWAS will calculate the near E- and H- 
fields. “Near-field arrays,” or “near fields,” are the matrices which contain the field strength 
values for each point within a near-field group and over the entire range of frequencies. 
Whereas a near-field group is a 2-D planar matrix containing x- and y-coordinates or y- and 


z-coordinates in space, a near-field array is a 3-D matrix containing the E- or H-field 


magnitude that AWAS calculated for each point within the near-field group and over the 
entire user-defined frequency range for which the antenna was tested. Although this fact is 
not critical to the user of the software, its knowledge will aid in understanding how the 


software is able to create the graphical analysis that is fundamental to the design process. 


B. INITIALIZING THE SOFTWARE 


After opening MATLAB to the correct directory and typing "nfldvis" at the 
command line prompt, the entire screen will fill with a figure window containing the GUI. 
The GUI’s purpose is to provide a media through which the antenna designer/software user 
may conduct all analysis of an antenna, as well as plotting, printing 3-D figures, and saving 
files entirely from within the GUI window and with only the use of a computer mouse. The 
GUI consists of two sets of axes for plotting both the data and an optional colorbar, the user 

interface controls (uicontrols) such as the push-buttons and pop-up menus, and the user 
| interface menus (uimenus) such as the "Plot" and "Print" pull-down menus, as depicted in 
Figure 9. It should be readily apparent that some uicontrols have darker text than others. 
Within uicontrols and uimenus, the text is dark when the uicontrol/uimenu 1s activated, and 
the user may press or use it. Light or nearly transparent text indicates the uicontrol/uimenu 
is deactivated, and it will not function, regardless of the user's attempts. In the case of the 
sliders at the bottom of the window, the sliding button disappears when these are 
deactivated. 

Immediately upon initializing the software, the user is provided with four options 


through which to proceed: 


l. Select "Start," and the user may start the analysis of a particular antenna; 


N 


Press the "Demo" pushbutton, and a complete demonstration (demo) of the 
software will run; 


3. Select "Movie," and either a new movie can be created using the data from the 
last analysis, or a previously-created movie may be replayed, or; 
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4. Use the "Plot" uimenu to plot the frequency-domain fields of the last antenna 
analyzed. 
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Figure 9. Initial Window on Start-Up of NFLDVIS.M 


The “Reset” and “Close” pushbuttons are available at all times, and while the 
function of the “Close” button should be obvious, the “Reset” button’s function 1s 
exclusively to restore all uicontrols to their default values while simultaneously returning 
the user to the screen shown in Figure 9. The new user of NFLDVIS.M should first select 
the "Info" button and then run the demo before beginning the analysis of a particular 
antenna in order to acquaint him/herself with the software, its capabilities, and the normal 
flow of the data and graphics subroutines. 

The "Info" button and its associated subroutine specifically serve to introduce the 


new user to the software and its capabilities and limitations. It is merely a quick synopsis 


and by no means a substitution for the User's Manual, reference [6], but it does provide the 
user with the requisite knowledge to get started on the nght track. 

Pressing the "Demo" button activates the demo package. This is provided to 
illustrate the inherent capabilities of the software and familiarize the new user with the 
general flow of analysis, data manipulation, and graphics presentation. It does not, 
however, actually conduct any calculations of the fields; rather, it simulates such data 
analysis to ensure the demo presents the entire graphical-analysis package without being too 
time-consuming. It is constructed ın a user-friendly manner that waits for cues from the 
user prior to moving forward to the next step. The entire demo may be run in a matter ofa 
few minutes, depending solely upon the delays desired by the user. The general flow of the 


demo consists of: 


l. File selection by the user; 

2. Simulated near-field array construction; 

3. Frequency-domain plotting of the E-field; 

4. Application of the gaussian impulse; 

5. Frequency-domain plotting of the E-field modified by the input impulse; 
6. Transformation of the frequency-domain E-field to the time domain, and; 


7. Creation of a movie representing the time-domain representation of the E-field, 
e.g., the propagation of the gaussian impulse through the antenna array. 


The specific methodology used to create the data arrays is contained in Section 


4.A.3, while the entire upper-level graphics package is covered in detail in Sections 4.A.5 


and 4.A.6. For a complete step-by-step discussion of the demo package, refer to reference 


[6]. 


U) 
u) 


C. DATA INPUT & CALCULATIONS 


The actual data input and calculations begin by selecting the "Start" button. 
Immediately, the user is tasked with selecting the appropriate 4WAS .GEO and .NFL files 
for the desired antenna, as shown in Figure 10. It is critical that both files should 
correspond to the exact same antenna. Normally, this means the two files will share the 
same name and differ only by their file extensions (e.g., TEMI.GEO and TEMI.NFL), since 
AWAS saves all input files in this manner and prompts the user to do likewise with the 


output files. In addition, both files must reside in the same Windows/DOS sub-directory. 
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Figure 10. Selection of .NFL File for Analysis 


Following selection of the appropriate files, the software extracts the necessary 
physical dimensions of the antenna and the near-field groups, and then asks the user to 


determine which near-field group to analyze and plot. It then reads the near-field data from 


34 


the .VFL file and creates separate 3-D matrices in MATLAB containing the field strength of 
the y- and z-components of the E-field and the x-component of the H-field in the frequency 
domain (Ey(f ), Ez(f ), and H,(f ), respectively). For the case of a near-field group 
constructed in the xy-plane, the first dimension of the 3-D array will be the number of points 
in the x-direction, the second dimension will then equal the number of points selected along 
the y-direction, and the third dimension will be the number of frequencies tested. The size 
of the final array constructed for the yz-plane near-field group, likewise, will be determined 
by the number of points in the y- and z-directions and the number of frequencies tested. 
However, if the ground plane was present ın the AWAS calculations, the software 
compensates for this in the construction of the matrices and creates the necessary data for 
the near-field array as if the ground plane were not present, 1.e., it constructs a mirror image 
of the AWAS-calculated array and merges the two together along the y-axis. Regardless of 
the type of near-field group, though, the 3-D near-field arrays will always be symmetrical 
about and include the y-axis. 

After the frequency-domain near-field arrays are assembled in MATLAB, the user 
is prompted to provide identifying information that the software will attach to the title of 
every plot or movie generated from the computed near-field arrays. The software then 
creates a Gaussian impulse input wave, Figure 11, whose maximum value in the time- 
domain is 1 Volt. The pulsewidth is a function of the minimum frequency and the total 
number of frequency steps. For this example, the pulsewidth as measured at the 3dB points 
is 1.0 nsec. This pulse is then convolved with the frequency-domain fields by taking the 
Fast Fourier Transform (FFT) of the time-domain Gaussian pulse, Figure 12, and then 
multiplying this result with each near-field component. As Figure 12 illustrates, the 
frequency-domain representation of the Gaussian pulse is similar to that of a low-pass filter, 
albeit a very broadband one [7]. The net results, as shown in Figure 13, are near-field arrays 
in the frequency domain that conform to the shape of Figure 12. These “filtered” frequency- 
domain E- and H-fields will, through use of the Inverse Fast Fourier Transform (IFFT), 


convert easily to an accurate time-domain representation of the near fields within and in 


close proximity to the antenna that result from the application of a Gaussian impulse 


voltage at the antenna port. 
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Figure 11. Gaussian Impulse Input Voltage 
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Figure 12. |V, (f), The Frequency Spectrum of V. (t) 
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Figure 13. Example of "Filtered" Frequency-Domain E-Field 


After the IFFT is applied along the 3” d dimension of the 3-D arrays, and only if the 
selected near-field group was in the yz-plane, the software will further filter the near-field 
arrays to eliminate the field values for points next to, on top of, or outside of the wires 
representing the flared portion of the TEM antenna structure. This last action completes the 
string of calculations and matrix manipulations necessary to create the variables that will be 
used for the various plotting and movie subroutines. The user is informed of the fact that all 
calculations are complete and then prompted to continue the analysis by either plotting the 
frequency-domain fields and the Gaussian input voltage or creating a movie of the near- 
fields in the time domain. The user may also elect to end the analysis there temporarily and 
return to complete his work at a later date, since all the variables needed to execute the 
plotting and movie functions were already saved automatically during the “Start” 


subroutine. 


D. OUTPUT FILES 


Rather than hold potentially large variables for the entire process of constructing the 
numerous near-field arrays, the software simultaneously saves all the critical variables 
within default .MAT files in the current directory and then deletes the variables from the 
workspace in order to free more memory for active matrix calculations. By saving the 
variables as .MAT files, they can then be readily accessed at any time in the future for 
plotting or movie construction. The eleven separate .MAT files are Eyf. mat, Ezf.mat, 
Hxfmat, Eyffmat. Ezffmat, Hxffmat, Hdr.mat, Eyt.mat, Ezt.mat, Hxt.mat, and 
Movihdr.mat. The first three contain the frequency-domain near-field arrays for each of the 
respective near-field components and are saved once the arrays are constructed from the 
AWAS .NFL data file. The second set of three files contain their respective “filtered” near- 
field arrays (the frequency-domain fields due to the input Gaussian impulse), and Hdr.mat 
contains a number of minor variables used to either support plotting of the frequency- 
domain fields or of the Gaussian impulse. These files are saved prior to transforming the 
arrays from the frequency-domain to the time-domain via the IFFT. The three time-domain 
near-field arrays are saved within Eyt.mat, Ezt.mat, and Hxt.mat, and Movihdr.mat contains 
the minor variables required to support construction of a movie. They are the last files 
saved prior to the “Start” subroutine concluding. 

The above .MAT files form the data sets that both the "Plot" and "Movie" 
subroutines use for graphics display. Since these files can easily exceed 5MB in size, it is 
much easier to load and clear them as needed, as opposed to creating them through the 
“Start” subroutine each time a plot must be rendered. Also, the .MAT files will hold the 
near-field arrays and the associated minor variables until either the user deletes them 
manually or another antenna is analyzed. Because of this, the user can then conduct the 
analysis at his or her own pace. Also, as will be seen in the next two sections, the user may 
manually rename any or all of the aforementioned .MAT files using the operating system’s 
file manager protocol, and thereby ensure the data never requires recalculation. This makes 


things easy when comparison of two antenna models is required/desired. 


E. GRAPHIC DISPLAYS 


Every plot of the frequency-domain near-fields or the input Gaussian function can 
be found beneath the "Plot" uimenu on the GUI's menubar, as seen at the top left corner of 
Figure 9. All of the plots are grouped within the three sub-menus “Frequency-Domain 
Fields”, “V”, and “Filtered Freq Domain Fields”. Immediately after completion of the 
“Start” subroutine or anytime there exists data within either the default output files or other 
user-defined files as described in Section 4.A.4, the user may access and display the data 
merely by moving the mouse-pointer atop the “Plot” uimenu, clicking the left mouse button 


once, moving the mouse-pointer to the desired plot’s sub-menu, and then clicking the left 


mouse button again, as illustrated in Figure 14. 
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Figure 14. |. Example of E-Field Frequency-Domain P Plot 


In the example shown, a frequency-domain plot of the E-field as a function of both 
frequency and the y-axis position within the near-field group is rendered. Specifically, this 


is a plot of the magnitude of the E-field along the boresight axis of the antenna and over the 


entire frequency range tested. The plot's colors range from blue at Ez-0 (cold) to red at 
E7-max(E;) (hot), and the exact color range and associated numerical values are duplicated 
in the colorbar to the right of the 3-D plot. If the user desires a different perspective angle 
to view the plot, he need only to place the mouse-pointer on top of the plot, press and hold 
the left mouse button, then move the mouse in any direction desired. As shown in Figure 
15,a3-D frame will first appear enclosing the original plot, and then it will move and resize 
according to the movements of the mouse by the user. Rotation is completely unlimited in 
all three dimensions; therefore, the user is able to view the 3-D plot from any perspective. 
This frame is merely a wire rendering of how the plot box will be reshaped once the user 
releases the left mouse button. Although not displayed in Figure 15, a clear text overlay 
will appear at the bottom left-hand corner of the GUI and provide real-time azimuth and 
elevation data regarding the view angle of the plot as the mouse and frame move. The 
azimuth and elevation text boxes at the bottom center of the window, however, are used 
within the movie subroutine only and will not provide any information to the user during 
any of the plot subroutines. 

All of the plot subroutines within either of the two frequency-domain sub-menus 
will generate their plots in the same manner as depicted in Figure 14 and 15. The only 
difference will be the type of data used to render the 3-D graphic. Also, every plot will 
have the user-provided identifying information appended to the title. In the example 
shown, the user wished to identify this plot as that for the yz-plane near-field group from 
an antenna called “TEM3C,”’. 

The user may also elect to view the input Gaussian impulse voltage function in 
either the time- or frequency-domains. Selection of either graphic is no different than 
before, but this time, the rendered plot is only 2-D and cannot be rotated. An example is 
shown in Figure 16. The pulse width and total length of V, the input voltage pulse, are 
determined by the maximum frequency value and the total number of frequencies defined 


by the user in AWAS. 
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Figure 15. Example of 3-D Plot As View Angle Is Altered 


In all, each of the plots allows the user to view the near fields along the centerline 
axis of the antenna for each frequency both before and after the Gaussian impulse is 
applied. The user may further view the input waveform directly, too. It should be noted 
that the fields viewed before application of the impulse are no different from the fields that 
would occur in the steady-state following application of a 1 Volt sinusoidal voltage signal at 
the antenna feed. 

The frequency domain, while handy in determining the effects the antenna has on 
the generation of the near fields at each particular frequency, still is not able to answer the 
questions of “how much of the pulse's energy actually couples and propagates beyond the 
antenna?”, and “to what degree is the pulse's waveform altered as it propagates through and 
out of the antenna?". These, and other questions, are best determined in the time domain. 


and that is where the “Movie” subroutine plays an integral role in the analysis process. 
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Figure 16. Example of Gaussian Input Voltage Plot 


F. MOVIE SUBROUTINE 


Once the “Start” subroutine 1s complete, or at anytime there exists valid data in 
either the default output .MAT files of Section 4.A.4 or other user-named files of the same 
configuration and design, the user may create a movie to visualize the near fields in the time 
domain. After depressing the “Movie” push-button, the user is asked to choose between 
creating a new movie of the time-domain fields or to replay a movie that had previously 
been created and saved. We will begin the discussion with the methodology for creating a 
new movie and hold off on the discussion of replaying previously-created movies until the 
end of this section. 

After electing to create a new movie, the user is asked to select the file which 


contains the desired component of the time-domain fields and the associated header file 
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that contains the minor variables needed to work with all three of the time-domain fields 
here. The software prompts the user with the correct default filenames for each .MAT 
file; however, if these files were manually renamed by the user earlier, then those files, 
too, can be loaded here in place of the default .MA7 files. 

After the files are loaded, the user 1s shown the values for three variables: 7, N, 
and df. The variable T is the total time duration of the movie, in usec. It does not 
represent the time it will take to view the movie, but rather the timeframe that the movie 
Is depicting. The movie is similar to a slow-motion instant replay of the effects of the 
Gaussian input pulse over the entire timeframe 7. N represents the maximum number of 
frames the user may use in designing the movie. It also represents the total number of 
time steps between 0 and 7. Finally, the variable dt represents the length of each time 
step (in nsec). At the same time the values of these variables are shown to the user, the 
software also prompts the user to first select the desired number of frames for the movie 
from the “Total Frames” pop-up menu, and then select the frame rate (in frames-per- 
second, fps) from the “Frames-per-sec” pop-up menu. 

After storing that information, the software then renders a sample frame from the 
movie, as depicted in Figure 17. Its purpose is to provide a “snapshot” from the movie for 
use in tailoring optimum/desired viewing angle and scale factors. Like the 3-D frequency- 
domain plots of Figure 14 and 15, the 3-D image projected here can also be rotated freely 
about itself in the manner described in Section 4.A.5. Unlike there, the azimuth and 
elevation windows at the bottom of the GUI will display the current view angle after every 
complete update cycle (after the “Continue” button is pressed). In addition, this figure can 
be scaled along the z-axis in order to provide better resolution or greater emphasis upon a 
specific area of the pulse of the antenna. The z-axis scaling “slider” is located at the bottom 
left corner of the GUI. It has a range of 0 to 1, representing 0 to 100% of the maximum 
value of the particular near field, since the default z-axis scale is set to the maximum value 
of the particular near field being plotted. This value will always occur near the feed, since 


the cross-sectional area of the antenna at that point is the smallest. Therefore, it is 


recommended that for use with any thin antennas, the z-axis should be scaled down. 


Reducing the scale of the z-axis does not change the data that will be plotted; it merely 


alters the range of the z-axis that will be plotted within the entire plot-box, or display. 
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Figure 18 displays this, and, as compared to Figure 17 


Gaussian pulse is much more clearly defined. Also, as with the azimuth and elevation 


windows at the bottom of the GUI, the rendering of the plot will only be updated with the 
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Figure 17. Example of Initial Sample Frame from Movie Subroutine 
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Figure 18. Sample Frame from Movie after Modification of Z-Axis and View Angle 


Another variable that can aid in better graphical analysis of the time-domain fields is 
the modification of the color-axis (c-axis). By adjusting the “Color-Axis Scaling” slider at 
the bottom of the GUI, the user ıs able to scale the range of the numerical values assigned to 
the colormap. The colormap, as depicted in the colorbar to the right of the plot in Figure 
18, consists of 128 separate colors which range from "hot" red at the minimum value to 
"cold" blue at 0 and back to "hot" red at the maximum value. The default values assigned 
to the colormap are the same as those assigned to the z-axis. In much the same manner as 
for the z-axis scaling, then, scaling of the c-axis will provide better resolution of a particular 
area. In the case of color, though, a subtle change in the amplitude of the pulse as it 
propagates can produce a dramatic color shift, depending upon how well the c-axis is 
scaled. Figure 19 is our original sample plot with the c-axis scaled for better emphasis of 


the pulse. By scaling the z-axis and c-axis separately, the user may keep a majority of the 


maximum value of the pulse within the figure window while also emphasizing an area of 
interest with the color palette. 

As further seen in Figure 19, after the user updates/changes the sample plot and 
depresses the "Continue" button, the software will automatically update the plot and then 
ask the user if he would like to modify the plot again. Until the user selects ^No", the 
process of modification and updating of the plot will continue ad infinitum. Once the plot 
is completely customized to the user's desires and the adjustments cease, the software then 
builds the movie frame-by-frame, plotting the near-fields for each time step. Once all 
requested frames are constructed, the movie will run automatically at 12fps and then at the 
frame rate selected by the user. After the movie runs through the playback mode (at both 
speeds) once, the user may select whether to replay the movie as many times as he/she 


wishes. 
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Figure 19. Sample Frame from Movie after Modification of the C-Axis 
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After the final replay, the user once again is given the option to make adjustments to 
any of the movie parameters and rebuild it again. If this option is selected, the movie is 
erased from memory, and the sample frame with current view angle and default scaling 
factors is again displayed for modification as outlined earlier. If the rebuild option is not 
selected, the user is then prompted to provide a filename under which to save the movie for 
replay at a later time. 

In order to replay a movie which has been previously created and saved, the user 
need only to depress the “Movie” push-button anytime it 1s not deactivated. The user is 
then prompted to select the filename under which the movie may be found and further 
asked to select a frame rate. Once that is done, the software loads the movie and runs it 
twice - once at 12fps, and once at the frame rate selected by the user. The user can replay 
the movie ad infinitum without reloading it each time, but no parameters may be changed, 
since the saved movie is merely a raster-scan image (like a photograph) of each frame. This 
prevents the movie axes or view angle from being modified. 

The movie subroutine is, by far, the more robust tool of this software package. Its 
importance should become even more evident as we move to analyze specific antenna 
structures on the way to designing an UWB antenna that will successfully reproduce a 
received impulse waveform with near-perfect fidelity in the frequency range of 10 MHz to 


10 GHz. 
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V. VALIDATION OF THE SOFTWARE AND ANALYSIS OF THE PROTOTYPE 
TEM HORN 


In order to first demonstrate the accuracy of the software and to secondly validate 
an UWB antenna design meeting certain parameters, various antenna geometries were 
developed and fully tested using the AWAS and NFLDVIS.M analysis package. A 
summary of each geometry, its purpose in the validation stage, and its performance as a 
wideband or UWB antenna is provided in the next section. Following that section, the 


prototype TEM Horn is then analyzed as an UWB antenna candidate. 


A. SIMPLE WIRE ANTENNA STRUCTURES 


To validate the accuracy of the software package, several basic two-wire antennas 
were created and tested. The specific behavior of each antenna structure presented in this 
section could be determined from known radiation characteristics of the simple two-wire 
dipole and the “vee” dipole, and from transmission line theory. The fundamental 
properties of dipole antennas are both well-known and easily demonstrated using basic 
antenna theory. As a result, they lend themselves well towards validation/verification of 
the software. 

As a point of reference for comparison between these antennas, the near-field 
values for each of the simple wire structures were calculated over a horizontal plane “cut” 
through the centerline of each antenna. The physical dimensions of the plane extended to 
two meters on either side of the antenna (along the x-axis) in one dimension, and from the 
antenna feed to a point 1-to-2 meters beyond the antenna aperture in the other dimension. 
Because of this, the perspective presented to the viewer of the time-domain movie is as if 
he/she is standing at the antenna feed and looking straight down the centerline axis of the 


antenna (the positive y-axis). 
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1. The Simple Dipole 


The dipole, as shown in Figure 20, consists of two equal length arms extending 
from the origin along the positive and negative z-axis. Dipoles are resonant structures, 
since the current wave travelling outward from the feed point to the end of the wire is 
reflected, supporting the formation of a standing wave. From basic antenna theory, it is 
well known that an ideal dipole of length ZL € A oriented along the z-axis will create a 
radiation pattern that resembles a donut shape surrounding the midpoint of the dipole, as 
illustrated in Figure 21 [2]. This pattern can be considered omnidirectional, since it is 
uniform in the xy-plane (the H-plane of this antenna). Further, it can be shown that 
whenever the length, L, of the dipole is an odd multiple of half the operating wavelength, 
pe 


D = (n=1,3,5,...), (5.1) 


the antenna will always radiate some portion of its power in the H-plane (the xy-plane for 
the antenna depicted in Figure 20). The symmetry observed in the radiation pattern of the 


dipole should then infer an expectation of symmetry in the time-domain E-field plot 


Current distribution 
Due ] = |, cos Bz 


Transmission 
line 


Dipole _- 


antenna 





Figure 20. A Half-Wave Dipole and Its Current Distribution 
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Figure 21. Radiation Pattern for an Ideal Dipole 


displayed by the software, and thus, the application of an impulse voltage to the dipole 
should result in a time-domain E-field which expands radially in the H-plane. 
For the dipole geometry shown in Figure 22, the theoretical z-component of the E- 


field can be expressed as 











c m e cam 
E, =- ponj e - 2cos(kl, ) - | (5.2) 


l 2 
and if we limit ourselves to measuring the field values along the y-axis, only, then 


Equation (5.2) simplifies to 








e eM 
E, =-j601, : - cos(kl, ) (5.3) 


s 
where 7, is the maximum amplitude of the sinusoidal current distribution on the dipole; 
k=w/c=2x/4 is the propagation constant, R= R, = R, = ae 1s the distance 
from the endpoints of the dipole to the arbitrary observation point P(0, y,0) along the y- 
axis; /, 2 L/2 1s the length of each arm of the dipole; and r = p is the radial distance 


from the dipole feed to the field observation point P [8]. Of note is the fact that the last 
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term of Equations (5.2) and (5.3) vanishes whenever L 2nA/2 for n 7 1,3, 5,... because 


cos(Kl,) 2 cos(nz /2) 20, forn- 1,3, 5,... 





Figure 22. Geometry Used in Dipole E-Field Calculations. From [8]. 


Created using the MATLAB code found in Appendix B, Figure 23 is a plot of the 
magnitude of E; as a function of both frequency and position along the y-axis. Most 
obvious is the periodic pattern of maximum and minimum field values that occurs over 
the entire frequency range at a distance of approximately 0.1 meters from the dipole feed 
point. If we solve for the wavelength, A, where the half-wave dipole condition is met, 1.e. 
where L = 4/2 , we find that for L=2.06156 meters (the total length of our dipole), a peak 
field value should occur at a wavelength of 4.12312 meters, or a frequency of 72.76 MHz. 
Not coincidentally, then, this 1s approximately the frequency where we observe our 
largest peak in the E-field plot of Figure 23. In addition, we are able to determine upon 
closer inspection of the plots that all of the successive peaks occur precisely 145.52 MHz 
apart. This frequency corresponds to the inverse of the roundtrip time it takes a pulse to 
travel from the feed of the dipole, to one of the endpoints, and back. This repetition 


frequency is equal to the resonant frequency of the dipole [2]. 
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Theoretical Near-Field IE ZG. f] for a Dipole of Length 2.0616m 


Normalized |E.(y,f)] 





1000 © y (meters) 


Frequency (MHz) 
Figure 23. The Theoretical Near Electric Field, |E,(y,f)|, of a Simple Dipole 


Returning to Figure 23, for a given fixed frequency, we note that an interference 
pattern develops as the distance from the feedpoint is increased. Figure 24 more clearly 
illustrates this interference pattern. Considering the fact that the simple dipole can be 
modeled as the superposition of three point sources (one at each end of the dipole and one at 
the feed point), it is clear that this pattern develops as a result of the constructive and 
destructive interference of the fields produced by the three sources. As an initial test of the 
visualization software, it then makes sense to construct a dipole of similar dimensions and 


compare the results to theory. 
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Figure 24. The Theoretical Interference Pattern in the Near Field of the Simple Dipole 


After creating the dipole in AWAS and calculating the frequency- and time- 
domain near-fields, Figure 25 was created by the software for the "steady-state" 
condition, 1.e. the representation of the fields after a unity amplitude sinusoidal voltage 
signal was applied to the antenna feed. Upon quick inspection, one should notice several 
"peaks" in the E-field plot. Not coincidentally, a distinct peak E-field value is at 
approximately 70 MHz, and several others appear approximately every 150 MHz after 
that. Theoretically, we expected the first peak to occur at 72.76 MHz, with each 
successive one precisely 145.52 MHz apart. So far, the fields generated by the software 
are consistent with theory. One primary difference, though, between the fields calculated 
by the Method of Moments and those derived from Equation (5.3) lies in the lack of 
distinct minima in the region 0.1 meters from the feed. This is due to AWAS's use of 


fewer discrete frequencies in solving for the fields and also the fact that the theoretical 


calculations ignore such real-life phenomena as fringe effects and assume 1deal current 


distributions along the entire length of the dipole. 
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Figure 25. Plot of E-Field Magnitude vs. Frequency and Position Along the Y-Axis 
for the Simple Dipole 


Also, when we look at this plot from the perspective shown in Figure 26, we are 
able to closely examine the interference pattern within the E-field for comparison with 
the theory illustrated by Figure 24. Over the entire frequency range, the interference 
pattern illustrated 1n Figure 26 is consistent with the theoretical pattern discussed earlier. 
All that remains, then, is a closer examination of the results generated by both AWAS 


and the theoretical model at specific frequencies. 
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Figure 26. Plot of E-Field Magnitude as a Function of Distance vs. Frequency 


Using the MATLAB code of Appendix B, a direct comparison of AWAS versus 
theory was made at frequencies of 36.38 MHz, 72.76 MHz, and 218.281 MHz. For the 
given dipole length of 2.06156m, these correspond to the frequencies where the dipole is 
electrically a quarter-wavelength, a half-wavelength, and three-halves of a wavelength 
long. The following figures are comparisons of the magnitude of £,, its real component, 
Re{E,}, and its imaginary component, Im{£E,}, as functions of position along the y-axis 

In Figure 27, the AWAS-derived solution for the E-field is almost exact in 
comparison with the theoretical values. This indicates that the theoretical approximation 
of the near-fields is obviously a good one for an electrically-short dipole. It is also 
indicative of the accuracy of AWAS in calculating the near E-fields. Additionally, it is 
evident that the two solutions converge as the observation point is moved further from the 
antenna. Figures 28 and 29 are the comparisons of the real and imaginary components of 
the field, respectively. As expected, they, too, indicate strong correlation between 


AW AS and theory. 
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Near-Field |Ez| for L= 2.0616 Dipole at z= 0 m & E 36.38MHz 
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Figure 27. Comparison Between AWAS vs. Theoretical |E;| for a A/4 Dipole 
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Figure 28. Comparison Between AWAS vs. Theoretical Re{E,} for a 4/4 Dipole 
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Figure 29. Comparison Between AWAS vs. Theoretical Im{E,} for a ?/4 Dipole 


For the half-wave dipole case, Figures 30 through 32 again clearly show good 
correlation between the AWAS model and theory. As the electrical length of the dipole 
increases, though, the margin of error between the two also increases. As noted 
previously, it is the theoretical solution which is losing accuracy, and not the Moment 
Method solution. The margin of error is even more significant in the 32/2 dipole case, as 


shown in Figures 33 through 35; however, most of the error is merely amplitude error, 


with the general “shape” of the theoretical and AWAS-calculated fields remaining well- 
correlated. 
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Near-Field |Ez| for L= 2.0616 Dipole at z= 0 m & f= 72.76MHz 
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Figure 30. Comparison Between AWAS vs. Theoretical |E,| for a 1/2 Dipole 
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Figure 31. Comparison Between AWAS vs. Theoretical Re{E,} for a 1/2 Dipole 
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Figure 32. Comparison Between AWAS vs. Theoretical Im{E,} for a 1/2 Dipole 
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Figure 33. Comparison Between AWAS vs. Theoretical |E;| for a 34/2 Dipole 
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Near-Field Re(Ez) for L= 2.0616 Dipole at z- 0 m & f= 218.281MHz 
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Figure 34. Comparison Between AWAS vs. Theoretical Re{E,} for a 3A/2 Dipole 
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Figure 35. Comparison Between AWAS vs. Theoretical Im(E;) for a 3A/2 Dipole 
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With the above observations validated by our theoretical model, we now turn our 
attention to the time-domain representation of the dipole's near E-field. Specifically, 
since we are viewing the propagation in the H-plane, we look for a radially symmetrical 
propagation pattern due to our Gaussian impulse. In lay terms, what we expect to see is 
an image similar to that seen by dropping a rock into a still pond: an ever-expanding 
circular propagation of a wave. As seen in Figure 36, this 1s most certainly the case. One 
thing to note about the figure 1s that the additional wave being generated behind the first 
one is not due to an additional input impulse, but rather, due to current reflections along 
the length of the antenna which partially launch into free space on each round trip along 
the structure. Note, too, that the period between “waves” is approximately 2 meters. A 
repetition frequency of 145.521 MHz means a period between reflections of 6.872 nsec. 
At the speed of light, an RF pulse travels 2.0616 meters in 6.872 nsec, which coincides 
exactly with the distance between the original Gaussian impulse and the first reflection in 


Figure 36. 
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Figure 36. A Snapshot of the Time-Domain E-Field of a Simple Dipole 


Transmitter at t=13.5nsec 
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To graphically illustrate the Lorentz Reciprocity Theorem discussed in Chapter II, 
we can look again at the time-domain E-field of the simple dipole as a receiving antenna. 
For this case, we use a double-Gaussian impulse as the incident wave, propagating in the 
—y direction (as shown in Figure 37). This incident field pulse shape has an amplitude of 
] V/m and contains no dc (a single Gaussian pulse does). The incident wave impacts the 
dipole and induces a current along the antenna structure. As the current proceeds along 
the dipole, it produces a reflected field (Figure 38) which is the same as in the 
transmitting dipole case (Figure 36), with the exception of a scale factor. In addition, the 
resonance of the dipole over such a large frequency range guarantees that we get the same 
current reflections, which partially launch into free space on each roundtrip along the 
structure. Thus, either the transmitting antenna or receiving antenna can be used 
interchangeably in any discussion of the antenna's near fields. 
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Figure 37. Incident E-Field Generated by a Double-Gaussian Impulse at t=15nsec 
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Figure 38. The Time-Domain E-Field of a Simple Dipole Receiver at t=20.5nsec 


Regardless of whether the simple dipole is analyzed as a transmitting antenna or 
directly as a receiving antenna, it is clear that both methods of analysis accurately 
determine that the simple dipole 1s a poor UWB receiving antenna. In the frequency- 
domain depictions of the dipole’s near E-field (Figures 25 and 26), the magnitude plot of 
the near E-field displays resonant peaks at a number of frequencies. These resonant 
modes will cause amplitude distortion of the input impulse signal. In addition, the time- 
domain depictions (Figures 36 and 37) clearly indicate the presence of current reflections 
throughout the dipole. These reflections are indicative of a nonlinear phase change, since 
it is the lower freequency components of the impulse that we view “ringing” back and 
forth along the dipole. Thus, the simple dipole causes phase distortion of the input 
impulse, too, which makes it an unsatisafactory choice as an effective UWB antenna. 

Using the simple dipole model, it has been shown that the software does, in fact, 
accurately display the electromagnetic phenomena associated with the application of 


input pulses to a dipole antenna. To ensure the software is robust and able to handle 
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complex antennas, we now move on to test 1t with some "special case" dipole antennas 


and modified two-wire transmission lines. 


2. The Vee Dipole 


The Vee Dipole is best thought of as a simple dipole whose ends are bent towards 
the y-axis to form an angle, y, as shown in Figure 39. It can also be thought of as an 
open-ended transmission line whose ends are bent toward the z-axis. Because of its 
unique shape, the directivity of such an antenna will be greater than that for a straight 
dipole, since the radiation pattern 1s merely the superposition of that from each straight- 
wire section [2]. Thus, the Vee Dipole is able to focus substantially more power down its 
boresight axis, and thereby reduce the sidelobe field levels. For the Vee Dipole shown, 
the corresponding directivity is 

D 2 294(h/A) - 115, (5.4) 
where Ah is the length of each arm of the dipole (^ 2 L/2 21.03078m for our example) 


Example of a Vee Dipole 
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Figure 39. Example of a Vee Dipole 


[2]. This equation is valid only over the frequency range 145.5 MHz to 873 MHz for our 
specific antenna. Regardless, it 1s expected that over the entire tested frequency range, 


more energy will be focused along the centerline axis than elsewhere. 
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To verify the above, a plot of the frequency-domain representation of the E-field 
along the centerline axis (y-axis) of the Vee Dipole was rendered. It is replicated here as 
Figure 40, and as a point of reference within the plot, the antenna aperture is at the y=1m 
point on the plot. What should be immediately noticed is that the maximum value of the 
field strength is found at the frequency meeting the half-wave dipole condition, just as it 
was for the simple dipole. The field strength in this case, however, is significantly larger. 
One reason is that within the 0-to-1m distance along the y-axis, the field is somewhat 
constrained by the physical structure of the antenna, and as such, the field is not able to 
expand as rapidly as in free space. The other is, of course, the directivity the antenna 
provides as discussed earlier. In addition, the resonant modes of this antenna are identical 


to those of the simple dipole tested previously. This 1s because the antenna arms of the 


tested Vee Dipole are exactly the same length as those of the simple dipole. 
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Figure 40. Plot of E-Field Magnitude vs. Frequency and Distance along the 
Centerline Axis for the Vee Dipole 


In Figure 41, we are able to better see the changes to the E-field for each 


frequency as a function of distance along the centerline axis. In this plot, it 1s clearly 
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evident that the field strength at 72.76 MHz falls off less quickly than for any other 
frequency along the entire length of the antenna. Also, the field strength at this frequency 
outside of the aperture (in free space) is, indeed, significantly larger than the simple 
dipole case. Of further note is that for the other frequencies, the field strength within the 
bounds of the antenna structure takes on an interference pattern with distance that is 
obviously a function of wavelength. Specifically, we know from transmission line theory 
that this 1s similar to the voltage standing wave of the transmission line and that the 
conditions (a maximum or minimum) repeat themselves every half wavelength. As an 
example, for a wavelength of 0.317m corresponding to the 946 MHz frequency, we note a 
distinct maximum at 0.26m along the y-axis. The next maximum then occurs at 0.42m, 


or approximately 0.16m (a half wavelength) after the first. 
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In addition, for those frequencies where at least one half-cycle (a maximum and a 
minimum) occurs within the antenna structure (between 0-to-1m along the y-axis), a large 
E-field component also denotes resonance (a reflection of the pulse at the aperture). The 


larger the value of the E-field, as 1s the case at 218MHz, the more energy is reflected back 
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from the aperture. In contrast, the lower the value of the E-field, as is the case at 
946MHz, the greater coupling efficiency should be observed in the time-domain. 

In order to determine whether the antenna’s directivity 1s specifically as theorized 
above, a time-domain plot is required, since that is where the field strength in more than 
one direction can be viewed. Figure 42 is a rendering of the time-domain E-field taken 
]3.5nsec after application of the Gaussian impulse wave. Clearly, the H-plane 
propagation pattern 1s of semicircular shape in free space with the greatest field strength 
displayed along the +y direction. Thus, the directivity enforced by this antenna is 
confirmed. Additionally, one should note the large disturbance near the origin and along 
the 0-to-1m point on the y-axis. This 1s caused by reflections generated by waves unable 
to immediately couple into free space at the very low frequencies. As they ring back and 
forth along the antenna, though, they eventually do couple and launch what appear to be 


additional impulse waves. 
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Figure 42. A Snapshot of the Time-Domain E-Field of the Vee Dipole 
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Like the simple dipole tested earlier, the vee Dipole suffers from the same 
amplitude and phase distortion problems. Thus, it, too is an unsatisfactory candidate for 


the role of an UWB receiving antenna. 


3. The Vee Dipole with a Parallel Extension 


After analyzing the Vee Dipole, a parallel two-wire transmission line was added 
to the antenna at the aperture, as shown in Figure 43. This structure was chosen next to 
evaluate the accuracy of the reflection coefficients determined by the software at each end 
of the transmission line model. Also, the change ın the directivity and frequency 
response of the Vee Dipole due to the extension was also investigated. The first 
difference that should be observed between this antenna and the first two dipoles is the 
shift in the wavelength meeting the L=A/2 condition due to the increased length of the 
extension. Using the same equations as in Section 5.A.1, the frequency that meets the 


half-wave dipole condition for this antenna 1s approximately 37 MHz. 


Example of a Yee Dipole with a Parallel-Wire Extension 
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Figure 43. Example of a Vee Dipole with 


Parallel Extension 


It can be shown that the current distribution on a two-wire transmission line will 


look very similar to the ideal dipole’s at certain wavelengths [2]. Because of this, the 
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time-domain view of the radiation pattern and directivity of this antenna should not differ 
greatly from the Vee Dipole's. Since we are again viewing the magnitude of the E-field 
in the H-plane, the wave front should appear less circular than for the Vee Dipole case, 
and begin to approach a flat planar front because of the extension. The reflection 
coefficients at the feed and aperture ends of the antenna are similar to those of the Vee 
Dipole, but because of the additional length, they should be much more discernible in the 
time-domain plot. Reference [9] predicts a reflection coefficient, L, of -1 at the feed and 


I =+] at the aperture, which, when multiplied by the value of the incident wave, will 


provide the value of the reflected wave. 


Looking first at the frequency-domain representation of the E-field in Figure 44, it 
is clear that the E-field at 37 MHz displays better strength than any other frequency 
through the entire length of the antenna, as was predicted. The resonant frequency of 74 


MHz is also clearly evident in the periodicity of the E-field maxima at a distance of 


y=0.1m from the feed. 
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Figure 44. The Frequency-Domain E-Field for the Vee Dipole with Parallel 


Extension 
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Additionally, Figure 45 confirms the periodicity of the voltage standing wave for 
each frequency, just as was seen for the Vee Dipole in Figure 41. Also, the coupling 
efficiency for every frequency above the limit defined by f 2 c/L , where c is the speed 
of light, can also be determined from the magnitude of the E-field within the limits of the 


antenna structure. 
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Figure 45. A Different Perspective of the Plot Shown in Figure 44 

The time-domain representation of the E-field for this antenna, Figure 46, clearly 
shows the initial impulse propagating out of the antenna into free space. Like the Vee 
Dipole, this, too, has a more directive nature than a simple dipole; however, in contrast to 
the earlier antennas, this one’s propagating wavefront is not at all semicircular. It is clear 
that the wavefront begins to flatten out and approach a planar front as it is forced to 
propagate further down a transmission line. Also, though not directly evident in this 
frame, but exceptionally clear in the movie non of the time-domain fields, the 
reflection coefficients at each end of the antenna do. indeed, match the theoretical 


predictions. 
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In order to eliminate, or at least reduce, the amount of energy reflected back to the 
feed from the aperture, a resistance would need to be added along the parallel section of 
this antenna. The resistance would serve to convert the energy from the uncoupled 
resonant frequency waves into heat, and thereby attenuate the reflected current each time 
it would pass through the resistance. This is done in the exact same manner as was 
discussed for the TEM Horn in Chapter II. In order to test this and additional 
transmission line phenomenology with the software, one final antenna structure was 
analyzed. 


4. The Vee Dipole with a Parallel Extension Terminated at a 1 KO Load 


The final two-wire structure tested utilizing the analysis software was the Vee 
Dipole with the parallel extension now terminated at a 1kQ load. As noted in Figure 47, 
this antenna differs from that last considered in only the addition of the vertical wire load 
segment at the aperture. The purpose of this modification was to ensure that the software 
would verify a change in both the reflection coefficient at the aperture and the attenuation 


of the reflected energy caused by uncoupled low frequency waves. 
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Vee Dipole with Parallel Extension Terminated at a 1k(2 Load 
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Figure 47. Vee Dipole with Extension 
Terminated at a 1k0 Load 


It can be shown that for a transmission line of this type, the reflection coefficient 
at the aperture end will approach -1 for those frequencies that do not produce coupling 
[9]. In addition, the fact that the load is along a wire segment that 1s 0.5 meters long in 
the vertical axis, means the loaded end should exhibit some radiation properties separate 
from the rest of the transmission line. 

Looking at Figure 48, the frequency-domain representation of the E-field for this 
particular antenna structure reveals a flatter frequency response of the antenna across all 
frequencies. This is very different from the “spikey” nature of the frequency response 
seen in the past antennas’ plots (Figures 40 and 44), and implies the load is, in fact, 
attenuating the uncoupled resonant frequency waves. There are certainly still peaks along 
the frequency band; however, the nulls are not as deep and the overall magnitude of the 
E-field is more constant across the entire band. Throughout the antenna length, the fields 
observed at the lower frequencies have been attenuated, and now more closely match the 
values for the frequencies able to couple and launch into free space. This does not mean 
that the resistive load will allow more frequencies to couple than for the unloaded 


antenna. Rather, it is merely attenuating the reflected energy at those uncoupled 
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frequencies. If any additional frequencies are able to couple more efficiently, it is most 
likely due to the 0.5 meter segment that was added to support the load. As stated earlier, 
this additional segment, although loaded, will not only add length to the original antenna, 
thereby improving the low-frequency characteristics, but also serve as an additional 
focused radiation point for the antenna. This additional radiation point by no means aids 
the radiation performance of the antenna; rather, it merely serves as another source point 


on the antenna structure. 
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Figure 48. Frequency-Domain E-Field of the Extended Vee Dipole Terminated at a 
1k0 Load 


The final area of interest is the 2-meter point in Figure 48. Specifically, this is 
where the aperture was replaced by the vertical loaded segment. Across all frequencies, 
there appears a ridge of slightly higher field strength. This is undoubtedly due to the 
radiant nature of the loaded segment. On one hand it is dissipating low-frequency energy, 
but on the other, it is most certainly serving as a source point. If this is true, then there 


should be evidence of this in the time-domain plots. 
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Figure 49 is a frame selected from the time-domain movie of this antenna's E- 
field. Chosen at a point after the initial pulse has reached the loaded end and is reflecting 
back to the feed point, this plot graphically depicts the predicted propagation pattern. 
Specifically, it should be noted that the wave propagating in the positive y-direction 
displays the same characteristics of the initial wave propagating through the antenna of 
Section 5.A.3, as depicted in Figure 46. In addition, the reflected wave and part of the 
leading edge of the initial wavefront are displaying propagation characteristics similar to 
those of the simple dipole. Also, the portion of the reflected wave propagating along the 
wire structure is definitely at a lower magnitude than the reflected waves from either of 
the last two antenna models. In hindsight, then, this antenna has functioned in many 
ways as the superposition of two different antennas upon each other. Although this was 


not directly predicted through theory, the end result was, at a minimum, enlightening. 
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Figure 49. Time-Domain E-Field for the Extended Vee Dipole Terminated at a 
1kQ Load 

It must be stated that the software not only was able to visually confirm that 

behavior which was theoretically predicted, but also provided heuristic insight. Overall. 


it is easy to see that the graphical analysis, while perhaps not as specific as a methodical 
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mathematical process, does, 1ndeed, allow the user to display certain antenna performance 
characteristics, and further present an accurate picture of the near-fields for further study 


and in-depth analysis. 


B. TEM HORN PROTOTYPE DESIGN 


Having validated the software’s capabilities as an analysis tool, we now turn our 
attention to an antenna for which the software was initially designed to support, the TEM 
Horn. Figure 50 is the wire model used to emulate the prototype TEM Horn. Its 
dimensions are 1.8288m (6ft) long, 0.1524m (6in) high, and 0.1026m (4in) deep. As 
discussed in Chapter II, resistive loading has been added along the parallel extensions. 
Resistive loads are placed between each of the dark nodes on the rıght half of the antenna. 
Figure 50 is a wire-based analog of the continuous surface structure shown in Figure 1. 
As such, it offers current paths which are the same as the actual antenna. Extensive 
validations [5] have shown that under the conditions observed here, such wire structures 


provide accurate emulations of solid surface antennas. 
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Figure 50. 3-D and Overhead Views of the Prototype TEM Horn 


76 


The goal of the analysis for this antenna, then, is first to develop an understanding 
of the E-field generated within the near field by the TEM Horn, and second, to 
demonstrate that the prototype TEM Horn either meets or fails to meet the design criteria 
of an UWB impulse receiving antenna operating in the 10MHz-10GHz bandwidth. To do 
so, we will analyze the frequency-domain E-field plot and garner whatever insight we 
may prior to investigating the time-domain response of the antenna. The analysis of the 
E-fields in the frequency and time domains should provide us with enough visual 
evidence of the antenna's success in meeting the design criteria set forth earlier. 

The first concern for the TEM Horn is a need for a relatively flat response across 
all frequencies. For this antenna, an initial frequency range of 0 to 2 GHz was analyzed. 
Figure 51 depicts the TEM Horn's E-field as a function of frequency and distance along 
the antenna's centerline axis, and Figure 52 provides the E-field at the antenna feed point, 
only. What should be evident immediately 1s that the response of the antenna, while not 
perfectly straight across the top of the plot, provides very little deviation from one 


frequency to the next. The "ridge" present at a very low frequency may indicate a 
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Figure 51. Frequency-Domain Plot of the E-Field for the Prototype TEM Horn 
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problem in the design; however, this plot does not provide a small enough scale to 
accurately identify a specific frequency. Although we expect to have poor response 
below the lower design limit of 10 MHz, we must ensure the response from 10-100 MHz 
is not causing fluctuations in the field, thereby causing the output to be distorted and lack 


fidelity. Thus, the 0-100 MHz Teque range will be EL in more Ma later. 
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Figure 52. Frequency Response of the TEM Horn at the Feed Point 


Looking at Figure 53, we now are able to view the response of the antenna vs. the 
distance along the centerline axis. Clearly, the field strength along the y-axis weakens 
rapidly yet smoothly for all frequencies above 200 MHz. Below this point, we notice a 
“bulge” beginning to form at the end of the antenna for the lower frequencies. This 
bulge, circled in white, is an indication that some of this low frequency energy is not 
completely coupled and launched into free space. As a result, the antenna may have 
degraded performance (i.e., less-than-desirable fidelity) at the very low end. To 


investigate this further, a low-frequency analysis is a must. 
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Figure 53. Frequency-Domain Plot of the TEM Horn's E-Field with Area of 
Interest Highlighted 


Prior to beginning the low frequency analysis, though, a study of the time-domain 
waveform was conducted to verify the assumptions made above. Figures 54 and 55 are 
frames taken from the time-domain propagation movie. At 9nsec after application of the 
Gaussian impulse, the waveform in Figure 54 is maintaining excellent shape as it 


propagates through the antenna, with no noticeable deformation occurring. 
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Figure 54. Time-Domain E-Field for the Prototype TEM Horn 

At 10nsec after application of the Gaussian impulse, the wavefront reaches the 
aperture, and as seen in Figure 55, most of the energy launches into free space. This 
wave, a double Gaussian, 1s approximately 2.5m from the antenna feed point. The low 
frequency energy that did not couple is then reflected back through the resistive elements 
and dissipated as heat. The large bulge between 1.25m and 1.8m is most of the reflected 
wave being "trapped" in the resistive load and dissipated. The strong field near the feed 
(at y-0.1m) is of no concern. While it may appear to be significant in this plot, one must 
remember that the magnitude of the E-field is dependent upon the area within which it is 
concentrated. Beyond the aperture, the area is growing rapidly, since the wave is in free 
space; whereas, at the feed, the field is confined within the few centimeters separating the 
top and bottom plates. Thus, the E-field at the antenna feed point, while small, appears 


amplified because it is so confined. 
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Figure 55. Time-Domain E-Field for the Prototype TEM Horn 


a, 
NE deae fio Ti tu d 


The reflection at the aperture is still cause for concern, since it validates the 
prediction made after viewing the antenna's frequency response. To determine the exact 
frequencies that are causing this, the antenna was next analyzed over the frequency range 
0 to 100 MHz. Figures 56 through 58 are plots of the antenna's frequency response over 
this range. What is immediately noticeable is the ridge of higher-than-normal field 
strength located at approximately 40 MHz, which is well within the antenna’s required 
bandwidth. Unfortunately, neither the frequency- nor time-domain plots give any 
indication as to the cause of this anomaly in the frequency response. Specifically, the 
answer can only be determined through in-depth analysis of the antenna’s key parameters 
at the frequency in question. While the software in this case may lack the depth to 
determine the cause for poor frequency response of an antenna, it does have the range and 
ability to at least display possible areas of concern for the designer. 

In Chapter II, it was determined the sensitivity of the prototype TEM horn could 
suffer by as much as 3 dB near 26 MHz. From Figures 56 through 58, it is clear that the 
frequency response of the prototype TEM horn experiences a maximum deviation of only 


1.5 dB in the low frequency range, which will have a negligible impact on performance. 
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56. Frequency-Domain E-Field for the TEM Horn at Low Frequenc 
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Figure 57. Prototype TEM Horn E-Field Showing Possible Resonant Frequency 
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Figure 58. Close-up of Possible Resonant Frequency 


As a final analysis step, the TEM horn was also set up as a receiving antenna with 
a double-Gaussian impulse functioning as the incident field. Recall that reciprocity 
dictates that the radiation pattern of an antenna is the same whether it is used as a 
transmitting antenna or receiving antenna. Although the excitation for our transmitting 
and receiving cases is not exactly the same, it makes sense that the general behavior of 
the TEM horn's near fields will be similar. In the ideal case, however, further validation 
of the software's accuracy could then be accomplished by matching the two near fields 
derived under the two different circumstances if the excitation in both cases were the 
same (1.e., the transfer impedances were equal). 

In Figure 59, the double-Gaussian incident impulse is seen propagating towards 
the TEM horn. At t=13.5nsec, the incident impulse field impacts the antenna feed, and a 
reflected wave begins to form and propagate in the +y direction along the antenna’s 
centerline axis. Figure 60 clearly illustrates the formation of the reflected wave at 


t=15nsec. Finally, at t=16.5nsec, the reflected wavefront reaches the aperture, and as 


seen in Figure 61, the resultant transmitted and reflected fields behave and look the same 


as they did in the transmitting antenna case (Figure 55). 
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Figure 60. Time-Domain Reflected E-Field within the TEM Horn Receiver 
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Figure 61. Resultant Transmitted and Reflected Fields at the TEM Horn Aperture 





VI. CONCLUSIONS 


The focus of effort of this research was in the creation of a high-level graphics 
package that would primarily aid in the design and analysis of an UWB antenna and 
secondarily function as an instructional aid in the understanding of basic radar and 
transmission line theory. In an effort to reach these goals, the MATLAB function 
NFLDVIS.M was written and designed to work hand-in-hand with the data output files 
created through AWAS. Specifically, with a general knowledge of both AWAS and 
NFLDVIS.M, a user could build a wire model of an antenna or transmission line and then 
actually view the magnitude of the associated Near Fields as either a function of frequency 
and distance, or time and distance. 

As a design tool, the software was demonstrated to perform well with both simple, 
classical antenna designs, as well as state-of-the-art TEM Horns. Although the designer 
does require the use of two separate sets of software, the synchronization between the two 
should make that issue somewhat transparent. VFLDVIS.M is designed to be user-friendly, 
and with the designer only required to use a keyboard for typing in the title for his plots, this 
goal has been met. 

As for the TEM Horn itself, the software indicated the prototype antenna design 
may have problems at the low frequency end of the spectrum. The one anomaly present at 
approximately 40 MHz may be indicative of a natural resonance in the antenna, but since 
the antenna tested in AWAS (Figure 50) is not exact in its representation of the resistive 
loading design, it is probable the reflection noted in Figures 56 through 58 will be 
attenuated more aggressively by the logarithmic loading scheme planned. In addition, 
recall that in Chapter IJ, it was theorized that the antenna may suffer a sensitivity 
degradation by as much as 3 dB near 26 MHz and 6 dB near 13 MHz; however, there was 
no indication of such in the frequency response curves created by the software. On the 
contrary, the fluctuation in the magnitude of the frequency response of the antenna never 


exceeded 1.5 dB over the entire tested range. To absolutely guarantee against any 
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degradation in the sensitivity, increasing the length, L, of the TEM hom to 4.775m will 
keep the lower frequency limit outside of the intended spectrum of the TEM horn. 

It was further theorized in Chapter II that the TM,, mode would be excited at some 
frequency above 6.225 GHz, thus causing distortion of the output waveform. 
Unfortunately, due to the limitations of AWAS discussed in Chapter III, the exact frequency 
where the TM,, mode was excited could not be verified. Assuming that were the case, 
though, the TEM horn would require a redesign to guarantee distortionless transmission 


throughout the entire frequency band. Specifically, if the antenna’s total aperture height 


was reduced to 94 mm, then fq, will definitely occur outside the frequency band of 


interest. 

At no time did the software package ever produce faulty data or display unjustifiable 
renderings, except for the occasional “garbage-in, garbage-out” operator errors. In other 
words, every plot or time-domain movie the software generated was completely in line with 
the theoretical expectations of the particular antenna or transmission line in question. In the 
long run, though, the true future for the NFLDVIS.M function may lie in its ability to 
provide a detailed glimpse into what most people never see. Beginning with the 
verification of the software's accuracy, 1t became increasingly evident that visualization of 
the near fields in either the frequency- or time-domain proved an invaluable educational 
tool. For even though one may hear or read about concepts, or take measurements to 
confirm theory, the actual ability to see not only the effects, but also the cause adds a new, 
logical dimension to the learning process. The old adage, “a picture 1s worth a thousand 
words” rings true. But even this software cannot tell you everything. As seen in Section 
5.B, it 1s capable of indicating potential problem areas in a design, but it can never tell you 


the “why”. That answer is left for the designer. 
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APPENDIX A. NEAR-FIELD VISUALIZATION SOFTWARE (NFLDVIS.M) 


function nfldvis(action) 
VoNFLDVIS This Graphieal-User-Interface (GUI) provides the antenna designer 


% with a tool to view the electric and magnetic fields generated 

% within and in elose proximity to an Ultra- Wideband Impulse Antenna. 
% This program uses AWAS-generated near-fteld (*.nfl) and geometry 
% (*.geo) files as input. and provides graphieal output of the 

% near-fields in both the frequeney and time domains. In the 

% frequeney domain. the designer ts able to view the fields in 3-D 

% as a funetion of frequeney along the antenna axis. In the time 

Yo domain. a movie is ercated. so that the designer may view the 

% antenna as it either transmits or reeeives an impulse waveform. 

% For the ease of the antenna operating as a transmitter. the user 

% may view the wave as it propagates through the antenna and ts 

% "Jaunehed" into free space and/or refleeted back to the feed point. 
% For the ease of the antenna operating as a receiver. the user mas 

% view the wave as it propagates through frec spaee. impacts the 

% antenna aperture. and refleets off the antenna or continues to 

% propagate within the antenna strueture. As an additional option. 

% the designer may simply plavbaek movies that were ercated earlier 
% and saved as .mat files. 

% 

% Follow the prompts provided in the text window at the top of the 

% figure to step through the analysis process. The "slides" at the 

% bottom of the figure window allow the experteneed user to control 
% the display of the data. Default settings are already built into 

% the program. The program also prompts the user when input from a 
% slide. popup menu. or text box ts required. 

% 

% Display of the frequenes - domain ficlds 1s eontrolled by the "Plot" 
% menu on the menubar at the top of the window, while display of the 
% time-domain representations is eontrolled bv the "Movie" pushbutton 
% on the right side of the window. 

% 

% A full demonstration of the software ts provided merely by pressing 
% the "Demo" pushbutton on the right side of the window. 

% 

% Although a "Print" menu is available on the menubar at the top of 
Yo the figure window, the only menu seleetion whieh ts currently 

% operable ts the "Save..." option. whieh allows the user to save the 
Yo entire figure as a .TIF file. The other options are eurrently 

% reserved tor future use. When operational, they will allow for 

% printing of either the matn axes separately or both the main axes 

Yo and the eolorbar together. Since MATLAB does not vet offer printing 
% of seleet "frames" (axes) within a figure window. these two options 


9] 


% have been deactivated. awaiting a future release of MATLAB. Ifa 


% print of the entire GU] (ic, buttons. graph, colorbar, and all other 
% uicontrols) is desired, merely use the normal "Print" submenu found 
% under the "File" menu at the top left corner of the figure window. 


% Written & designed by Dave Adamiak 
% Copyright (c) 1997 
% Last mod: 10/10/98 


% Set the required initial conditions & create the GUI: 
if nargin<l, 

action='initialize’; 
end; 


switch(action) 
0/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/00/09/09/09/00/09/09/09/09/09/69/69/09/69/09/09/09/99/99/0 
0/99/09/09/09/09/09/09/o9/o 9/09/00 9/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/69/69/69/09/69/69/o 
case 'initialize' 
0/0 /69/0 0/6 /09/69/69/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/9/09/69/09/09/09/69/69/69/09/09/69/69/09/69/99/99/6 
O/69/69/o9/69/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/69/o 

oldFigNumber=watchon; 


figNumber=figure( 'Color',[0.6 0.6 0.6], ... 
'Name','Transient Field Visualization for Ultra-Wideband Antenna Design by Dave 
Adamiak, ... 
'NumberTitle','off, ... 
'Units';'normalized',Position',[0 0.036 1 0.9], ... 
Narben 
'BackineStore','off'); 
axHndl-axes('Color,[1 1 1], 'Position', [0.075 0.2 0.675 0.6], ... 
'Visible','off,"Tag','Axes1*); 


text(0,0. 1, Press the "Start" button to begin the analysis, ... 
'HorizontalAlignment’,'center’); 

text(0,0,'Press the "Demo" button to see a demonstration of the software;', ... 
'HorizontalAlignment','center'); 

text(0,-0.1,'Press the "Movie" button to replay a saved movie or create a new one. ... 
'HorizontalAlignment','center'); 

text(0,-0.2,'Or click on the "Plot" menubar to create various plots using data from your last 

session.', ... 

'HorizontalAlignment','center); 

text(0.-0.4, All commands will be issued in the text window at the top of this screen." ... 
'HorızontalAlignment','center'); 

axis([-1 1 -1 1]) 
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% The INFO button 
callbackStr = 'nfldvis("info")': 
infoHndl = uicontrol('Units' 'normalized', ... 
'Position'; [0.795 0.8168498168498168 0.08 0.0787545787545 7875]. ... 
‘String’, Info’, 'Callback',callbackStr, 'Tag’,’Pushbutton!'); 


% The START button 
callbackStr 7 'nfldvis("start")'; 
startHndl = uicontrol('Units','normalized', ... 
'Position',[0.89125 0.8168498168498168 0.08 0.07875457875457875], ... 
'String','Start', Callback',callbackStr, 'Interruptible','on', ... 
"Tag, 'Pushbutton2); 


a —— — —À — —— — — — — —— — — — — — — 


% The DEMO button 
callbackStr = 'nfldvis("demo")'; 
demoHndl = uicontrol('Units','normalızed), ... 
'Position',[0.89125 0.7142857142857143 0.08 0.07875457875457875], ... 
"String", Demo’, 'Callback',callbackStr, '"Tag','Pushbutton3'); 


% The MOVIE button 
callbackStr 7 nfldvis("movie")'; 
movieHndl = uicontrol('Units' 'normalized", ... 
'Position',[0.89125 0.6105006105006106 0.08 0.07875457875457875], ... 
'String‘,'Movie', 'Callback',callbackStr, 'Tag',’Pushbutton4'); 


% The CONTINUE button 
continuHndl = uicontrol(‘Units', normalized’, ... 


'Position', [0.89125 0.5073260073260073 0.08 0.07875457875457875], ... 


'String','Continue', 'Enable','off, "Tag',Pushbutton5!); 
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% The RESET button 
callbackStr 7 'nfldvis("reset")'; 
resetHndl = uicontroi('Units', normalized’, ... 
"Position',[0.89125 0.4047619047619048 0.08 0.07875457875457875], ... 
'String‘,'Reset', 'Callback',callbackStr, ... 
‘Tag’, Pushbutton6'); 
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9o | he CLOSE button 
callbackStr = 'close(gcf)’; 
closeHndl = uicontrol(‘Units', normalized’, ... 


Position',[0.89125 0.3003663003663004 0.08 0.07875457875457875], ... 
'String' 'Close', 'Callback',callbackStr, "Tag','Pushbutton 7’); 


—————————————————————ÓÁ——————————————————— 


96 The Instructional Text Box 
edittextHndl = uicontrol('Units' 'normalized', 'BackeroundColor',[1 1 1], ... 
'Position', [0.09875 0.9065934065934066 0.60375 0.0695970695970696], ... 
FontSize 10] Max PI Styles edits sini se 
'Look Here for All Commands & Prompts", 'Tag';'EditText1'); 


% The Near-Field Group popup menu 
popupStrl = ('Near-Field Group”); 
grppopupHndl - uicontrol('Units','normalized', ... 
'BackgroundColor',[0.75294] 0.75294] 0.752941], ... 
'Position',[0.09875 0.8589743589743589 0.12625 0.0347985347985348], ... 
‘String',popupStr], 'Enable','off‘, 'Style','popupmenu', 'Tag’, ... 
'PopupMenul', 'Value',1); 


% The Total Frames popup menu 
popupStr2 = (‘Total Frames’); 
frampopupHndl - uicontrol('Units','normalized'", ... 
'BackgroundColor',[0.752941 0.75294] 0.752941], ... 
'Position', [0.3375 0.8589743589743589 0.12625 0.0347985347985348], ... 
'String',popupStr2, 'Enable','off', 'Style','popupmenu', 'Tag., ... 
'PopupMenu?', 'Value',1); 


% The Frames-Per-Second popup menu 

for n=1:5, fpsstr(n,:)=sprintf(% 13.0f,2*n); end 

popupStr3 = char('Frames-per-sec',fpsstr); 

fpspopupHndl 7 uicontrol('Units','normalized', ... 

'BackgroundColor',[0.752941 0.75294] 0.752941]. ... 

'Position', [0.57625 0.8589743589743589 0.12625 0.0347985347985348)], ... 
'String',popupStr3, 'Enable’,'off, Style’, popupmenu’, 'Tag’, ... 
‘PopupMenu3’, 'Value',1); 


% The Z-Axis Scaling slider & its associated text boxes 
zscalslidHndl 7 uicontrol('Units','normalized', ... 
'BackgroundColor',[0.752941 0.752941 0.752941], ... 
'Position',[0.1 0.06227106227106227 0.1475 0.03663003663003663], ... 
"Style" slider”, ‘Enable’, off, ‘Tag’,'Slider1’); 
zscaltextHndll — uicontrol('Units','normalized', ... 
'BackgroundColor' [0.75294] 0.75294] 0.752941], ... 
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'Position', [0.11875 0.102564103 0.11 0.024420024], ... 
'String‘, Z- Axis Scaling, 'Stvle';'text', "Tag','StaticText1); 
zscaltextHndl2 = uicontrol(‘Units', normalized’, ... 
'BackgroundColor',[0.752941 0.752941 0.752941], ... 
'Position',[0.116666666667 0.0347985347985348 0.0175 0.02564102564102564], ... 
tine 0. Style text. las Static lext2 ) 
zscaltextHndl3 = uicontrol(Units' 'normalized', ... 
'BackgroundColor',[0.752941 0.752941 0.752941], ... 
'Position',[0.215 0.0347985347985348 0.0175 0.02564102564102564], ... 
Steine (style text. Tag. Static Texts); 


% The Azımuth View Editable Text Box 

azıeditHndl = uicontrol('Units''normalized', 'BackeroundColor',[1 1 1], ... 
'Position',[0.3266666666666667 0.06227106227106227 0.05 0.03663003663003663], ... 
Stan, >70. Emable,oif, SwileJedir, Tac, Editlext2 ), 

azitextHndl - uicontrol( Units','normalized', ... 
'BackgroundColor’,[0.752941 0.752941 0.752941], ... 
‘Position'’,[0.3229 166666666667 0.102564103 0.0575 0.02564102564102564], ... 
‘String’, Azimuth’, 'Style’,'text’, 'Tag’,’StaticText4’); 


% The Elevation View Editable Text Box 

eleditHndl = uicontrol(Units' 'normalized", 'BackgroundColor',[1 1 1], ... 
'Position',[0.4191666666666666 0.06227 106227106227 0.05 0.03663003 663003663], ... 
Stine 507 Enable. of? styles edi Tas" Edit Texts’) 

eltextHndl = uicontrol('Units' 'normalized", ... 

'BackgroundColor',[0.75294] 0.752941 0.752941], ... 

'Position',[0.4154166666666667 0.102564103 0.0575 0.024420024], ... 

Sine Elevation. style text Tas Static lext5); 


% The Color-Axis Scaling slider & its associated text boxes 
cscalslidHndl 7 uicontrol('Units','normalized', ... 

'BackgroundColor' [0.75294] 0.75294] 0.752941], ... 
'Position',[0.555 0.06227106227106227 0.1475 0.036630036630036063], ... 
'Style';'slider, 'Enable','off, 'Tag','Slider2) 

ecscaltextHndll = uicontrol('Units''normalized', ... 

'BackeroundColor',[0.752941 0.752941 0.752941], ... 
Position',[0.5725 0.102564103 0.1125 0.024420024], ... 
'String','Color-Axis Scaling’, 'Stvle';'text', 'Tag’,’StaticText6’); 

cscaltextHndl2 7 uicontrol('Units','normalized', ... 

'BackgroundColor',[0.75294] 0.752941 0.752941], ... 
'Position',[0.57125 0.0347985347985348 0.0175 0.02564102564102564], ... 
"Sinis 0. sivic text, Tae Sram Text? 

cscaltextHndl3 - uicontrol( Units','normalized', ... 
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'BackgroundColor',[0.75294] 0.752941 0.752941], ... 
'Position',[0.67 0.0347985347985348 0.0175 0.02564102564102564], ... 
"String", 1. Style text. Tas Static Texts 


% The Plot Options menu & sub-menus 

plotHndl = uimenu('Label','Plot‘, ... 
Tag’, PlotMenu'); 

plotimpuls = uimenu(’Parent',plotHnd|, ... 
'Label','Frequencv-Domain Fields', "Tag','Plotuimenul'); 

plotEyl = uimenu('Parent',plotimpuls, 'Callback','nfldvis("plotEv 1"), ... 
'Label','Ev(f)', "Tag',Plotuimenu2); 

plotEzl 7 ummenu( Parent',plotimpuls, 'CallbacK','nfldvis("plotEz1")', ... 
ltabel -EzCi)iealac, Eletuimentions 

plotHx1 = uimenu(Parent',plotimpuls, 'Callback','nfldvis("plotHx1"), ... 
"Label" 'Hx(f), "Tag", Plotuimenu4”); 

plotVin = uimenu('Parent',plotHndl, 'Separator''on', ... 
'Labe!l','V in', 'Tag','Plotuimenu5’); 

plotVint 7 uimmenu( Parent',plotVin, 'Callback','nfldvis("plotV int"), ... 

'Label''Vin(t), 'Tag", Plotuimenu6”); 
plotVinf = uimenu( Parent',plotVin, 'Callback' 'nfldvis("plotVinf”), ... 
'Label', Vin(f)|', "Tag", Plotuimenu”7”); 

plotfiltrd = uimenu('Parent',plotHndl, 'Separator','on', ... 
'Label',""Filtered" Freq Domain Fields’, "Tag';,Plotummenu$); 

plotEy2 7 ummenu( Parent'plotfiltrd, 'Callback','nfldvis("plotEv2" y, ... 
'Label','Ev(f)', 'Tag','Plotuimenu95); 

plotEz2 ^ uimenu( Parent ;plotfiltrd, 'Callback','nfldvis("plotEz2");, ... 
'Label','Ez(f)', 'Tag’,’Plotuimenu 1] 0’); 

plotHx2 = uimenu('Parent',plotfiltrd, 'Callback''nfldvis("plotHx2"), ... 
'Label''Hx(f), 'Tag'Plotuimenul 1”); 


% The Print Options menu  sub-menus 
printHnd! = uimenu(‘Label’,’Print’, "Enable','on', ... 
"Iag','Print Menu); 

prntscreen 7? uimmenu( Parent. printHndl, 'Callback', 'nfldvis("save"), ... 
'Label','Save Figure as .TIF', "Tag','Printuimenul'); 

prntfig = uimenu('Parent',printHndl, 'Callback', 'nfldvis("printaxis1"), ... 
"Label" 'Print Figure Only”, 'Tag','Printuimenu2”); 

prntfigncb ^ uimenu( Parent';printHndl, 'Callback', 'nfldvis("printallaxes"),, ... 
'Label','Print Figure and ColorBar', 'Tag’,’Printuimenu3’); 


% The Colorbar Axes 
cbaraxes=axes('Color',[0.6 0.6 0.6], "Units','normalızed', ... 
'Position'.[0.795 0.3003663003663004 0.04 0.4926739926739927], ... 
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'Visible','off, 'Tag';,'Axes2); 
axes(axHndl) % Reset the current axes to axHndl for plotting/renderine 


% Uncover the figure 

hndlList=[infoHndl startHndl demoHndl movieHndl continuHndl resetHndl closeHndl, ... 
edittextHndl grppopupHndl frampopupHndl fpspopupHndl zscalslidHndl, ... 
zscaltextHndll zscaltextHndl2 zscaltextHndl3 azıedıtHndl azıtextHnd], ... 
eleditHndl eltextHndl cscalslıdHndl cscaltextHndl1 cscaltextHndl2, ... 
cscaltextHndl3 plotHndl plotEy! plotEz1 plotHx1 plotVint plotVinf, ... 
plotEy2 plotEz2 plotHx2 printHndl cbaraxes prntscreen |; 

set(figNumber,'Visible''on", UserData',hndlList) 

set(continuHndl,' Value',1) 

watchoff(oldFigNumber); 

figure(figNumber); 


YY V/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/0/00 00/0/0000 
9/00/09/09/69/09/09/09/09/09/09/99/99/09/69/09/09/09/09/09/0/o9/09/09/60/09/09/09/09/09/09/09/69/69/09/09/09/09/69/09/09/09/o9/o 

case 'info' 
0/59/69/69/09/99/69/09/69/69/69/09/09/69/09/09/09/09/09/09/69/09/09/09/09/09/09/69/09/09/09/69/99/09/09/69/09/09/99/9/69/09/69/09/9/69/0 
YSL VV %%V%%0%%/%0% 


helpwin(m filename) 


0/99/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/o 


W%%%%%%%0%0%%0V0%0%%0%% 
01% VV%%% 0% %0%0%0%0%%%%%%%%%%% Yo 


MMMM %%%%% 
case ‘start’ | 
0/59/99/69/69/09/09/69/69/o9/o9/99/09/09/09/09/09/o9/09/09/09/09/o9/09/09/09/09/09/09/09/09/09/99/09/09/09/09/09/09/09/99/09/09/09/o9/09/o 
0/59/6 9/69/69 9/99/99/99/69/69/69/09/69/69/69/69/99/09/99/69/99/69/99/09/69/99/69/69/09/o 

axHndl = sca; 

figNumber = gcf; 

hndlList = get(figNumber, UserData'); 

infoHndl=hndlList(1); startHndl-hndlList(2); demoHndl-hndlList(3); 

movieHndlI-hndlList(4); continuHndl-hndlList(5); resetHndl-hndlList(6); 

closeHndlI-hndlList(7); edittextHndl-hndlList(8); grppopupHndl-hndlList(9); 

plotHndl-hndlList(24); printHndl-hndlList(33); cbaraxes-hndlList(34); 


0 
/ 
0 
/ 


cla 

set(startHndl,'Enable','off') 

set(demoHndl,'Enable','off') 

set(movieHndl,'Enable','off') 

set(resetHndl,'Enable','on') 

set(plotHndl,'Enable','off') 

text(0,0,'Select the .zeo file containing the AWAS antenna geometry.., ... 
'HorizontalAlignment’, center’); 

pause(le-20) 

set(edittextHndl,'String','Select the .geo file containing the AWAS antenna geometry...) 

[geofile,filepath] = uigetfile(’*.geo', AWAS Geometry Files’); 
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% “*** Reading AWAS Antenna Geometry Data *** 


% WARNING: Ensure the following parameters are met in order to read the correct data 


% from the .GEO file. The # of nodes below assumes that the angle between 
% the v-axis and a line drawn between Nodes | & 5 accurately duplicates the 
% angular dimension of the TEM horn. If this is not so, then the following 
% section will require changes (as will the "filtering" code located 

% immediately prior to the movie function). 


geofid = fopen(geofile,'rt'); 

geohdr = fscanf(geofid,'%76c',1); % Skips the header data on the .GEO file 
gl = fscanf(geofid,'%3d',1); % Skips the first integer value 

GndPlane 7 fscanf(geofid,?o5d',1); 9o Indicates presence of a Ground Plane 
g2 = fscanf(geofid,'%3d',3); % Skips the next 3 integer values 

nodel = fscanf(geofid,'%14e',3); % X. Y. & Z coordinates for Node #1 
node2 - fscanf(geofid, 9o14e',3); 9$ X. Y. & Z coordinates for Node #2 
node3 = fscanf(geofid,'%14e',3); 96 X. Y. & Z coordinates for Node #3 


h *** Readinp AWAS Near-Ficld Data *** 


cla; 
text(0,0,'Select the .nfl file containing the AWAS Near-Field data., ... 
'HorizontalA lignment','center'); 
pause( le-20) 
set(edittextHndl,'String','Select the .nfl file containing the AWAS Near-Field data...’) 
[filename,filepath] 7 uigetfile( *.nfl',' AWAS Near-Field Files"); 
fid = fopen(filename, rt’); 


Ne = fscanf(fid,963d',1); % Number of Position Groups 
Fp = fscanf(fid,'% 1 2e',2); Sol 12) Sta Stop Tregikits 
Ns = fscanf(fıd,'%4d',1); % Number of AWAS Freq Steps 


fl =Fp(1/le6; f2=Fp(2)/le6; df-(f2-flYNs; f= 0:df:f2; 
cla;  set(edittextHndl,'String',") 
text(0,0.1,[ Number of Groups: '‚int2str(Ng)], HorizontalAlignment','center') 
text(0,0,['Number of AWAS frequency steps: ',int2str(Ns)], ... 
'HorizontalAlignment','center') 
text(0,-0.1,['f. 1: ,num2str(f1), MHz: df-: ',num2str(df), ... 
"MHz; f 2:', num2str(f2), MHZ ],'HorizontalAlignment.'center') 
axis([-1 1-1 1]) 
AntennaType = fscanf(fid,'%3d',1) % The Antenna Type (Transmitter/Revr) 


Nd = fscanf(fid,'%3d',2); % Skips the next 5 integer values. 
if AntennaType--2, % For Revr Antennas Only! 


E theta real ^ fscanf(fid,9612e',1); 95 Real Incident E-Field fm Theta Direction 
E theta_imag = fscanf(fid,'%12e',1); % Imag Incident E-Field fm Theta Direction 
E_phi_real = fscanf(fid,'%12e',1); % Real Incident E-Field fm Phi Direction 

E phi imag = fscanf(fid,'"%12e',1); 96 Imag Incident E-Field fm Phi Direction 
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theta = fscanf(fid,'%8f',2); % (Thetal. Theta?) Start/Stop Theta Angeles 
Ntheta = fscanf(fid,'%4d',1); % Number of Steps from Thetal to Theta? 


phi = fscanf(fıd,'%8f',2); % (Phil. Phi2) Start/Stop Phi Angles 
Nphi 7 fscanf(fid,'?64d', 1); % Number of Steps from Phil to Phi? 
end 


xp = zeros(Ng,2); yp=zeros(Ng,2); zp = zeros(Ng,2); 


for ng=1:Ng; % Reading Group Coordinates 

xp(ng,:) = fscanf(fid,'%l 2e',[1,2]); 9^ (xl. X2) Near Fld Pts 

Nxs(ng) = fscanf(fid,'%4d',1); % No.x-steps 

yp(ng,:) = fscanf(fid,'%ol 2e',[1,2]); % (v1. v2) Near Fld Pts 

Nys(ng) = fscanf(fıd, %4d',1); 9o NO. v-steps 

zp(ng,:) = fscanf(fid,?o12e',[1,2]); 96 (zl. z2) Near Fld Pts 

Nzs(ng) = fscanf(fid, ?o4d',1); % No. z-steps 

erpstr(ng,:) = sprintf('%13.0f.ng); % Create a string array of nfld grps 
end 


popupstrl = char('Near-Field Group',grpstr); % Creates the final near-ficld group list 

set(grppopupHndl,'Enable','on''Strine',popupstrl) 

set(continuHndl,'Enable','on') 

set(edittextHndl,'String',' Select the Near-Field group to analyze from the "Near-Field Group" 
pop-up menu below. Press the "Continue" button once you have made your selection...') 

waitfor(continuHndl, Value") 

set(edittextHndl,'String','Calculating the Frequency-Domain E- and H-Fields...") 

watchon; 

set(continuHndl,'Value',1.'Enable' 'off') 

pause(1e-20) 


NG = get(grppopupHndl, Value’); 
set(grppopupHndl,'Enable’, off’) 

if NG>1, NG=NG-1; end 

xl = xp(NG,1); x2 = xp(NG,2); 

yl = yp(NG, 1): y2 = yp(NG,2), 

zl = zp(NG,!); z2 = zp(NG,2):; 

Nx = Nxs + ones(size(Nxs)); Ny = Nys + ones(size(Nys)); 
Nz = Nzs + ones(size(Nzs)); % Total # of x.y & z points 
dx = (x2-x1)/Nxs(NG); x = x1:dx:x2; 

dy = (y2-y1)/Nys(NG); y = yl:dy:y2; 

dz = (z2-z1)/Nzs(NG); z 7» zl:dz:z2; 


if z1==z2, % AWAS near-field data is ın the xy -plane 
N1 =Nx(NG) Nls=Nxs(NG); N2=Ny(NG); N2s= Nys(NG); 
rowposit = 1:N2; colposit = (N1+1)/2; xdim = x; xlabl 7 'x'; ydim 7 y; 
ylabl = 'v'; minxdim 7 x1; maxxdim 7 x2: minydim 7 yl; maxydim 7 y2; 
else % AWAS near-field data is in the vz-plane 
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NI=Ny(NG); Nls=Nys(NG); N22 Nz(NG); N2s » Nzs(NG); 
rowposit = (N2+1)/2; colposit = 1:N1; xdim = y; xlabl = 'y'; ydim = z; 
ylabl = 'z'; minxdim = yl; maxxdim = y2; minydim = zl; maxydim = z2; 
end 
Nf = Ns+2; % Number of points to f2 folding freq (needs to be 2^m-« 1) 
Nrow = 2*Nx.*Ny.*Nz; EH = zeros(6,Nrow(NG)); 
Eyf = zeros(N2,N1,Nf); ifGndPlane==1 & z1~=z2, Eyf = zeros(2*N2s+1,N1,Nf); end 
Ezf -Eyf; Hxf-Ezf, 





for m=2:Nf; % Frequency Stepping 
for n=1:Ng; % Reading Group Data 
S — fgetl(fid); % Read blank line between EH data arrays 
EH = fscanf(fid,'% 12e',[6, Nrow(n)]); 
He = ING % Saving Data at Selected Group 


Eyf(1:N2,:,m) = reshape(EH(3,1:2:Nrow(NG)-1)'+j*EH(4,1:2:Nrow(NG)-1)',N2,N]); 
Ezf(1:N2,:,m) = reshape(EH(5, 1:2: Nrow(NG)-1)'+j*EH(6, 1:2: Nrow(NG)-1)',N2,N1); 
Hxf(1:N2,:,m) = reshape(EH(1,2:2: Nrow(NG))'t+j* EH(2,2:2:Nrow(NG))',N2,N 1); 


end 
end 
end 
if GndPlane== 
Eyf=.5*Eyf; Ezf=.5*Ezf; Hxf=.5*Hxf; % Scale factor to comp for endplane use 
if z1~=z2, % Changes the appropriate variables duc to 


N2s-2*N2s; N2=N2st+1; zl—z2; % the presence of the ground plane. 
dz=(z2-z1)/N2s; z-zl:dz:z2; rowposit-(N241)2; ydim-z; minydim-zl; 
Eyf((N2+1)/2:N2,:,:) = Eyf(1:N2+1)/2,:,:); % Reshape the Evf matrix to 
Eyf(1:N2s/2,:,:) = flipdim(Eyf((N2+3)/2:N2,:,:),1); % fit the end plane "image" 
EZIGQIN2-3:DY2:N2. y EzICI(N2ZSS DASS % Reshape the Ezf matrix to 
Ezf(1:N2s/2,:,:) = flipdim(Ezf((N2+3)/2:N2,:,:),1); % fit the end planc "image" 
Axf((N2+1)/2:N2,:,:) = Hxf(1:(N2+1)/2,:,:); % Reshape the Hxf matrix to 
Hxf(1:N2s/2,:,:) = flipdim(Hxf((N2+3)/2:N2,:,:),1); % fit the gnd plane "Image" 
end 

end 

fclose(fid); fclose(geofid); clear EH 


save Eyf Eyf 

save Ezf Ezf 

save Hxf Hxf 
watchoff 


% Retrieve the uscr's ID.info for all plots: 
set(edittextHndl,’String’,’ Type brief 1D information for all plots in the input box below...") 
IDinfo 7 inputdlg( ('Enter brief ID information for all plots:’}, ... 
ID Info for All Plots',1); 
IDinfo 7 char(IDinfo): 
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watchon; 


?7/o Create the Gaussian Impulse Input Voltage. Vin(t). & |Vin(f)i. the Fourier 
% Transform of Vin(t): 
set(edittextHndl,'String’,'Creating the Gaussian impulse input voltage...') 
pause(le-20) 
Nt = 2*Nf-2; % Number of radix-2 frequency and time points (22m) 
Tr= 1/df; dt=Tr/Nt; T=Tr-dt: t= 0:dt:T; 
Vint = zeros(1,Nt);  Vinf = Vint; 
tl = 0:29; 
if AntennaType == 
sigmal 21; sigma2=2; Al =sıgma2/(sigma2-sıgmal); A2=1-Al; 
vin 7 Al*exp(-(t1-14).^2./(2*sigmal^2)) ^ A2*exp(-(t1-14).^2./(2*sigma2^2)); 
else 
vin 7 exp(-(t1-14).^2./2); 
end 
Vint(1:30) 7 vin; 
fp= f(l:Nf);  Vinf- le-O*fft(Vint); Vif= abs(Vinf(1:Nf); 


% Frequency Scaling of AWAS Computed Fields Due to Gaussian Input Voltage 
set(edittextHndl,'String','Scaling the E- and H-Fields due to the Gaussian Input Voltage...") 
pause(1e-20) 


Vif2 7 zeros(1, 1, Nt); % Transform Vif from a vector to a 
Vif2(1,1,:) = Vinf; % 3-D array of N2 rows. ! column. & 
Vin 7 repmat(Vif2,[N2,1,1]); 9o Nt "tiles" (the 3rd dimension) 


Eyf=-Eyf,  Ezf--Ezf, Hxf--Hxf; % Corrects for a 180 degree phase shift 


if AntennaType == 
k = 2* pi/3e2* fp; % Create a time delay for the case 
d = zeros(1,1,Nf); % of the antenna acting as a receiver. 
yo = 3; % Initial location of incident wave. 


d(1,1,:) = exp(-j* yo*k); 
delay = repmat(d,[N2,1,1]); 


eta = 120* pi; % Intrinsic Impedance (~377 ohms). 

Einc = zeros(N2, 1,length(fp)); % Create the incident fields for the case 
Einc(:,1,:) = exp(j*y'*k); % of the antenna acting as a receiver. 
form >= 1:N] % Calculate the total fields. 


Eyf(:,m,:) = Vin(:; 1, 1: Nf).*delay.* Eyf(:,m,:); 
EX uam = Vin LIND, "delay Ezi m. we 
Hxf(:,m,:) » Vin(:; 1, 1T: Nf).*delay.* Hxf(:,m,:); 

end 

else 

for m =1:N] 
Eyf(:,m,:) = Vin(:,1,1:N5.*Eyf(:,m,:); 
Ezf(:,m,:) = Vin(:,1,1:Nf).*Ezf(:,m,:); 
Hxf(:,m,:) = Vin(:,1,1:Nf).*Hxf(:,m.:); 
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end 
end 


save Hdr f y rowposit colposit Nf IDinfo T t Vint fp Vif 


save Eyff Eyf % Saving variables for all "Plot" functions 
save Ezff Ezf 
save Hxff Hxf 


% Convert the Frequency-Domain Data to a Time-Domain Representation 
set(edittextHndl,’String’,’Calculating the Time-Domain Fields...’) 
pause( 1e-20) 
Al = zeros(N2,N1,Nt); A2=Al; A3=Al; 
Eyt- Al; Ezt=Al; Hxt=Al; fscl=2e6*f; 
Al(:,:,1:Nf) = Eyf;  A2(:,:,1:Nf) = Ezf;  A3(:,:,1:Nf) = Hxf; 
clear Einc Eyf Ezf Hxf Vif Vint Vinf Vink ddelay % Free up some memory 
Al(:,:,NE:Nt) = Al(,:,Nf:Nt) + conj(flipdim(A 1(:,:,2:Nf),3)); 
A2(:,:,Nf:Nt) = A2C:,:,Nf:Nt) + conj(flipdim(A2(:,:,2:Nf),3)); 
A3(:,:,Nf:Nt) = A3(:,:,Nf:Nt) + conj(flipdim(A3(:,:,2:Nf),3 )); 


for m=1:N2 
for n=1:N 1 
Eyt(m,n,:) 7 fscl*real(ifft(A 1(m,n,:))); 
Ezt(m,n,:) 7 fscl*real(ifft(A2(m,n,:))); 
Hxt(m,n,:) 7 fscl*real(ifft( A3(m,n,:))); 
end 
end 
clear Al A2 A3 % Free up more memory 


% "Filter" the Eyt. Ezt. & Hxt Data for the YZ-Plane Casey This 
% 1s done to remove the high values of the fields directly on top 
% of the wire and those outside the wire structure: 
if x] == x2, 
a = (node3(3)-nodel(3))/(node3(2)-nodel(2));  9?$ a-dz/dw. slope of "V" 
b=-a*yl; % Offset of line to ensure passage thru lsty point of ncar-field 
Z] =a*y + b*ones(sıze(y)); 


Z2 --Zl; 
M = zeros(size(Eyt)); M1 = zeros(N2,N1); 
for N = 1:Ny; 

M1(:,N) = flipIr(z)'>Z1(N) | fliplr(z)'<Z2(N); 
end 


M 7 repmat(M1,[1 1 Nt]): 

Eyt = Eyt.*~M; Ezt- Ezt.*-M; Hxt = Hxt.*-M; 
end; clear M MI M2 
watchoff 


save Movihdr T Nt dt minxdim maxxdim minvdim maxvdim xdim ydim xlabl ylabl 


save Eyt Eyt 7o Saving variables for use with "Movie" function 
save Ezt Ezt 
save Hxt Hxt 


clear Eyt Ezt Hxt % Free up more memory 

set(edittextHndl,'String',") 

set(grppopupHndl,' Value',1l) 

set(figNumber,'UserData',hndlList); 

cla; 

text(0,0.1,'All calculations are complete.' 'HorizontalA lignment''center') 

text(0,0,'Press the "Movie" button to create a movie of the time-domain fields.', ... 
'HorizontalAlignment','center') 

text(0,-0.1,'or plot the frequencv-domain representation of the fields using the "Plot" pull-down 

menu at the top of this window.”, ... 

'HorizontalAlignment','center') 

axis([-] ] -1 1]) 

set(movieHndl,'Enable','on') 

set(plotHndl,'Enable','on') 


0/59/99/99/59/69/99/99/69/99/99/99/99/99/69/69/99/99/69/69/9 9/59/69/59/59/, 9/5 9/59/5 9/595 9/5 9/6 9/69/59/69/,9/5 9/59/6976 %% W% %% 
0/99/69/99/59/69/69/69/69/09/99/o09/09/99/99/09/09/09/09/09/09/09/99/99/09/09/09/o9/o9/99/99/69/69/09/69/o 9/09 99/9 9/09/69/99/69/o 
case 'demo 
0/99/99/99/99/09/99/09/99/909/09/09/09/09/09/09/09/09/99/99/09/09/09/09/09/09/09/09/o9/09/o 9/o9/09/09/09/09/69/09/09/99/0 9/09/09/09/o 9/0 9/0 
0/99/09/09/99/09/09/09/o9/99/09/09/99/09/09/09/09/09/09/09/o9/09/09/09/69/09/99/09/09/o9/o 
axHndl = gca; 
figNumber = gcf; 
hndlList = get(figNumber,'UserData'); 
infoHndl-hndlList(1); startHndl-hndlList(2); demoHndl-hndlList(3); 
movieHndl-hndlList(4); continuHndl-hndlList(5); resetHndl-hndlList(6); 
closeHndl-hndlList(7); edittextHndl-hndlList(8); grppopupHndl-hndlList(9); 
frampopupHndl-hndlList(10); fpspopupHndl-hndlList(11); zscalslidHndl-hndlList(12); 
azieditHndl-hndlList(16); eleditHndl-hndlList(18); cscalslidHndl-hndlList(20); 
plotHndl-hndlList(24); printHndl-hndlList(33); cbaraxes-hndlList(34); 


cla 

set(startHndl.'Enable','off") 

set(demoHndl,'Enable','off) 

set(movieHndl,'Enable','off") 

set(resetHndl,'Enable','on') 

set(grppopupHndl,'Enable','on') 

set(plotHndl,Enable','off') 

text(0,0,'Select the .geo file containing the AWAS antenna geometry.', ... 
'HorizontalA lignment' center”); 

pause(1e-20) 

set(edittextHndl,'String','Select the .geo file containing the AWAS antenna geometry...) 

[geofile,filepath] = uigetfile('demo.geo';, AWAS Geometry Files’); 


if -stremp(geofile, Demo.geo') 
er — errordlg( Y ou may onlv select the "Demo.geo" file for this demo!‘ ERROR’); 
waitfor(er, ButtonDownFen') 

end 

geofile ='Demo.geo'; 


% *** Reading AWAS Antenna Gcometry Data *** 


% WARNING: Ensure the following parameters are met in order to read the correct data 


% from the .GEO file. The # of nodes below assumes that the angle between 
% the y-axis and a line drawn between Nodes | & 3 accurately duplicates the 
% angular dimension of the TEM horn. If this is not so. then the following 
% section will require changes (as will the "filtering" code located 

% immediately prior to the movie function). 


geofid = fopen(geofile, rt’); 

geohdr = fscanf(geofid,'%76c',1); % Skips the header data on the .GEO file 
gl = fscanf(geofid,'%3d', 1); % Skips the first integer value 

GndPlane = fscanf(geofid,'%3d',1); % Indicates presence of a Ground Plane 
22 = fscanf(geofid, ?03d',3); % Skips the next 3 integer valucs 

node! = fscanf(geofid,'%14e',3); 96 X. Y. & Z coordinates for Node #1 
node2 - fscanf(geofid,9o14e',3); 96 X. Y. & Z coordinates for Node #2 
node3 7 fscanf(geofid,?9614e',3); 95 X. Y. & Z coordinates for Node #3 


% *** Simulated Reading of AWAS Near-Ficld Data *** 


cla; 

text(0,0,'Select the .nfl file containing the AWAS Near-Field data.', ... 
'HorizontalAlıgnment','center'); 

pause(le-20) 

set(edittextHndl,'String','Select the .nfl file containing the AWAS Near-Field data...) 

(filename, filepath] = uigetfile('demo.nfl',' AWAS Near-Field Files”); 

cla; 

text(0,0.1,"' Number of Groups: l','HorizontalAlignment','center') 

text(0,0,'Number of AWAS frequency steps: 127", ... 
'HorizontalAlignment','center') 

text(0,-0.1,f_l:4 MHz; df=: 4 MHz; f_2: 512 MHz, ... 
'HorizontalAlıgnment','center') 

axis([-1 ] -1 1]) 


grpstr 7 sprintf('9o13.0f,1); % Create a string array of nfld erps 

popupstrl = char('Near-Field Group',grpstr); 

set(grppopupHndl,' String',popupstrl) ^6 Creates the final near-field group list 

set(continuHndl,'Enable','on') 

set(edittextHndl,'String','Select the Near-Field group to analyze from the "Near-Field Group' 
pop-up menu below. Press the "Continue" button once you have made your selection...") 

waitfor(continuHndl,' Value") 
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set(edittextHndl,'String','Calculating the Frequency-Domain E- and H-Fields...') 
watchon; 

set(continuHndl,' Value',1,'Enable','off) 

set(grppopupHndl,'Enable','off') 

pause(5) 

fclose(geofid); 

load DemoHdr. load DemoEzf 

watchoff 


set(edittextHndl,'String'," Type brief ID information for all plots in the input box below...) 
IDinfo = inputdlg({'Enter brief ID information for all plots:’}, ... 
'1D Info for All Plots', 1); 
IDinfo = char(IDinfo); 
pause(1e-20) 


9o Plot a sample ot the "Impulse Response" Fields in the frequency domain: 

cla; set(axHndl,'V isible','on') 

mesh(f,y,squeeze(abs(Ezf(rowposit,colposit, 1 :Nf)))) 

title(['Near-Field E z(v.f)| for x-0 m; z-0 m; ';IDinfo]) 

ylabel('y (meters), xlabel('Freq (MHz)), zlabel(E z| (V/(m-Hz))) 

axis([0 max(f) 0 max(y) 0 1.1 *max(max(squeeze(abs(Ezf(rowposit,colposit, | :Nf)))))]) 

rotate3d on 

colorbar(cbaraxes); 

Data=1:64;Data=(Data'* Data)/64: 

h-msgbox([' Place the mousepointer on the plot, then press ';'and hold the left mouse button 
to rotate the plot.'],"To Rotate Plot:','custom',Data,hot(64)); 

set(continuHndl,'Enable','on') 

set(edittextHndl,' String',char('This is a plot of the E-Field in the frequency-domain.','Press the 
"Continue" button once you are through with this plot...')) 

pause(1e-20) 

waitfor(continuHndl,' Value") 

set(continuHndl,Value',1 'Enable' 'off) 

rotate3d off 

set(cbaraxes, Visible' 'off) 

set(axHndl,'Visible','off') 

axes(axHndl) 


% Plot the Gaussian Impulse Input Voltage, Vin(t): 
cla reset; ,plot(t, Vint), title([' Antenna Input Voltage V_(iny(t): ',IDinfo]); 
xlabel('Time (uSec)^), ylabel(' Volts), axis([0 T -0.1 max(Vint)+0.1]) 
set(continuHndl, Enable','on') 
set(edittextHndl,'String';char( This is a plot of the Gaussian impulse input voltage, 
Vin(t).',Press the "Continue" button once you are through with this plot...) 
waitfor(continuHndl,' Value") 
set(continuHndl,' Value',1,'Enable','off) 


105 


?o Plot |Vin(f). the Fourier Transform of Vin(t): 
cla; plot(fp, Vif), title(['Antenna Input Voltage |V_{in}(f)]: ',[Dinfo]) 
xlabel('Freq (MHz)), axis([0 f(Nf) 0 max(Vif)]) 
set(continuHndl,'Enable','on') 
set(edittextHndl,'String,char('This is a plot of |Vin(f)|, the frequency-domain representation of 
Vin(t).', Press the "Continue" button once you are through with this plot...')) 
waitfor(continuHndl,' Value") 
set(continuHndl,' Value, l,'Enable','off) 
set(edittextHndl,'String','Scaling the E- and H-Fields due to the Gaussian Input Voltage...) 
watchon; 
pause(5) 
watchoff 


% Displaying a sample of the "Filtered" Frequency-Domain Fields vs v: 

load DemoEzff 

cla; set(axHndl,'Visible','on"), set(cbaraxes, Visible','on') 

mesh(f,y,squeeze(abs(Ezf(rowposit,colposit, 1 :Nf)))) 

title(['Near-Field |E z(y,f)| for x20 m; z-0 m; ',IDinfo]) 

ylabel('y (meters) ), xlabel(Freq (MHz)), zlabel(|E z| (V/(m-Hz))) 

axis([0 max(f) 0 max(y) 0 1.1 *max(max(squeeze(abs(Ezf(rowposit,colposit, | :Nf)))))]) 

rotate3d on 

colorbar(cbaraxes); 

h=msgbox([' Place the mousepointer on the plot, then press ';'and hold the left mouse button 
to rotate the plot.'], To Rotate Plot:','custom',Data,hot(64)); 

set(continuHndl,'Enable','on") 

set(edittextHndl,'String'.char( This is a plot of the E-Field due to the Gaussian input 
voltage.','Press the "Continue" button once vou are through with this plot...')) 

pause(1e-20) 

waitfor(continuHndl,' Value"); 

set(continuHndl,'Value',1,'Enable', off”) 

rotate3d off 

set(axHndl,'V isible','off‘) 

axes(axHndl) 

set(edittextHndl,'String','Calculating the Tıme-Domain Fields...') 

watchon; 

pause(5) 

watchoff 


% Generate € Display the Near-Fields Ex(t). Ez(b). € Hx(t) Using 
% the Movie Format: 
Clean Ze 
load DemoMovihdr, load DemoEzt 
Tns = 1000*T; dt=Tns/(Nt-1); t=O0:dt:Tns; % Time tn nanoseconds 
frames = 0:50:Nt; 
for n=1:length(frames)-] 
framestr(n,:) 7 sprintf('9613.0f .frames(n-1)); 
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end 

popupstr3 = char(' Total Frames',framestr); 

set(frampopupHndl,'String',popupstr3) % Creates the Total Frames pop-up list 

cla; 

text(0,0.1,[T = 'Ànum2str(T),' uSec']) 

text(0,0,['N =',int2str(Nt),’ Samples']) 

text(0,-0.1,['dt = '¡num2str(dt),' nSec']) 

axis([-1 1 -1 1]) 

set(frampopupHndl,'Enable','on','String',popupstr3) 

set(continuHndl, Enable','on') 

set(edittextHndl,' String',' Select the number of frames for the movie from the "Total Frames" 
pop-up menu below. Press the "Continue" button once you have made your selection...') 

waitfor(continuHndl,' Value") 

set(continuHndl, Value, 1,'Enable','off') 

Mt = 50*(get(frampopupHndl,' Value')-1); 

if Mt==0, Mt=25; end 

set(fpspopupHndl,'Enable','on') 

set(continuHndl,'Enable','on') 

set(edittextHndl,'String','Select the frame rate of the movie from the "Frames-per-sec" pop-up 
menu below. Press the "Continue" button once you have made vour selection...') 

waitfor(continuHndl,' Value") 

set(continuHndl,' Value',1,'Enable','off') 

fps = get(fpspopupHndl,' Value’); 

if fps>=3, fps=(fps-1)*2; end 

set(frampopupHndl,'Enable','off) 

set(fpspopupHndl,'Enable','off‘) 

set(edittextHndl.'String',' Now, create a movie to illustrate the E- or H-Field in the time domain. 
Use this frame from your movie to create the axes, coloraxis, and view you desire...) 

pause(5) 

cla; set(axHndl,'Visible','on') 

rotate3d on % Creates a moveable mesh plot for scaling axes prior to movie creation 

Fieldmax = max(max(max(Ezt))); Fieldmin = min(min(min(Ezt))); 

mesh(xdim,ydım,Ezt(:,:,25)) 

vmax = max(abs(Fieldmax),abs(Fieldmin)); 

vaxis = [minxdim maxxdim minydim maxydim -vmax vmax]; 

axis(vaxis), cmatrix = [flipud(jet(64)); jet(64)]; 

colormap(cmatrix), caxis([-vmax vmax]) 

set(cbaraxes, Y Lim',[-vmax vmax]), colorbar(cbaraxes), axes(axHndl) 

set(zscalslidHndl,'Enable','on',' Value',1) 

set(azieditHndl,'Enable','on') 

set(eleditHndl,'Enable','on') 

set(cscalslidHndl,'Enable','on', Value',1) 

while 1 % Sets up the loop to customize display of data 

set(continuHndl,'Enable','on') 
set(edittextHndl,'String';'To customize your movie, adjust the z-axis & c-axis sliders at the 

bottom of the window. Press the "Continue" button once you have made vour adjustments...) 
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waitfor(continuHndl,' Value") 
set(continuHndl,'Enable','off'," Value',1) 
set(edittextHndl,'String',") 
zscale = get(zscalslidAndl,' Value'); if zscale==0, zscale=0.01; end 
cscale = get(cscalslidHndl, Value; if cscale==0, cscale=0.01; end 
[azı el] = view; % Reads the adjusted azi & el 
set(azieditHndl,'String',¡num2str(azi)) 
set(eleditHndl,'String',num2str(el)) 
vaxis = [minxdim maxxdim minydim maxydim -zscale* vmax zscale*vmax]; 
mesh(xdım,ydım,Ezt(:,:,25)) 
view([azi el]), axis(vaxis) 
cmatrix = [flipud(jet(64)); jet(64)]; colormap(cmatrix) 
axHndl=gca; caxis([-cscale*vmax cscale*vmax]) 
set(cbaraxes,'Y Lim',[-cscale*vmax cscale*vmax]), colorbar(cbaraxes); axes(axHndl) 
mulligan 7 questdlg( Would vou like to make further adjustments to the plot?', 'Continue 
Adjustments?’,'No’); 
if stremp(mulligan,'No’), break, end 
end 
rotate3d off 
set(zscalslidHndl,’Enable’,'off’) 
set(azieditHndl,' Enable’, off) 
set(eleditHndl,'Enable','off") 
set(cscalslidHndl,'Enable';'off") 


ela: 
FieldMovi = moviein(Mt); % Create the movie 
set(gca,'NextPlot','replacechildren') 
for nt = 1:Mt 
cla; 
mesh(xdim,ydım,Ezt(:,:,nt)) 
view([azi el]), axis(vaxis) 
caxis([-cscale*vmax cscale*vmax]) 
xlabel([xlabl,' (meters)' ]), ylabel([ylabl,' (meters) ]) 
zlabel('Ez(t) (V/m)), title(['Frame Time: ',num2str(t(nt)),' nSec']) 
FieldMovi(:,nt) = getframe; 
end 
title(" ^), 
movie(FieldMovi,1,fps); 


while 1 
mulligan = questdlg(' Would you like to view the movie again”, ... 
‘Instant Replay’,'No’); 
if stremp(mulligan,'No'), break, end 
movie(FieldMovi,1,fps); 
end 
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set(grppopupHndl,' Value',l) 

set(fpspopupHndl, Value',1) 

set(frampopupHndl,' Value',1) 

set(figNumber,' UserData',hndlList); 

cla; set(axHndl,'Visible','off") 

text(0,0," That concludes the demo...!, ... 
‘'HorizontalAlignment’,'center’) 

axis([-1 1 -1 1]) 

pause(2) 

clear Ezt FieldMovi 

eval('nfldvis(" reset"); 


0/99/09/09/09/09/09/09/09/09/09/09/o9/09/60/00/09/09/09/09/09/69/09/69/09/o9/09/09/09/09/o9/o9/o9/o0/09/09/09/09/09/00/09/09/99/o9/69/o9/o 
0/00/9/o9/99/o9/09/09/09/09/o9/09/09/09/09/99/909/00/69/09/09/09/09/09/09/69/09/09/09/09/09/09/09/00/09/09/09/09/00/09/09/99/09/o9/o 
case 'movie' 
0/99/69/99/99/99/69/99/99/69/09/69/69/09/99/69/99/99/99/99/09/69/99/09/99/99/69/09/09/09/09/09/09/09/09/09/09/09/09/09/99/99/69/69/99/99/o 
0/09/06 9/6 9/69/69/09/09/69/09/09/69/09/09/09/09/09/09/09/09/09/09/09/o9/09/09/09/09/69/09/0 

ela: 

axHndl = gca; 

figNumber = gcf; 

hndlList — get(figNumber, UserData ); 

infoHndl-hndlList(1); startHndl-hndlList(2); demoHndl-hndlList(3); 

movieHndl-hndlList(4); continuHndl-hndlList(5); resetHndl=hndlList(6); 

closeHndl-hndlList(7); edittextHndl-hndlList(8); grppopupHndl-hndlList(9): 

frampopupHndl-hndlList(10); fpspopupHndl-hndlList(11); zscalslidHndl-hndlList( 12); 

azieditHndl-hndlList(16); eleditHndl-hndlList(18); cscalslidHndl-hndlList(20); 

plotHndl-hndlList(24); printHndl-hndlList(33); cbaraxes-hndlList(34); 

set(startHndl,'Enable','off') 

set(demoHndl,'Enable','off') 

set(movieHndl,'Enable','off) 

set(plotHndl,'Enable','off) 

movieMode = questdlg('Do you want to create a new movie or replay a saved one?', ... 

‘What kind of movie?" Create Movie", Replay','Replay'); 


if stremp(movieMode, Replay’) 
cla: 
text(0,0,' Select the .mat file containing the recorded movie.,, ... 
'HorizontalAlignment','center') 
pause(le-20) 
set(edittextHndl,'String','Select the .mat file containing the recorded movie...") 
[filename,filepath] — uigetfile( *.mat',' Movie Files’); 
cla; set(edittextHndl,'String',[ Loading the movie file, ';filename]), pause(1e-20) 
eval(['load ',filepath,filename,' -mat ]); 
set(fpspopupHndl,'Enable','on') 
set(continuHndl,'Enable','on') 
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set(edittextHndl, String','Select the frame rate of the movie from the "Frames-per-sec" pop-up 
menu below. Press the "Continue" button once you have made your selection...) 

waitfor(continuHndl,' Value") 
set(continuHndl,' Value',l,'Enable','off') 
fps = get(fpspopupHndl,' Value"); 
if fps>=3, fps=(fps-1)*2; end 
set(fpspopupHndl,'Enable','off) 
set(azieditHndl, String',num2str(azi)) 
set(eleditHndl,'String',num2str(el)) 
axis(vaxis), axisCon'), view([azi el]) 
xlabel([xlabl,' (meters)']), ylabel([ylabl,' (meters)']), zlabel(zlabl); 
colormap(cmatrix) 
caxis([-cscale*vmax cscale*vmax]), set(cbaraxes,' Y Lim',[-cscale*vmax cscale*vmax]) 
colorbar(cbaraxes), axes(axHndl) 
wnumber = watchon; 
movie(axHndl,FieldMovi, 1 ,fps); 
watchoff(wnumber); 
while | 

mulligan = questdig(' Would you like to view the movie again”, ... 

Instant Replay','No); 

if stremp(mulligan,'No'), break, end 

movie(FieldMovi,1,fps); 
end 
cla; set(axHndl,'Visible','off'‘); 
set(fpspopupHndl,'Value',1); 
set(figNumber,'V isible','on', "UserData',hndlList); 
clear FieldMovi 
eval('nfldvis("reset”)); 


else 

text(0,0,'Select the .mat file containing the header data (default is MOVIHDR.MAT)...', 
'HorizontalAlignment','center') 

axis([-1 ] -1 1]) 

pause(1e-20) 

set(edittextHndl,'String','Select the .mat file containing the header data (default is 
MOVIHDR.MAT)...') 

[filename,filepath] — uigetfile('movihdr.mat',' Movie Header Data File"); 

eval(['load ',filepath,filename,' -mat' |); 

set(edittextHndl,'String',") 

cla; 

filetype = questdlg('Select the E- or H-Field you want to view:','Field 
Type, Eyt. Ezt Hirt Ezi); 

text(0,0,'Select the .mat file containing the time-domain E- or H-Fields (EYT.MAT. 
EZT.MAT, or HXT.MAT)...', 'HorizontalAlignment’, center’); 

pause(1e-20) 
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set(edittextHndl,' String','Select the .mat file containing the time-domain E- or H-Fields 
(EYT.MAT, EZT.MAT, or HXT.MAT)...”) 

[filename, filepath] = uigetfile([filetype '.mat']Ey(t). Ez(t). or Hx(t) Data File”); 

eval(['load ',filepath,filename,' -mat ]): 

set(edittextHndl,'String',") 


% Generate & Display the Near-Fields Ey(t), Ez(t), & Hx(t) Using the Movie Format: 
Tns = 1000*T; dt=Tns/(Nt-1); t=O:dt:Tns; % Time in nanoseconds 
frames = 0:50:Nt; 
for n=1:length(frames)- 1 
framestr(n,:) = sprintt('%13.0f',frames(n+1)); 
end 
popupstr3 = char('Total Frames',framestr); 
set(frampopupHndl,'String',popupstr3) — ?6 Creates the Total Frames pop-up list 


while 1 % Loop for adjusting the parameters of the movie 

cla; 

text(0,0.1,['T 7 ''num2str(T),' uSec']) 

text(0,0,['N 7 'int2str(Nt),' Samples']) 

text(0,-0.1,['dt ^ ,num2str(dt),' nSec']) 

axis([-1 ] -1 1]) 

set(frampopupHndl,'Enable','on','String';popupstr3) 

set(continuHndl,'Enable','on') 

set(edittextHndl,'String','Select the number of frames for the movie from the "Total Frames’ 
pop-up menu below. Press the "Continue" button once you have made your selection...') 

waitfor(continuHndl,' Value") 

set(continuHndl,' Value',1,' Enable','off) 

Mt = 50*(get(frampopupHndl, Value')-1); % Total # of frames to use in Movie 

if Mt==0, Mt=25: end 

set(fpspopupHndl,'Enable','on') 

set(continuHndl,'Enable','on') 

set(edittextHndl,'String','Select the frame rate of the movie from the "Frames-per-sec" pop- 
up menu below. Press the "Continue" button once you have made your selection...) 

waitfor(continuHndl, Value") 

set(continuHndl,' Value',1,'Enable',' off) 

fps = get(fpspopupHndl,' Value"); % Frames-per-sec for the Movie 

if fps>=3, fps=(fps-1)*2; end 

set(frampopupHndl,'Enable'.'off") 

set(fpspopupHndl,'Enable','off') 

set(edittextHndl,'String',"This is a sample mesh plot from your movie with the default 
settings..." 

pause(2) 

cla; set(axHndl, Visible','on') 

rotate3d on % Creates a moveable mesh plot for scaling axes prior to mov ie creation 

ıf stremp(filetype,'Eyt') 

Fieldmax = max(max(max(Eyt))); Fieldmin = min(min(min(Eyt))); 


t 
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mesh(xdim,ydim,Eyt(:,:,37)), zlabl 2'E v(t) (V/my; 
xlabel([xlabl,' (meters)']), ylabel([ylabl,' (meters)]) 
zlabel([zlabl]), title(['Frame Time: '‚num2str(t(37)),' nSec']) 
elseif stremp(filetype,'Ezt') 
Fieldmax = max(max(max(Ezt)));  Fieldmin = min(min(min(Ezt))); 
mesh(xdim,ydim,Ezt(:,:,37)), ^ zlabl-'E z(t) (V/m); 
xlabel([xlabl,' (meters) ]), ylabel([ylabl, (meters)']) 
zlabel([zlabl]), title(['Frame Time: '‚num2str(t(37)),' nSec']) 
else 
Fieldmax ^ max(max(max(Hxt));  Fieldmin = min(min(min(Axt))); 
mesh(xdim,ydim,Hxt(:,:,37)), zlabl 2'H x(t) (A/m); 
xlabel([xlabl,' (meters) D), ylabel([ylabl,' (meters)']) 
zlabel([zlabl]), title(['Frame Time: ';Ànum2str(t(37)),' nSec']) 
end 
azi — str2num(get(azieditHndl,'String")); 
el = str2num(get(eleditHndl,'String')); 
view([aziel]) % If rebuilding the movie, this holds the last view paramcters 
vmax = max(abs(Fieldmax),abs(Fieldmin)); 
vaxis = [minxdim maxxdim minydim maxydim -vmax vmax]; 
axis(vaxis), cmatrix = [flipud(jet(64)); jet(64)]; 
colormap(cmatrix), caxis([-vmax vmax]) 
set(cbaraxes, Y Lim',[-vmax vmax]), colorbar(cbaraxes), axes(axHndl) 
set(zscalslidHndl, Enable','on',' Value',1) 
set(azieditHndl,'Enable','on') 
set(eleditHndl,'Enable','on') 
set(cscalslidHndl,'Enable','on', Value',1) 
while 1 % Sets up the loop to customize display of data 
set(continuHndl,'Enable','on') 
set(edittextHndl,'String',"To customize your movie, adjust the z-axis & c-axis sliders at the 
bottom of the window. Press the "Continue" button once you have made your adjustments...) 
waitfor(continuHndl,'Value") 
set(continuHndl,'Enable','off,'Value',1) 
set(edittextHndl,'String',") 
zscale 7 get(zscalslidHndl,' Value); if zscale==0, zscale=0.01; end 
cscale = get(cscalslidHndl,' Value’); if cscale==0, cscale=0.01; end 
[azi el] = view; % Reads the adjusted azi & el 
set(azieditHndl,'String',num2str(azi)) 
set(eleditHndl,'String',num2str(el)) 
vaxis = [minxdim maxxdim minydim maxydim -zscale*vmax zscale*vmax]; 
ıf stremp(filetype,'Eyt') 
mesh(xdim,ydim,Eyt(:,:,37)) 
xlabel([xlabl,' (meters) ]), vlabel([ylabl,' (meters) ]) 
zlabel([zlabl]), title(['Frame Time: ',num2str(t(37)),' nSec']) 
elseif stremp(filetype,'Ezt') 
mesh(xdim,ydim,Ezt(:,:,37)) 
xlabel([xlabl,' (meters)']), ylabel([ylabl,' (meters)']) 


zlabel([zlabl]), title(['Frame Time: '‚num2str(t(37)),' nSec']) 
else 
mesh(xdim,ydim,Hxt(:,:,37)) 
xlabel([xlabl,' (meters) ]), vlabel([ylabl,' (meters) ]) 
zlabel([zlabl]), title(['Frame Time: 'Ánum2str(t(37)),' nSec']) 
end 
view([azi el]), axis(vaxis) 
cmatrix 7 [flipud(jet(64)); jet(64)]; colormap(cmatrıx) 
axHndl=gca; caxis([-cscale*vmax cscale*vmax]) 
set(cbaraxes, Y Lim',[-cscale*vmax cscale*vmax]), colorbar(cbaraxes); axes(axHndl) 
mulligan 7 questdlg( Would vou like to make further adjustments to the plot?', Continue 
Adjustments”',No'); 
if stremp(mulligan,'No'), break, end 
end 
rotate3d off 
set(zscalslidHndl,'Enable','off‘) 
set(azieditHndl,'Enable','off‘) 
set(eleditHndl,'Enable','off") 
set(cscalslidHndl,'Enable','off') 


cla; % Create the movie 
FieldMovi = moviein(Mt); 
set(gca,'NextPlot',replacechildren') 
for nt = 1:Mt 
cla: 
if stremp(filetype,'Eyt') 
mesh(xdim,ydim,Eyt(:,:.nt)) 
elseif stremp(filetype,'Ezt') 
mesh(xdim,ydım,Ezt(:,:,nt)) 
else 
mesh(xdim,ydim,Hxt(:,:,nt)) 
end 
view([azi el]), axis(vaxis) 
caxıs([-cscale*vmax cscale*vmax]) 
xlabel([xlabl,' (meters) ]), ylabel([ylabl,' (meters)']) 
zlabel([zlabl]), title(['Frame Time: '‚num2str(t(nt)),' nSec']) 
FieldMovi(:,nt) = getframe; 
end 
title( '), 
movie(FieldMovi, 1,fps); 
while | 
mulligan = questdlg('Would vou like to view the movie again, ... 
‘Instant Replay’,'No'); 
if stremp(mulligan,'No'), break, end 
movie(FieldMovi,1 fps); 
end 


rebuild = questdlg( Would vou like to adjust any of the movie parameters and build it 
again?,'Rebuild?','No'); 
cla; set(axHndl,'Visible','off); 
set(fpspopupHndl,' Value',1) 
set(frampopupHndl,'Value',1) 
if stremp(rebuild,'No'), break, end 
end % Now save the movie for replay later 
[newmoviename, newmoviepath] = uiputfile('*.mat', 'Save Movie As’); 


save(newmoviename, FieldMovr','vaxıs','xlabl','vlabl','zlabl','cmatrix','zscale','cscale','vmax’,'azi', 
el”) 
clear Eyt Ezt Hxt FieldMovi 
eval( nfldvis("reset")); 
end 


0/0 6/0 /09/09/09/09/00/09/09/09/o9/09/09/09/09/09/09/09/09/o9/0 0/69/00 6/0 9/09/69/09/09/09/09/09/00/09/09/09/69/9/69/69/69/99/99/99/o 
0/0/06 0/09/o9/09/o/o/09/09/o9/o9/09/09/09/69/09/09/09/09/09/09/69/09/09/09/69/09/69/09/69/09/69/69/6 0/0 0/69/69/69/69/69/ 
case 'plotEy ]' 
LIVIANA 
0/9 /69/09/09/09/09/09/09/09/09/6 9/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/o 

axHndl = sca; 

figNumber - gcf; 

hndlList = get(figNumber,'UserData'); 

startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 

continuHndl-hndlList(5); edittextHndl-hndlList(8); plotEy1-hndlList(25); 

cbaraxes=hndlList(34); 

set(startHndl,'Enable','off') 

set(demoHndl,' Enable','off) 

set(movieHndl,'Enable' 'off') 


% Plot the "Impulse Response” y-component of the E-Field tin the frequency domain: 

load Eyf, load Hdr 

set(edittextHndl,'String',char('This is the frequency-domain representation of the E- 
Field.',’Press the "Continue" button once you are through with this plot...')) 

cla; set(axHndl,'Visible','on') 

mesh(f,y,squeeze(abs(Eyf(rowposit,colposit,1:Nf)))) 

title(['Near-Field |E y(y,f)| for x0 m; z-0 m; , IDinfo]) 

ylabel('y (meters), xlabel('Freq (MHz)), zlabel(|E v| (V/(m-Hz))) 

axis([0 max(f) 0 max(y) 0 1.1*max(max(squeeze(abs(Eyf(rowposit,colposit, 1:Nf)))))]) 

rotate3d on 

colorbar(cbaraxes); 

set(continuHndl.'Enable",'on') 

waltfor(continuHndl,'Value") 

set(continuHndl 'Value',1,'Enable','off') 

rotate3d off 

set(axHndl,'Visible','off') 
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axes(axHndl) 
clear Evf 
eval(nfldvis("reset")'); 


0/6/69/69/69/69/09/99/69/89/09/99/69/09/09/69/09/09/09/09/69/69/09/69/69/09/69/09/09/09/09/99/09/09/69/99/69/99/09/99/99/909/99/69/6 9/99/o 
0/69/69/69/99/99/99/99/09/69/09/09/69/69/09/09/09/09/69/09/69/09/09/69/69/09/09/09/69/09/09/09/09/69/09/09/69/69/09/99/99/99/09/0 
case 'plotEz1' 
0/99/99/69/99/99/09/99/09/o9/09/o9/o9/09/69/09/09/o9/o9/o0/o9/09/09/09/09/09/09/09/o9/o9/09/09/09/09/o9/09/09/09/09/09/09/909/09/99/99/99/o 
0/69/09/09/09/09/o9/09/09/0 9/o9/09/09/09/09/09/09/09/09/09/69/09/09/09/69/69/09/0 

axHndl = gca; 

figNumber = gcf; 

hndlList = get(figNumber,'UserData'); 

startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 

continuHndl-hndlList(5); edittextHndl-hndlList(8); plotEzl-hndlList(26): 

cbaraxes=hndlList(34),; 

set(startHndl,'Enable','off') 

set(demoHndl,'Enable’, off’) 

set(movieHndl,'Enable','off‘) 


% Plot the "Impulse Response" z-component of the E-Field in the frequency domain: 

load Ezf, load Hdr 

set(edittextHndl, String',char( This is the frequency-domain representation of the E- 
Field.','Press the "Continue" button once vou are through with this plot...')) 

cla; set(axHndl.'V isible','on') 

mesh(f,y,squeeze(abs(Ezf(rowposit,colposit, 1:Nf)))) 

title(['Near-Field |E. z(y.f)| for x20 m: z-0 m; '.IDinfo]) 

ylabel('v (meters)'), xlabel('Freq (MHz)'), zlabel(|E z| (V/(m-Hz))') 

axis([0 max(f) 0 max(y) 0 1.1 *max(max(squeeze(abs( Ezf(rowposit,colposit, 1 :Nf)))))]) 

rotate3d on 

colorbar(cbaraxes); 

set(continuHndl,' Enable','on') 

waitfor(continuHndl,' Value") 

set(continuHndl,'Value',1,'Enable','off') 

rotate3d off 

set(axHndl,' Visible','off) 

axes(axHndl) 

clear Ezf 

eval(nfldvis("reset") ); 


0/59/69/99/09/99/69/99/09/99/09/09/09/60/09/69/09/09/09/09/09/09/09/09/69/09/o9/09/69/09/69/99/99/99/69/09/69/09/99/69/09/69/69/99/09/09/o 
0/09/o9/09/09/09/o9/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/o9/09/09/09/09/09/o 
case 'plotHx1' 
0/5 9/09/99/69/99/99/69/69/99/69/69/09/99/0/99/09/09/a9/9/09/09/09/09/09/09/09/99/69/69/69/69/99/09/99/a9/a 9/6 9/6 9/69/99/99/59/59/69/99/« 
0/59/99/69/99/69/69/69/69/09/o9/09/69/09/99/09/09/09/09/09/09/09/09/09/09/99/0 9/0 

axHndl = gca; 

figNumber = gcf; 
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hndlList = get(figNumber,'UserData’); 

startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 
continuHndl-hndlList(5); edittextHndl-hndlList(8); plotHx1-hndlList(27); 
cbaraxes-hndlList(34); 

set(startHndl,'Enable','off') 

set(demoHndl,'Enable','off) 

set(movieHndl,'Enable',' off) 


? Plot the "Impulse Response" x-component of the H-Field in the frequency domain: 

load Hxf, load Hdr 

set(edittextHndl,'String’,char(‘This is the frequency-domain representation of the H- 
Field.','Press the "Continue" button once vou are through with this plot...")) 

cla; set(axHndl,'V isible','on") 

mesh(f,y,squeeze(abs(Hxf(rowposit,colposit, | ND) 

title(['Near-Field IH_x(y,f)| for x=0 m; z=0 m; '‚IDinfo]) 

ylabel('y (meters)'), xlabel('Freq (MHz)'), zlabel('\H x] (A/(m- -Hz))) 

axis([0 max(f) 0 max(y) 0 1.1*max(max(squeeze(abs(Hxf(rowposit,colposit, 1 :Nf)))))]) 

rotate3d on 

colorbar(cbaraxes); 

set(continuHndl,'Enable','on"). 

waitfor(continuHndl,' Value") 

set(continuHndl, Value',1,'Enable','off) 

rotate3d off 

set(axHndl,'Visible','off") 

axes(axHndl) 

clear Hxf 

eval('nfldvis("reset'')'); 


0/0 9/99/99/99/99/99/69/99/o9/99/o9/o 0° YY 0 %0% %% o» O / 0/5 99/99/o9/o9/o % 0% %%%%%%%%MM%%%V% 
Y WWW / /99/99/09/o9/o9/09/69/o9/ 0%%%%% "o Y % %MM%MM%%%% 09/99/99/09/99/09/99/69/69/9 9/9 9/a9/69/6 
case 'plotVint' 
0/09 99/99/99/09/99/99/69/99/69/69/99/69/5 9/99/99/69/99/99/o9/99/o 0/0 /99/99/o9/99/69/99/69/69/9 9/59/69 9/99/6 % HNy% VA WA WL WALA 
WLWSPVVV%%V%%V%%%%%%V%D%%%%V%O%M% 
axHndl = gca; 
figNumber = gcf; 
hndlList = get(figNum ber,’ User Data’); 
startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 
continuHndl-hndlList(5); DE plotVint-hndlList(28); 
set(startHndl,'Enable','off') 
set(demoHndl,'Enable' 'off”) 
set(movieHndl,'Enable','off') 


% Plot the Gaussian Impulse Input Voltage, Vin(t): 
load Hdr 
cla reset; set(axHndl,'Visible','on') 
plot(t, Vint), title(['Antenna Input Voltage V (inj(t): '.I[Dinfo]); 
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xlabel( Time (uSec)), ylabel( Volts), axis([0 T min(Vint)-0.1 max(Vint)*0.1]) 

set(continuHndl,'Enable','on') 

set(edittextHndl,'String',char(' This is a plot of the Gaussian impulse input voltage, 
Vin(t).','Press the "Continue" button once you are through with this plot...")) 

waitfor(continuHndl,'Value') 

set(continuHndl,'Value',1,'Enable','off‘) 

set(axHndl,'Visible','off‘) 

eval('nfldvis("reset")'); 


0/59/99/9/09/69/99/99/69/69/09/69/59/69/09/69/69/69/09/99/69/a9/o9/09/99/99/09/90/99/o9/90/o 9/o9/o9/o9/09/09/09/09/o9/99/9/o9/0 9/o9/69/o 
0/69/69/99/09/69/09/09/09/09/09/09/09/09/09/09/09/99/09/o 9/o9/o 9/00/09/09/09/09/09/09/09/09/09/09/09/09/09/o9/99/09/99/99/09/o9/o 
case 'plotVinf 
0/5 0/99/69/99/99/60/99/69/69/09/99/69/99/69/69/69/69/69/99/69/69/69/69/69/99/99/69/60/69/69/69/69/59/99/9/69/99/99/99/69/69/59/59/9/9/0 
0/50/59/5 9/99/69/69/99/99/69/99/oS/o9/o9/99/99/o 9/69/o9/o 9/o 9/99/69/69/99/99/99/o 
axHndl = gca; 
figNumber 7 gcf; 
hndlList 7 get(figNumber,' UserData'); 
startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 
continuHndl-hndlList(5); edittextHndl-hndlList(8); plotVinf-hndlList(29); 
set(startHndl,'Enable', off) 
set(demoHndl,'Enable','off') 
set(movieHndl, Enable','off)) 


% Plot |Vin(£).. the Fourier Transform of Vin(t): 
load Hdr 
cla reset; set(axHndl,'Visible','on') 
plot(fp, Vif), title(['Antenna Input Voltage |V 1inj(f):: '.JDinfo]) 
xlabel('Freq (MHz)), axis([0 f(Nf) 0 max(Vif)*0.1*max(Vif)]) 
set(continuHndl,' Enable','on') 
set(edittextHndl,'String',char(‘This is a plot of |Vin(f)|, the frequency-domain representation of 
Vin(t).','Press the "Continue" button once you are through with this plot...')) 
waitfor(continuHndl, Value") 
set(continuHndl,'Value',l,'Enable'.'off) 
set(axHndl,'Visible','off') 
eval( nfldvis("reset")); 


0/5 0/99/59/99/99/59/69/69/69/69/99/09/99/99/09/09/69/69/99/69/69/09/69/09/69/99/69/99/99/69/59/99/69/6 9/9 9/99/69/69/69/9 9/995 59/696 9/q 
0/09/69/09/0/09/09/69/9/99/09/69/09/09/09/0'/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/69/69/09/o9/0 9/9/00 
case 'plotEy2' 
0/09/09/69/09/09/99/09/09/o9/09/09/09/09/09/09/09/o9/09/o9/09/09/09/09/09/o9/00/09/09/09/o9/o9/o9/o9/0 /09/o9/09/09/09/o9/09/09/o9/o9/o9/o 
0/09/09/09/09/99/99/09/09/09/09/09/09/09/09/09/0 9/0 /00/0 9/0 /09/09/09/0 00/0 
axHndl = gca; 
figNumber = gcf; 
hndlList = get(figNumber,'UserData’); 
startHndl-hndlList(2): demoHndl-hndlList(3); movieHndl-hndlList(4); 
continuHndI-hndlList(5); edittextHndl-hndlList(8); plotEy2-hndlList(30); 
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cbaraxes-hndlList(34); 
set(startHndl,'Enable','off) 
set(demoHndl,'Enable','off') 
set(movieHndl,'Enable','off") 


% Plot the "Filtered" y-component of the E-Field in the frequency domain: 

load Eyff, load Hdr 

set(edittextHndl,'String', The frequency-domain representation of the E-Field due to the 
Gaussian input voltage. Press the "Continue" button once you are through with this plot...') 

cla; set(axHndl,'Visible','on') 

mesh(f,y,squeeze(abs(Eyf(rowposit,colposit, 1 :Nf)))) 

title(['Near-Field |E y(v.£f)| for x20 m: z-0 m; ',.IDinfo]) 

ylabel('y (meters)), xlabel('Freq (MHz)), zlabel("|E v! (V/(m-Hz))) 

axis([0 max(f) 0 max(y) 0 1.1 *max(max(squeeze(abs(Eyf(rowposit,colposit,1:Nf)))))]) 

rotate3d on 

colorbar(cbaraxes); 

set(continuHndl,'Enable','on') 

waitfor(continuHndl,' Value") 

set(continuHndl,' Value’,1,’Enable','off) 

rotate3d off 

set(axHndl,'Visible’,'off) 

axes(axHndl) 

clear Eyf 

eval( nfldvis("reset")); 


0/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/69/09/09/09/09/99/09/09/09/09/09/69/09/09/09/09/09/09/69/09/89/99/909/99/99/99/o 
0/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/99/09/09/09/09/09/o9/09/09/09/09/09/69/09/09/09/69/69/09/09/09/o9/o9/09/o 
case 'plotEz2' 
0/99/99/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/00/09/09/09/09/09/09/69/09/69/69/69/09/09/09/09/09/99/09/09/0 
0/99/09/99/09/09/09/09/09/09/09/09/09/09/09/09/09/09/00/69/09/09/09/09/09/09/09/0 

axHndl = gca; 

figNumber = gcf; 

hndlList = get(figNumber, UserData'); 

startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 

continuHndI=hndlList(5); edittextHndl-hndlList(8); plotEz2-hndlList(35); 

cbaraxes=hndlList(34); 

set(startHndl,'Enable','off') 

set(demoHndl, Enable','off') 

set(movieHndl,'Enable'.'off') 


% Plot the "Filtered" z-component of the E-Field in the frequency domain: 
load Ezff, load Hdr | 
set(edittextHnd]l,'String','The frequency-domain representation of the E-Field due to the 
Gaussian input voltage. Press the "Continue" button once vou are through with this plot...) 
cla; set(axHndl,' Visible';'on') 
mesh(f,y,squeeze(abs(Ezf(rowposit,colposit, | : Nf)))) 
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title(['Near-Field |E z(v.f)! for x20 m: z-0 m; ,IDinfo]) 

ylabel('y (meters)), xlabel(Freq (MHz), zlabel(|E z| (V/(m-Hz))) 
axis([0 max(f) 0 max(y) 0 1.1 *max(max(squeeze(abs(Ezf(rowposit,colposit, 1 :Nf)))))]) 
rotate3d on 

colorbar(cbaraxes); 

set(continuHndl,'Enable''on') 

waitfor(continuHndl, Value") 

set(continuHndl,'Value',1,'Enable','off') 

rotate3d off 

set(axHndl,'Visible','off) 

axes(axHndl) 

clear Ezf 

eval('nfldvis(" reset" )); 


0/50/69/99/99/99/99/09/99/99/69/99/09/09/09/09/09/09/69/09/69/09/09/09/09/09/69/69/09/99/69/09/99/69/69/69/69/69/69/99/69/99/99/99/6 9/69/o 
0/09/09/09/09/09/99/09/o9/09/69/09/09/09/09/09/09/09/09/09/09/09/09/0 9/0 9/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/o 
case 'plotHx2' 
0/99/99/09/09/09/09/09/09/09/09/09/09/09/09/09/69/o9/09/09/09/09/09/09/09/09/69/09/09/09/09/09/o9/09/69/09/9/09/o9/09/09/09/69/09/09/o9/o 
0/69/99/99/69/09/60/09/09/6 9/69/69/09/09/09/09/69/09/09/09/09/99/99/09/09/09/09/o 

axHndl = gca; 

figNumber = gcf; 

hndlList = get(figNumber,'UserData’'); 

startHndl-hndlList(2); demoHndl-hndlList(3); movieHndl-hndlList(4); 

continuHndlI-hndlList(5); edittextHndl-hndlList(8); plotHx2-hndlList(32); 

cbaraxes=hndlList(34); 

set(startHndl,'Enable','off‘) 

set(demoHndl,'Enable','off') 

set(movieHndl,'Enable’,'off) 


% Plot the "Filtered" x-component of the H-Field in the frequency domain: 

load Hxff, load Hdr 

set(edittextHndl,'String',"The frequency-domain representation of the H-Field due to the 
Gaussian input voltage. Press the "Continue" button once you are through with this plot...") 

cla; set(axHndl,'Visible','on') 

mesh(f,y.squeeze(abs(Hxf(rowposit,colposit, 1 :Nf)))) 

title(["Near-Field |H_x(y,f)| for x=0 m; z=0 m; '‚IDinfo]) 

ylabel('v (meters), xlabel('Freq (MHz)), zlabel(]H x| (A/(m-Hz))) 

axis([0 max(f) 0 max(y) 0 1.1 *max(max(squeeze(abs(Hxf(rowposit,colposit,1:Nf))))))) 

rotate3d on 

colorbar(cbaraxes); 

set(continuHndl,'Enable','on') 

waitfor(continuHndl, Value") 

set(continuHndl,' Value',1,'Enable','off‘) 

rotate3d off 

set(axHndl,' Visible’, off) 

axes(axHndl) 
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clear Hxf 
eval( nfldvis("reset")); 


0/0 0/00/09/0/0 /09/09/09/09/0 0/0 /09/09/6 /09/09/09/00/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/69/09/09/09/59/09/09/99/99/9/0 
0/00/09 09/09/09 /0/09/ /6/69/o/o/69/99/9/0/09/09/09/99/09/09/09/9/69/09/69/09/09/09/09/09/69/09/09/09/09/09/69/09/09/ 
case 'save' 
7/0 /o/09/0 9/0/0000 /09/09/09/09/09/09/07/09/09/09/09/09/09/09/09/0/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/69/69/09/9/6 
0/90/09/09/09/09/09/0/09/09/09/09/09/09/09/69/00/09/09/09/09/09/09/09/09/09/69/09/69/69/6 

[X,map]-capture( 1); 

plotname = inputdlg( ('Enter the filename under which you want to save this figure:'?, ... 

Eile Damen): 
plotname = char(plotname); 
imwrite(X,map,[plotname,'.tif']) 


0/99/09/09/09/09/09/09/69/09/09/09/09/09/09/69/69/09/09/09/09/09/09/09/69/09/09/09/09/09/69/09/09/09/09/09/09/69/09/69/69/99/99/99/99/09/0 
0/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/69/09/69/09/09/09/69/09/09/09/09/69/09/09/99/0 
case 'printaxis]' 
0/99/09/09/09/69/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/909/09/09/09/09/00/69/69/09/09/09/09/99/09/09/09/09/99/99/09/99/99/0 
0/99/99/09/99/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/0 
Data-1:64;Data-(Data'*Data)/64; 
h=msgbox(['Sorry, this function is not supported by MATLAB 5.1.'],'Reserved for Future 
Release", 'custom',Data,hot(64)); 


0/69/99/99/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/69/69/09/09/09/69/09/69/09/69/69/0 9/99/69 9/6 
0/69/69/69/69/09/69/09/69/09/09/69/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/69/09/09/09/09/09/6 
case 'printallaxes' 
0/69/69/09/69/69/99/09/09/09/09/09/09/69/69/09/09/09/09/09/09/09/09/09/09/09/09/09/60/09/69/09/09/09/09/09/69/69/09/69/69/09/09/99/69/69/0 
0/99/99/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/o 
Data=1:64;Data=(Data'*Data)/64; 
h=msgbox(['Sorry. this function is not supported by MATLAB 5.1 .'] Reserved for Future 
Release' 'custom',Data,hot(64)); 


0/99/99/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/69/09/09/09/69/69/o 
0/99/99/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/o 
case ‘reset’ 
0/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/0 0/0 9/09/09/0 
0/0 9/6 /6/09/09/69/09/09/09/09/09/09/09/09/69/09/09/09/09/09/09/09/09/09/09/09/09/09/0 
axHndl = gca; 
figNumber = gcf; 
hndlList = get(figNumber,'UserData'); 
infoHndl-hndlList(1); startHndl-hndlList(2); demoHndl-hndlList(3); 
movieHndl-hndlList(4); continuHndl-hndlList(5); resetHndl-hndlList(6); 
closeHndl-hndlList(7); edittextHndl=hndlList(8); grppopupHndl-hndlList(9); 
frampopupHndl-hndlList(10); fpspopupHndl-hndlList(1 1); zscalslidHndl-hndlList(12); 
azieditHndl-hndlList(16); eleditHndl-hndlList(18); cscalslidHndl-hndlList(20); 
plotHndl-hndlList(24); printHndl-hndlList(33); cbaraxes-hndlList(34); 
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9o Reset the Handles & variables to their default values: 

cla; set(axHndl,'Visible','off,"NextPlot',' replace") 

colormap(jet(64)); 

rotate3d off 

set(startHndl,' Enable','on') 

set(demoHndl,'Enable','on') 

set(movieHndl,'Enable''on') 

set(continuHndl,' Value',1,'Enable','off') 

set(grppopupHndl,'Value',l,'Enable','off','String','Near-Field Group") 

set(fpspopupHndl,' Value',1,'Enable','off) 

set(frampopupHndl,' Value',1,'Enable','off,'String',"Total Frames") 

set(edittextHndl,'String';'Look Here for All Commands & Prompts!) 

set(azieditHndl,'String','-37.5') 

set(eleditHndl,'String','30') 

text(0,0.1,'Press the "Start" button to begin the analvsis., ... 
'HorizontalAlignment'’,’center’); 

text(0,0,'Press the "Demo" button to see a demonstration of the software,’, ... 
'HorizontalAlignment’,'center’); 

text(0,-0.1,'Press the "Movie" button to replay a saved movie or create a new one.', ... 
"Horizontal Alignment''center'); 

text(0,-0.2,'Or click on the "Plot" menubar to create various plots using data from your last 

session.', ... 

"Horizontal Alignment''center'); 

text(0,-0.4,' All commands will be issued in the text window at the top of this screen. ... 
Horizontal Alignment','center'); 

axis([-1 ] -1 1]) 

set(plotHndl,'Enable','on') 

axes(cbaraxes), cla 

set(cbaraxes, Visible''off'), axes(axHndl) 

set(figNumber, Visible';'on', UserData',hndlList); 


end 


heal 





APPENDIX B. Theoretical Near E-Field Pattern of a Simple Dipole 


% DipoleEzf.m 

% 

% Caleulates the theoretical z-component of the near eleetrie field of a simple 

% dipole transmitung antenna and eompares it to the same AWAS-ealculated near 
% field. The dipole is assumed to be along the z-axis. and measurements are 

% taken along the y-axis. Also. the eode is optimized for single-frequenev 


% analysis within a single "near-field group", as defined by AWAS. 


clear all 
% **' Readine AWAS Antenna Geometry Data *** 


[geofile,filepath] ^ uigetfile( *.geo0',' AWAS Geometry Files); 

geofid = fopen(geofile,'rt'); 

geohdr = fscanf(geofid,'%76c',1); % Skips the header data on the .GEO file 
gl = fscanf(geofid,%3d',1); % Skips the first integer value 

GndPlane = fscanf(geofid,'%3d',1); % Indicates presence of a Ground Plane 
g2 = fscanf(geofid,'%3d',3); % Skips the next 3 integer values 

nodel = fscanf(geofid,'%14e',3); % X. Y. & Z eoordinates for Node #1 
node2 = fscanf(geofid,'%14e',3); % X. Y. & Z coordinates for Node #2 
node3 = fscanf(geofid,'%14e',3); % X. Y. £ Z eoordinates for Node #5 
*** Reading AWAS Near-Field Data *** 


/ 
0 


«o 


[filename,filepath] = uigetfile(’* nfl’, AWAS Near-Field Files’); 
fid = fopen(filename, rt’); 


Ne 7 fscanf(fid,963d',1); % Number of Position Groups 
Fp = fscanf(fid, %12e',2); Zot 2) Stary Stop Freq Pts 
Ns = fscanf(fid,'%4d', 1); % Number of AWAS Freq Steps 


fl =Fp(1/le6; f2=Fp(2)/le6; df=(f2-f1)/Ns; f= 0:df:f2; 
AntennaType = fscanf(fid,'%3d',1); % The Antenna Type (Transmitter/Revr) 
Nd = fscanf(fid,'%3d',2); % Skips the next 3 integer values. 


xp = zeros(Ng,2); yp-zeros(Ng,2) zp = zeros(Ng,2); 


for ne=1:Ng; % Reading Group Coordinates 
xp(ng,:) = fscanf(fid,'%ol 2e’,[1,2]); % (xl. x2) Near Fld Pts 
Nxs(ng) = fscanf(fid,' %4d',1); % No. X-steps 
yp(ng,:) = fscanf(fiıd,'12e',[1,2)); % (y1. y2) Near Fld Pts 
Nys(ng) = fscanf(fid,'%4d',1); % 


o No. y-steps 


zp(ng,:) = fscanf(fid, %12e',[1,2]); % (z1. z2) Near Fld Pts 


Nzs(ng) = fscanf(fid,'%4d',1); % No. z-steps 
end 


NG=l; 

x1 =xp(NG,1); x2 = xp(NG,2); 

ANA 

ZII 22 = 2p Nee) 

Nx = Nxs + ones(size(Nxs)); | Ny = Nys + ones(size(Nys)); 

Nz = Nzs + ones(size(Nzs)); % Total # of x.y & z points 

if Nxs(NG)==0, dx = 0; else dx = (x2-x1)/Nxs(NG); end, x = x1:dx:x2; 
if Nys(NG)==0, dy = 0; else dy = (y2-y1)/Nys(NG); end, y = y1:dy:y2; 
if Nzs(NG)==0, dz = 0; else dz = (z2-z1)/Nzs(NG); end, z = z1:dz:z2; 


if z1==z2, % AWAS near-field data is in the xy-plane 
N1=Nx(NG); Nls=Nxs(NG); N2=Ny(NG); N2s =Nys(NG); 
rowposit = 1:N2; colposit = (N1+1)/2; xdım = x; xlabl = 'x'; ydım = y; 
ylabl = "y"; minxdim = x1; maxxdim = x2; minydim = yl; maxydim = y2; 
else % AWAS near-field data ts in the vz-plane 
N1=Ny(NG); Nls=Nys(NG); N2=Nz(NG); N2s = Nzs(NG); 
rowposit = (N2+1)/2; colposit = 1:N1; xdim = y; xlabl ='y'; ydim = z; 
ylabl = 'z'; minxdim = yl; maxxdim = y2; minydim = zl; maxydim = 22; 
end 
Nf = Ns+2; 9o Number of points to f2 folding freq (needs to be 2^m-1) 
Nrow = 2*Nx.*Ny.*Nz; EH = zeros(6,Nrow(NG)); 
Eyf = zeros(N2,N1,Nf);  1fGndPlane=1 & z1——z2, Eyf 7 zeros(2*N2s* 1,N1,Nf); end 
Ezf=Eyf; Hxf=Ezf; 





for m=2:Nf; % Frequency Stepping 
for n=1:Ng; % Reading Group Data 
S = fgetl(fid); % Read blank line between EH data arrays 
EH = fscanf(fid,'%12e',[6,Nrow(n)]); 
ifn == NG, % Saving Data at Selected Group 


Eyf(1 :N2,:,m) = reshape(EHG,1:2:Nrow(NG)-1)+3*EH(4,1:2:Nrow(NG)-1),N2,N 1); 
Ezf(1:N2,:,m) = reshape(EH(5,1:2:Nrow(NG)-1)+3*EH(6,1:2:Nrow(NG)-1),N2,N1); 
Hxf(1:N2,:.m) = reshape(EH(1,2:2:Nrow(NG))'+j* EH(2,2:2:Nrow(NG))',N2,N 1); 
end 
end 
end 


if GndPlane==1 


Eyf=.5*Eyf; Ezf=.5*Ezf; Hxf=.5*Hxf; % Scale factor to comp for endplane use 
if zl--z2, % Changes the appropriate variables due to 
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N2s=2*N2s; N2=N2s+1;  zl=-z2; % the presence of the ground plane. 
dz=(Z2-zZ1)/N2s; z=zl:dz:z2; rowposit=(N2+1)/2; ydım=z; minydim>z]; 
Eyf((N2+1)/2:N2,:,:) = Eyf(1:((N2+1)/2,:,:); % Reshape the Evf matrix to 
Eyf(1:N25/2,:,:) 7 flipdim(Eyf((N243)2:N2,:,:,1); 96 tit the end plane "image" 
Ezf((N2+1)/2:N2,:,:) = Ezf(1:(N2+1)/2,:,:); % Reshape the Ezf matrix to 
Ezf(1:N2s/2,:,:) = flipdim(Ezf((N2+3)/2:N2,:,:),1); % fit the end plane "image" 
EISEN SE V2END2 Vz HR l(N2 12): % Reshape the Hxf matrix to 
Hxf(1:N2s/2,:,:) = flipdim(Hxf((N2+3)/2:N2,:,:),1); % fit the gnd plane "image" 
end 

end 

fclose(fid);  fclose(geofid); clear EH 


% *** Reading AWAS Input Impedance Data *** 


[filename, filepath] 7 uigetfile( *.zpm',' AWAS Z Parameter Files); 

Zfid = fopen(filename, rt); 

Np = fscanf(Zfid,'%03d',1); | 9o Number of Ports 

Fp = fscanf(Zfid,%12e',2); %(fl. (2) Start/Stop Freq Pts 

Ns = fscanf(Zfid, vo4d',1); | 9o Number of AWAS Freq Steps 

zp = zeros(Ns+1,2); 

for ns=1:Ns+1; % Reading Group Coordinates 
zp(ns,:) 7 fscanf(Zfid, 9612e',[1,2]); 96 (Zreal. Zimag) Zin points 

end 

Zp = ZEW 1290). 

if length(Zp)>1, Zp =[1;Zp]; end 


EZ = Zp EZ 22) % The final AWAS Near-Field Values 


“o Calculating the theoretical near-field values of Ez 
c=3e8; lin=1; 
f 7 input( Enter the frequency (MHz): '); 
=2e6* pi*f./c; 
% L = input(Enter the dipole length (meters): '): 
L=2.06156; 
H=L/2; kH=k*H; coskH=cos(kH); sinkH=sin(kH); 
% zpt = input('Enter fixed z for calculation (meters): '): 
zpt — 0: 
R1=sqrt((zpt-H)2 + y'.52); kR1=R1*k; 
R2-sgqrit((zpt-H) 2 + y 02); KR2=R2*k; 
R =sqrt(zpt^2 + y'.^2); kR=R*k; 
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Ez2-repmat((-j *30*Tin./sinkH),length(y).1).*(exp(-J "KR 1)./repmat(R1,1,length(f)) - 
exp(-]"KR2)./repmat(R2, 1 ,length(f)) - 2*repmat(coskH,length(y),1).*exp(- 
J*kR)./repmat(R,1,length(f))); 

figure, plot(y,abs(Ez2(:,1)),r--,y,abs(Ez)); 

xlabel('y (meters)”), title(['Near-Field |Ez] for L=',num2str(L),... 

' Dipole at z- ',num2str(zpt),' m & f- ',num2str(f), MHz ]) 

legend( Theoretical Values',' AWAS Calculations") 

figure, plot(y,real(Ez2(:,1)),'r--',y,real(Ez)) 

xlabel('y (meters)’), title({"Near-Field Re(Ez) for L= ';num2str(L).... 

' Dipole at z= ',num2str(zpt),' m & f- ',num2str(f), MHz]) 

legend( Theoretical Values, AWAS Calculations 

figure, plot(y,imag(Ez2(:,1)),'r--',y,imag(Ez)) 

xlabel('y (meters)”), title(['Near-Field Im(Ez) for L= '¡num2str(L).... 

' Dipole at z=',num2str(zpt), m & f- ')num2sr(f; MHz] 
legend( Theoretical Values',' AWAS Calculations) 


126 


u) 


Un 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center ....................... eene 


8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, VA 22060-6218 


Brillo A 


Naval Postgraduate School 
411 Dyer Rd. 
Monterey, CA 93943-5101 


Director, Training and Education sas meneren rir 


MCCDC, Code C46 
1019 Elliot Road 
Quantico, VA 22134-5027 


DPitector Marine Corps Resear M OCeEnter amire een 


MCCDC, Code C40RC 
2040 Broadway Street 
Quantico, VA 22134- 


Director, Studies and Analysis Divisionen... rasen 


MCCDC, Code C45 
3300 Russell Road 
Quantico, VA 22134-5130 


Marine Corps Tactical Systems Support Activity......................... 


Technical Advisory Branch 

Attn: Maj J.C. Cummiskey 

Box 555171 

Camp Pendleton, CA 92055-5080 


Manne Corps NeDrtesenfdlve o ne E 


Naval Postgraduate School 
Code 037, Bldg. 234, HA-220 
699 Dyer Road 

Monterey, CA 93940 


12% 


No. Copies 


.06..00000000000000000. A 


Chairman, Code EC uuu RU T EEMMME 
Department of Electrical and Computer Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5121] 


Professor Michael A. Morgan., Code EC MV 
Department of Electrical and Computer Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5121 


Professor R. Clark Robertson, Code EC/RG.....— on ea 
Department of Electrical and Computer Engineering 

Naval Postgraduate School 

Monterey, CA 93943-5121 


Maj David V. Adamiak, Code EC/Ad....... 5 "TTE 
Department of Electrical and Computer Engineering 

Naval Postgraduate School 

Monterey, CA 93943-512] 


128 








p PTA V 
6/02 22521-200 "* 











En 




























































































































































































































































E 
Lr E 
La 
L pa 
b odd .. E .2 E 
E = m - E "I » x 
e -æ senu sao - PR = E = rS 
= E . € . ^ a . =- E eo. 
rm m. Dr c "- m E "ED LE b . "m 
- L Lo - ? EERE 5 . E i .*^ Hr Į E O . 
= E -- = o ee "c 
A LN - se. A A IS "1 H E m y E - o. B - - 
^ A = “.. , E E E 2 > .- z " PES 
- .» . E . . - q LI ” . "] 
« De er =P” LE . - Ld LI FLELI 1 e . . . 
e. E ~. ax E Lem - .. - "s g E 
- Mn LEGE en Ld 2 Ld Ld LI LJ Ll sr 
ae - ES . LES P o^. -.. n . PS 
- a OA o E a o - P E 
-- - "- . e .._ LAT. i » LED LINE ` 
E Lr - >» -x = . `% o» L E - . . 
P LM P e Ll . . LI ."9 "- - E We * = a P 
e a C . Py » . x LE E LA - " LINE 
rar pe n » AS IT E LES ms .. - "ID 
= - - - p O - ES -s - E L LL" - "T 
c T o. E ma... - LITE T E E = "€, M " a 
“oo. E E -. E ~~ = © - .. E e o cl LACUS ES A g . 
TE - "S ~ " LN lol LLL LL T a f=- LEM D "n .. z E ry - mo. 
y E E - m es ns US. 0 gu EM - E 3 as Son D u a 
- E > » LI " B ee = .- Li = > ..... " M € v aT 
mE > r L „u -.. .. "a^ E = “eo. = A 
EE -- - P - 1 ss~ DL Er AA E 
E PES - 7] a 0 mo 5 er ys .. " * — 9. wu 
Ld » A mm mo = 0 - ee LE ` rd == . . 
E Eee een lm LP E - ET E > 
"Nt ..o a p Le . .. ELM P e ~ a E "n M e = 
- ama „me... - L] P" O ec ey D u... = .. "n -~ " 
mw a^ == am E OT Cees .- ... A .. LM CE rs 
a LE LIS LI "2 "- o» - o ape - ==st [IX LIE En to Ly LI 
rm .. POL - eo. m “=u "i9 " E T a -. LI mo =; e. .. - r 
>e- a LM ee P *- -- en nn H a 
= = E P A AA - . eae = Ji .. E > e 
sure D - ”. 
x LI E »oo.4 -. 
* H i^ 
* * Bee, 2 
s E] E ^ > E 
AO etd "uu pS A ds <<) ae 
A A ODE On r . = . u .- ` 
ECT .- .r E ze = pa E " er, 
we fF LLLA x ur o LLL "I 
2 =P = AR o. <. Mae = TX E num og 
pw P RESULT LE S B ^ SO m CE E e ta" 
"9.5 REL 90 -..-.. E % G . " E "E m T D . 5 A oo "cet a 
a... LI + ae Om E =" es EI . PE L] [nm f = qe... CINES -— = 7 ae ae he. 
M ES ML PET e See aee hn: er Pe -ron Er i" NU A AA de > A km M m FL ae T 
-— v. LII NEN A ey P 7 LII IT € m I LE. m an E D "P OPOR pv . 0» A PE A . "P "e. 
P "ns EE e. Ey we P ve .. "E - Fr m [EE PES C | LEM Li IT u. 
P "o g^ e-q na 2. o- o roa p- [TRA s » m E LI ^w ws oy LN T a s ELI A = .. "earn Los 
= E ~ LI . E aha ` » ner IM LT LI $e»25 979 . - eo. de eng ^ bor be 
. "- P" . IL P ve ql Uy .. DIDI TI - qn AND „un nun nA p^ L- 
en E Y pr pm " A mm Dd UA M yrsa "59  -. Fh! aaa. MNAE e A ap sagose tug "o 9 
= wry. © LU oe . "^. A "^ M o Li P q. MINE LII P LEES - u ee LII E MM LIN ded - 
Lob LUPO Pe tare, - a -. LM r ev = ar LI en E LJ -.. e. - Ye 2.1 w 5 
e A A TA TS RS = e EO EI a PI A " en De ees TT AE ten 8 
MW a wae dee FE A LS » Le a "TI "Eee rent m > et 
A A ie Are Beten np Lr = Loa . A PU u 
bel et Pe Ll O AAA te UX ry - "e Oe ae ee % E =o mw 5 sw. unter Ld 
AA AE ARA TA AA e A " a. a alt - » A . } .. 
A A RS o LPS LN PI nen n = Lux ni S Seale oe LENA mr cra Y a 
ao M Pvt DLL a 2 099 MD on "o. ^ ^ 32 E OO ec a > m La... sor y o... PO P D PP P. M PM DELE 
Pen = ed . os - yu » ~ T E LET y P "m LITE J SAS A “e 
- PE DEL ES WM oath n un. ere ele x " Hnnc d 
E G " E E Lm T LN LETT: P r _ we a e 
> $ LU METER TR x Do 2 T" er 2 me y A A < rm en 0 
" ^9 . es Ths E ange Bey DET " LET JL PP ... ... 4 ye E " ra e . wu. 
ER Pr k a o O E " wwd- " It og ° aso - AT 5 T E ia BEE 3 7 BE Be 
SH Wir Er LET = ep rl T .- D mr A m E mé. a ee ` 7 ow Im À La zu Lm DN Ar 
- na o Em TITS da Y qx yo em. y "mo Fe ” +. LR " 7 an y. w^ y P P LIP DA ey LI M M WM A A 
Pé A A ON A ae F " ee ee "^ 9 5 ore D COM LE LR P bed ILU T 
a tell ES e. » L a Ao . PV AA a T © sewr PAN a a a TP klein, 
SEL Pm LER Ap ed hen a L] e... ya m a * P » “a . 8 . - iilis LAE LUCI NL nn — "ipd idi 122 Y A 2 0, 
e CE ma e "m "v. m nn *" 9^» e Mees , Ep = a... .. 9*5 s PER ^ * wa ^A. 9 vm - e. bh 
O NN NAA s q : RETR dumm o LET AAA A as eT AIT wu LN | 
Lal" ... vat Pop =P papp = E. E m EE Ba a LE a ET TR Gr Ve oe LIN I — ] “e 
ZI o o LAA o A ene LM V LIC Nr 3 s 2 p "n I PLE E [o PRAET Medios Loop P MS E A T E 
ku (mm u e In I LIT att Mis nd TE FC Sou =, oe LAT Y A YT A na " LE IA u LEX 3 dd n NN ^p* o» 
ot) an ee ’ Be E rM LM we SL a Pr pe w wg c einai bate telat ht E E ^ev v» E y 
nd a bibis E 20.7 ^ * EL 03 7 Lesbo WE £ o Pi - na as s * m LIN 7 EJ Sehr way 1 Au LIS Lh SS x 
Du nm | eu. 212 ^ e q A B O A AA ey “\ Fests we ves ML en . 0. ahd M „Rn 
BL aiia Ar er nn AAA LE D e n ptas Puede " Ep ar Eae em Jia Lo A pepe phar. GRANDI I 
aule el etis, MEL E ER Du 2 EL SLT t eset E A A e e P PART fu > uer E P CE e. LIS a... . A AA 
A d. - beled bile eT oe eee RIP NP LLL pu - „. a = »* 5 veaa ^ » > "EY. oe P a IL y ALES 2 0*9 — Aem - 
A E — e ú e ce t A EA oe re c. Ll TM LI yn E M he ET EA af .. » USA E LP 3 L albo LO T *-.9 9 £z A PA 
m... —— nn un Ir a LM m aT LM LP e o. EA O se gt AS Orr Om ae Lal. ML E Momm, wien wenn mores 
"im. of = DIM E CTET E n TM "e bI y "MES "e "Een DS E TTA S 13 rn A TS P gar A amd nn 19 
a ps diu, -—— ...ú-» a I ET = ner Dur TE Tan m on „eo. E - EI 999,29 .a EL Tr e a E E PA A A on lacs Lr TI mn He m on 
bue a aa a —— A LL: E" M S nn PT m E uM m . eg 0. 0 E -— T AA weeny E A ET ^ PLI LENLI MS De AAA U EAN ap ST 
u en ——— Lei Lr mad PX” e EA oT ATT A A v.n. e e. >» os v A a, as . 00, er oo ER 5 A Se IA MA 8, eas Cott Re Pe een LEID 
bad LIE ELT E— — — — 4 O) sarr ar; a vaso a uss A M T LIN M" -s. * 1988 - - * “<A 2 ., 5 A se «* "o. s MEAT Ere eS IA TAIANA t PAN e rl Ld iid] 
A SU AI oh, Let ere P ch ae . E A A O A AS Py MI ZIEL OE "Us ELT, r T aba LR nn na en NER EP UE ULL Ta 
utr tm => er Te Eee e pa w^ rn. s OTT E P to g s a... wen rt u m rn ns neue . A ee Pe we am MIRI 
I ———— un Bene sr D ^s ?— "91536. 9.99* 9 LED "E WIS O A "m EN ed A TS A P m A m dotar LAT VI u C be ad od ee ee A ae = see - . um 9 en d). 
rne Am ^ — PH€— ON A O O ' IIA -- n ee [LENT M CE PT" LIB c 100 q mo eer ef tenes mPa e LLL PREIS er “e. em 
A aR c= > y ae?) AI O EN er eee E LS . .. E ET RO A Lor» AA Aa. que 2 8 LII Ars c PI gunt ge A Panem - 
wu... A en. Tr . EE e " AT PO e a ee ^ HIT LITT o m P P PIT dis es^ hte ETA E o pide. "ur" PUN FUEL VA o e a rs lA Ii irre in 
LE Sy Se AS XL FS A ARA bed, E vob , Lao equ PS E DADA EE n E (T EM EN id LI NS . & a" -» 5a wa LU" n Su A E A RT A -" * - nn en ra I 
rr O ble sva» : Ra TIL ed NT m LT y má n POE “u MB ie P - e E AN A AT AA E er te a... ) A ee A CU en 
LE —— LI LECT ... tens ge mei, "T A PLA * ras LP - " nn. oes # Upon 9? m^». E "T .. Don MM e -9 *P "*-*-» ab e, a A TS Pince LI ME IT S bid LL IE 
T) ——. co v LS p ELIT. UaTÜ?W ue C, or o CENE CELL [D E es y o wer i Ty ` 1% A . EANDEM MM. "wot p MER ble. beak oh en oe > 000/00, se 
e "ase "Er LT "n ry et. $ RON Un 000 m. .- wow we $* 13 ALET Lehre hy Ce PED Lia T ELS on ger bd ee ah dad 
C ME bees dod 0 gen un. BO rs [t tT = EI e... e T lio LLL PPS ILLI P Qa "I9 AR e oe9 r^ aues o5. 568 mati A eth dh lean oe 
-—- — —9 O RE TETT ee E A T eure LS men E or -* rn... ..os $ Morin dett MO nen E e Ld io M ee ores 
== UM dioc M^ "m nn RT RENT" AREA MER io y a A aL me A Lr 
A —n <T Id SE Zur DE SE DE Te MILI A A eee os IT A E - E PT ee Lv E $ Ams & OA IP  T - 299099929, 1l Cet aea. sua. a m rh gru, 
— — > LET SE u Pre e EE "IM Ir - *M oz. nme Y» »^» 7 Le Fee D e Tre PP rn. = A cake! abit bitte nT) P m RI 
„= O un, Fon MELLE TIC A pe t HUE P VEM are *- “æ. =u E rv LAN tn een ee a TT 
a LEE N Pin tal th Ne" PL O MO A SEER ea - m 9.5 mem IL 0 s. IER IET S rH 
a O eto Gh Cs 5 9 92.99 Y diese PPM rfe Fairness a Se re Denen Iruya 
ME 3 = Lexi A RO 4 a PM v^ EP enge Sauger id “nn Shen erg wende bett 0 75. NN 
— - a a an A A ee es dti P 4 Z AER LE CET Tr ee ee NELLE TE RET [Qe e... n E LO a os DAA ld a de ee 
— Es FAUNA e A ehe EI OL PL aaa aT pe A 1 O NTE Fo dia, E D AAA 2”. y E IDO LL SS O A A er 
- = CO En antenne dic ED PaL RETE PLI PM ELT RIT MET Sd 9 Lv TID M V FER P DEEP E M w 1> 2 De a TA P gia 
— e ely rta re má es A IE ee A Oi e " LLL SES CN So do Mere Chi - Mes "m TET ts o A PRA 
a) mm = co i-o S MN LT eek Led nu dica =g aaae e MET Pant Test ahi a a 7 IR A TE OO m ON 1 E es my Mri ee Me a d D a ee al 
L| AAA Lm s Cr AP nr serap E A t Aso h te Lan ot E ree Pm "RO A AN o PIO tuu . ames > ehe ht Restate oat) od eat ear itis 
- u bo dare TA A "1 - Pe. V lo y hd o PM Air ETS Cd he ee OS eee DA NP TT o n bi ee Oe 73 LE. $9 eia usn Ld d 26 p 70. UA pe am Lil A TL 
Ga m rv" BIET A A E y um 9*9 x. E “EEE u N ee ee CUA IA et RT msn tan ng; UA 0 ma 
> 0d A nn yo . er) I 40. ele Foe de TR bik Le Yee en tyne UT A T TT a LT FT T ke the le ee ee 
> E O do hel. LIT ord JT dL s. Lll etree ee bie A halen Pen disti Er a a vt am t oc o 
LJ ~ rra aa 25 ILS "a. “we ar > lo” dem o agree ap Ibl LI Ka IL" w""os 0-9 ".« qu mn 
MILITE wann AL fv mos **es et tact PCS ARE N ap obe Me LIE T ee N nm Er ie er TP TEMO 
——— LSe peoa A E AA * "piane me ss O ae %, IA E Dur Lee SAS E Fie Ge = 
— € r ) „2... . E a LIT LIT IIT EI A  * aequa “v ae eT £ q. ar, 
- —— E " 999-65 6. 94*$*8 "aram. x Lus ls, "A CENT) A 6 M LE ET TS .. IS A A AT 
= Saran LU ee tal ed Pld AI DAL PPS [EUN E ML rro o o — IM DIPL UM ere Dur Tu n 
nn, a m MOLLIS wm OA T d TT Ce tet or er : A A Peng CA Rae LINDE ade Y 
— heh LIES aw BIEN TN LA I De FI ET ‘ve Pasta Fi) 
== Jar A s e« Mel dii SUP 4 
—_ - 


Sa YS et 
Selle o LAT. EM T ET f^ 


E 





DL EL SLUT TL ee .1$ 
= 











id LLL on A att nt, | 
ei MILITE. ne en e] En SZ I 
BE O A rr 199659, r^t egi ^n. waap . 
A CIL Ten 8 36* 928 
use LEM etos . 


Eee nie a WE ER 
ze) n 





M NIIS A 
-- 4 Da TL LES 










































































APIO Ps e ato inel o Sur od binae d 13] 
UA UNI FT e t LANI eadeni Ad aL OPI un 
ARA VO LOU un au, damen Pea” aegea, „s tc A S MALI Lore" 
"té*r*9 59 .*8€ aa c, «e Y ^ e E EE OA ANI IIA “e etra Elea ara d 9 ee E 
Bil SLE TL a BILL MAL e. fault Lo ptao més ,9, tur mora A ber a alata hen debili ran I ee rc SS . e w 
FEM X ty a MESI LI T eat atom we es ne eet heed ^ ET aL ee ee te te ds ELE T bad ARES bodie leg Beim 
45959 na d, m COTELAR [TT] A Pen ELIT hd LEA ia TS AA Ah LONA ya Boc A O ERAS ADA Ballo LET 
ee A Ld idc ied hoi Tk E OT . ie a a PEA AT emm LL et a ee a "me ibd Lm qa Led Ad pl de et 
A ee bed toed BT ot e P LN" Too m A RTT p ra us A LET Dee M ws a EI U Tn | bid ola uL er i ihn 
wat E LIT Pda rrr" Bills LL 9h ws es ENT UA sss . mne IS 8 ws a : [DIT Bibi Re ELT BE RN me. A ad eei E ana ek cse etz, wee 
d Cahad ARTEL AT O PR E ET Id TI T lom "EM PL PP IA ¡UU ts DIA ELITE ZT Tr ee RR ihe] Eg ela 
p **a d did bey ts E um. DA EA A DS TE MIT b ehe bec, AE P me LT KIT De a Le She u ed 7 
e E .. E . A CoL PR hdl dnd at eto S E IATA AO Re hd edt Reelin Pn LL Re IT We ge 25596 Se" PLUS vade ema S oL PROP SN 
iere ate eii EL T E re A Bi EL SETS Ettore MUR MET SPetr. Kb et BE Re UTI s a, Pen Ll 
ps. ida 3 DEN LT T LLL Core E taht bata sata A Nara ied T7 CE DL T PON re ehe) Ltd babel 
ee Sep A Poole dad. didi ee EL RG ^ae 5 E eh ee PN 
E [ J LI rue o TP te A TP S rT Do Za PI ES ¿AN aa a hito Là re 
d eae E Porq 00019 $ IS DO MEI PIT up" Meis 0 MUA Pe as Ld 
DEAL TE Mi O Vr BE Lr udi dodi. LN Em 
rer in : PL IM MSS a LP x i 
Al en ER ne ee w ici a en 5 b ato PA 
iip dre deett TDI ARA TTT TS ene 
AA aT = bie ts eee nn tn 






OF Os Ld S NV" 


LIED Um rM Du a 2 


uw 





LALA AAA e a a eek 
CN Y aes EIS! 

ELLE TY 
A AT 


AD 
PLE C LITE ee 















DE LS PELLE pp" 
i 2 br De he ek ee 
Ai IN TA ARO YN. en 
ML TES E Male oll Media TTP hr, CP m" LITT II LIMITI Las Lit eT eee IE PUT s 
d L] abb M LOT Y Seales we von A pe FIL ME UY EP Pr Tm er 
* 9* *"5afóta^ ve Setups marcus P. ok HN P n ET 

CLASE ETT LESERT E Ba A 

Dosto tator eb on tours 


La TOT A 
s. aaa T PT T T a TT 
EA CETA A A 
AAA A ann A AAA PUB nre AS 
TS sung Lic ITI er 79.0 9*45e6* 7:3447 4,90 AOS AS 
AA Gen A OA LI 
"ehe hod hdd Le FR FE EI rr 
LAIT ULIS ELLA LP. LII « EE ATN 
AAA LR ol hd hd 


RII LAE] 






Pee rd 7 
ML T 

LENDERS 

d “avd poto 
e kaea 





"e" 
20? dar 
AAA 





m 


* bla EM ci vor TE Uy 
AMAIA a 


Da gm Sti. gu 








tél x iin eir es 
bd NT Em ed A E 
bi bid Lt cete Depart e 
e O e ak be 
lied teh Met he he et tien NL LLL tl 






Ce ar DU Sr h OA I 
lam ELETT TEF 
wae d E ET 


. € . Sr 






































Mob LLL RETO 
iet ¿TU a 
LII m aroma 4^ "anat “asien ai, baa ILIUM 
Le LLL: P ECET ELA LL Pte LP. LIEU 
e ms N, mr SML Mm re MANLIO wb 
W. 100% ya AS A E IA a 
t OIT id Gransee d Il » Maianen a te 
CO ve wa C VILI PLI u... Pr o SORTE Fr te 
E bled et TTR Ka LU IT 1 SE SI pe TI PPM" Be a I a n LAA a e k aaisa aa pe 
TS A IT NT - ..ws ds A LITE Street at BS et a eel eat eee Re 
beth Lae tT Th TS dat HILL MALE tall Pappe Aiboa La DR LTE EEE TE np PO DE H 
MELLE Te ee AS TE" al hl - bild Ris o ad A E T AA 
ae te oe | A Bl TIAM = A IT ARE Ms - 
O ees LLLI PME TM mu urn, 


bleed bel cont at Ur e 





St tte Lt he LIN M] 
LJ OO S 


LE LM 


b Landa n Te Pe 
bet a Paten het Restos 

yey e MATT PE prs 
ILLI A biet d ALL PT PN y] 
a Lb A ATEN 
p | 


BITS Y 
$ 99*e1089,09 ae 

eh LET TEO 
A II TT jos mos 
DM NIIS Tiree tye ae 
CA I FE wer oe LUE uud. rl 
ALI PTS Mic A OS heer es 
AS ATA TOA O Pr 


FE LIT ETET A 











CIA a A 
m mer IR er en ae de 
TS eit ee ED eee 





. 












E A YU 
MIT PT 
AA A TTA EU 
mea A hee Oe 





I eb lb Li Twp . bidet ca) 
ee eh heh te ET den eae tlt Lhe 

















E MILI b dera tela ee ee 
Lr CE =... Dv ¿ag ra ie, 4. y 
» NAE las p Fer 
ALI LA f E * «32 C39 yd Pese a iL Tre P: Er ridi Quá rot) 
Mis DE LA p AU entries ln AA Pay T TT “eee On * 
HALTE rete Erbe ELEME en LEE s A er ee PIO La ET Qr A. Ld dee nae Ld 
MIB IPMT MA TAO edi . EN la ss Tt) Phe lode EL E T ur 
L LII LEE MP v do. Sra id bela Ld H LPEN E OS AS Dd bade tne bh el ee el LES II Hr 
HE ALIE EM LES se LLLI rn N T w*std jew « S 
iA PLATEIS te A 0 LUE PI “ Pao 19 ^ 
LRIPLE TES LAE INE TI LEITM ROLIT 
f € HV Li 


A a T T 
Mall E ETT 
e Y et 
A? 













heb ek teh re) 
Pp pe lege E Li 
AA 

EIA 
"E ye DE AN 
A Veep bind ow Pee Es | nn ADI EN eT 
LS E KEPET AA LETTY .. AS AA a 
Lies cee Laadi TLA EL] 4 014” Pm. a. 


TA CTO 
odds M LLL T E Aa A » 
POLA yo de a yo ud rr A ee ETA A LIT 
LL MEL Te LT A E Po Sen O ORTEN 
KB LIE wee Le Pa Lt ee ees und CPP PIC A DIL S E 


TA © Se er we je PL IT II Zu DIT LTE RR 
DENM Tnm e LLL ette A 
tera, ANT 


O tee eT 
i ee 


IL SE ey 


Da ee ere rete 
e Se bee the oe 


















dd Like beth La) pe eee ne AAA 
DAA ARI rya Ima ane td Il mr 
LL o a DS TR ed 
TTV m IE "TI DELE, NET TT hell A LARA tid 
u... u. 0 N I n raiat Te TRY MISA 
.»e um T “Palo sun gung TELE Il T Pr ELLE C] dh A LITT mort Ge, i du o EL ee St Td et Lat 
HU PL ECL PSP" Abl 0 LL de dE TE ELLO IS LLLI Teu a zer re Pn retry) AL. ae een LL ty ent Gey 
ett etd ln A Mal LIT IP rare TA E E at Le a) *es" AA ere Mh SMT ae T be Ub heheh des E Y PTT) 
ELLO vd hehe PU hale, Sele O T E TTTS LT [DET i A 7*9 94 aho. : tas. a o TC bitke ea a, a Free 
De a id LT PET E L TE PIE a LA ZEOT pen ure rr T % . 14 ETTR e "i 
NH ES I ve LIT TM ball al ide “ate D S" ao THÉ regis t a ondak 
ben ed ede ET ST eT Pe 9 99" * , eo ala LITT U TIP Tee sag 5 
T $96 (*9* s^ 8: Put, Ye a got TIT 










Ce tt) MUT TS Medios ME EUR ah PE 
KT} 
SAS 
O AS 
AS SE 































* f 
A DLL ALINE UE A a ee hee rt LI EL el a Dee MILES 
" Lio M oe er d He Y AO ANAIS a il cod TTT A iia dd 
* “urn uL Pak Pie NI IMMER rrt eng Dd t Le au EM 
IC a A ^ PIECE MALLEM CEP S AM TE Fre A eame LLLA TS ET T We Cir) 
A AA LLL E AAA HARO ICO abre e o TA pls Pre e ML ehe 
TA DT] LAIT e... y .. “u... 9.9 s»sag bdo TES MAN UID .. - LE Jeyo Mareo y PITA Pr A al | EIU ec 
OLLEET TELTE] SS eT Ee e. y o tm IL SIR IR tre DI re 
LT e vi ho ABRIR TM SN , -s e dL M TIS EL FETTE no tms pte uy e. 
E "isti SL EL TI T ETA Ms Fd Vb Lot TCR AE TT arie hi. 
SS ca o te eb hen a TEE u 
i ] FEM Cr A) P n HELL "I LII T 
wh bhi’ did To et 2 a I; uni LL TT REI 1 LA IIT. Y E REP = 
A O eee ee > [PTUS dh. A 
taii LER A ake E erm z 
AA DO ELIO e y q LL TT 
ern. heed 






